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Abstract 

One-way measurement based quantum computations (IWQC) may describe 
unitary transformations, via a composition of CPTP maps which are not all uni- 
tary themselves. This motivates the following decision problems: Is it possible to 
determine whether a "quantum-to-quantum" IWQC procedure (having non-trivial 
input and output subsystems) performs a unitary transformation? Is it possible to 
describe precisely how such computations transform quantum states, by transla- 
tion to a quantum circuit of comparable complexity? In this article, we present an 
efficient algorithm for transforming certain families of measurement-based com- 
putations into a reasonable unitary circuit model, in particular without employing 
the principle of defeiTed measurement. 

1 Introduction 

The one-way measurement model of quantum computation (IWQC), presented by |[T] 
and subsequently elaborated in [T T], allows descriptions of unitary transformations in 
terms of operations, such as measurements, which are not themselves unitary. This may 
be regarded as being due to an appropriate combination of the observables which are 
measured on the one hand, and the stabilizer group of the state-space (in the "quantum- 
to-quantum" setting where we admit an input subsystem which may support an arbi- 
trary input state) on the other hand. This intuition may be used in certain cases to prove 
that particular measurement-based computations perform unitary transformations 
using a slight modification of the stabilizer formalism [5J. 

As we are often interested in describing computation in terms of unitary evolutions, 
measurement-based computing presents us with natural decision problems. Is it pos- 
sible to efficiently determine whether a IWQC procedure performs a unitary transfor- 
mation? Is it possible to describe precisely how such computations transform quantum 
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states, by translation to e.g. a unitary circuit of bounded complexity? If not, is it at least 
possible to efficiently verify whether it performs a given transformation? 

These questions arise because simple unitary circuit models — such as uniform 
circuit families generated over the gate-set {H, T, ^Z} consisting of the Hadamard 
gate, the "/s phase-gate, and controUed-Z gate respectively (with no classically con- 
trolled operations and measurement deferred to the end of a computation) — are the 
de facto standard models of quantum computation. Such circuits provide the (thus 
far) best-developed set of idioms for uniformly expressing algorithms independently of 
architecture. Other mathematical representations of quantum computation may be use- 
ful as abstractions of proposals for implementations of quantum computers, as in the 
case of the measurement-based models |6, 1 1 and the adiabatic model |7|; and "purely 
physical" idioms seemingly divorced from uniform circuit families may suggest new 
quantum algorithms, as in the case of the 0{V iV)-time algorithm for evaluating NAND- 
trees by quantum walks |8|. However, simple unitary circuit models are currently the 
preferred means for expressing quantum computations. 

Because of the present uncertainty as to which implementation proposal will prove 
successful, the unitary circuit model represents a natural candidate for an intermediate 
language for representing algorithms independently of the target architecture, even if 
alternative models suggest new algorithms. Therefore, one may argue that how effec- 
tively one can translate from a given model of computation to unitary circuit models, 
and back, is an important measure of how well-understood the model is. Unless some 
one proposal for implementation comes to dominate over the others, or an unlikely 
breakthrough (e.g. the discovery of an efficient classical algorithm for simulating quan- 
tum circuits) is made, translation to a "simple" model of quantum computation such as 
the unitary circuit model is the best that can be reasonably hoped for as a uniform way 
of obtaining low-level semantics of a quantum computation. 

We are generally not interested in arbitrary translations to and from unitary circuit 
models: we should require a translation scheme which does not substantially increase 
the complexity of the operation being translated. For instance, provided a IWQC pro- 
cedure *p for performing a unitary transformation, we may decompose each classically 
controlled measurement into a classically controlled change-of-basis rotation followed 
by a standard basis measurement; applying the principle of delayed measurement, we 
may then replace the classically controlled rotations with coherently controlled ones. 
We may then decompose the gates of this circuit with respect to a desired gate-set, such 
as {H, T, ^Z}. However, in constructions for IWQC procedures such as those defined 
in |l2l|3l, each single-qubit transformation requires at least one additional qubit, one 
two-qubit entangling operation, and a classically controlled measurement which may 
depend on the results from several measurements. If one applies the principle of de- 
ferred measurement to such a IWQC procedure Cp, which "represents" a unitary circuit 
£ using one of these constructions, the result is a circuit which requires many more 
qubits than €, and has at least one additional two-qubit operation for each single-qubit 
operation in £. This represents a constant factor inflation in both the number of opera- 
tions of the circuit; and depending on the ratio of single-qubit to either the number of 
qubits or the number of two-qubit operations in the original circuit £, the asymptotic 
increase in the either the memory required or the number of two-qubit gates in £' over 
£ is in principle unbounded. 
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In contrast, we may consider translation procedures S from the IWQC model, 
possibly with a restricted domain, which yield circuits of "reasonable" complexity. 
In this case, we may establish a reference for the complexity of the resulting circuits 
by considering (as in the previous paragraph) IWQC procedures which arise from a 
translation routine TZ from unitary circuits to the IWQC model. We may then require 
that img(7?.) C dom(5), and that S oTZ maps each circuit in dom(7?.) to one which is 
no more complex (and which performs the same unitary transformation). In particular, 
we will be interested in the case where dom(5) contains a wide range of operations. In 
this article, we require that dom(7?.) consists of the set of unitary bijections generated 
over {H, T, AZ} (without measurements or classical control). Then dom(5) must be 
similarly general, as it must "recognize" different representations of arbitrary unitary 
bijections in the IWQC model, and produce unitary circuits which are equivalent both 
in function and efficiency to the "original". We encapsulate these requirements as 
follows 

Definition 1.1. Let Ce be a model of unitary circuits generated over a gate-set B. 
A semantic map for the IWQC model in terms of unitary circuits is an efficiently 
computable partial function S from IWQC procedures to Cb, for which membership 
in dom(5) is efficiently decideable, and for which there exists a corresponding efficient 
translation routine TZ from Cb back to the IWQC model such that 

• dom(7?.) = Cb and img(7?.) C dom(5); and 

• for any circuit C E Cb, C = {S o 7V){C) performs the same transformation as 
C, and has complexity bounded above by that of C: in particular, 

- C acts on at most as many qubits as C; 

- C has at most the same depth as C; and 

- the number of one-qubit and two-qubit gates in C are each at most the 
number of such gates in C. 

We say that 7?. is a representation map of Cg in the IWQC model, for which S provides 
the semantics. 

Summary of results. We present a single semantic map iS with respect to two re- 
lated constructions for IWQC procedures: one based on the work of Raussendorf, 
Browne, and Briegel |2| (which we will refer to as the simplified RBB construction) 
for producing cluster-state based procedures; and the simplified procedure defined by 
Danos, Kashefi, and Panangaden 1 3 1 (which we will refer to as the DKP construction) 
for producing IWQC procedures whose entanglement graphs are not required to be 
subgraphs of rectangular grids. We will restrict our attention to the case of unitary bi- 
jections, i.e. unitary circuits generated over {H, T, ^Z} which do not introduce fresh 
qubits. In doing so, we present the following: 

''The following definition is a refinement of the conditions described in pages 81-82 of |9| to include 
pairs of maps S and TZ for which, for some IWQC procedures ?p, the procedure {Ti o S){^) differs in 
complexity from 5p. 
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A simple alternative notation for quantum circuits. 

In order to facilitate descriptions of the construction of IWQC procedures, as 
well as the task of determining when two circuits are congruent, we define in 



Section 2.1.3 a variant notation for unitary circuits (and products of linear oper- 



ators generally) having convenient combinatorial properties. 

An explicit account of two constructions of IWQC procedures. 

In order to explicitly describe a semantic map, we present in [Section 23] the 
construction of ifTOI ("the DKP construction") and a simplified version of the 
construction of [jS] ("the simplified RBB construction) for IWQC procedures, 
including an explicit account of the time complexity of these constructions. 

A succinct characterization of measurement dependencies in these constructions. 

Using adjacency matrices of directed graphs, we characterize the classical mea- 
surement dependencies of the IWQC procedures arising from the DKP construc- 



tion and the simplified RBB construction in Section 3.1 



Circuit constructions from combinatorial decompositions. 

In [Section 33] we illustrate how we may efficiently decompose the structures 
underlying a IWQC procedure obtained from the DKP or simplified RBB con- 



structions. Using these decompositions, we show in Section 3.4 how such a 
decomposition may be used to describe unitary circuits which may perform the 
same operation as the original IWQC procedure. 

These results allow us to define an efficient algorithm for computing a semantic map 
S which produces circuits for unitary bijections from an efficiently identifiable class of 
IWQC procedures. 

Related work. 

A study into interpreting the effect of IWQC procedures in terms of unitary circuits 
is performed by Childs, Leung, and Nielsen fTHl via single-qubit "teleportation-based 
quantum computation". Their results may be described as providing alternative proofs 
of validity of the constructions of Raussendorf, Browne, and Briegel ||2] of particular 
unitary operations in the IWQC model, but do not present conditions for when such a 
reduction to the circuit model may be performed. 

The first algorithm to convert IWQC procedures to unitary circuits was presented 
by Kashefi and Danos fT2l, who investigated conditions under which IWQC proce- 
dures produced by the DKP construction perform unitary transformations. These con- 
ditions are formulated in terms of combinatorial structures, which were shown to be 
efficiently detectable in lfT3lfT4ll . Using these structures, fVS\ outlines a procedure to 
generate a unitary circuit using those structures. Building on the results of 1 13 , 15 .fT4||. 
we extend those results to include procedures produced by the simplified RBB con- 
struction. We also characterize the measurement dependencies which arise in such 
IWQC procedures when they are in a particular canonical form, in order to verify 
whether such a IWQC procedure performs the unitary described by the constructed 
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circuit; and demonstrate that the complexity of the circuits which are constructed are 
no more than the complexity of the original circuits. 

The results here represent an extension of results in Chapter 3 of |9| concerning the 
DKP construction, by proving remarks made briefly there concerning the simplified 
RBB construction. In the few places where the terminology here differs from |9 | (or 
from prior results), the differences are explicitly noted. 



2 Preliminaries 

We begin by presenting conventions and notational devices used throughout the article. 
Among these are a concise operational notation for procedures in the IWQC model, 
and a transparent notation for unitary circuits which has convenient combinatorial prop- 
erties. These will prove instrumental to simplify the description of the constructions 
for IWQC procedures in the following section. 

We take for granted a number of basic concepts in graph theory; the interested 
reader may refer to Diestel's text Iil6j for elementary definitions and proofs. 



2.1 Unitary circuits 

We will consider a particular model of quantum circuits: circuits without classical 
control, and with measurement (and trace-out operations) deferred to the end. We refer 
to such circuits as unitary circuits. As is common practise, we will ignore the terminal 
measurements and trace-out operations, and consider only the component involving 
unitary transformations and unitary embeddings. 



2.1.1 Sets of approximately universal gates 

As we noted in the introduction, we must consider a particular unitary circuit model 
with a fixed set of unitary gates. We consider a model closely related to most common 
gate sets; the set {iJ, T, ^Z} consisting of the Hadamard gate H, the ^/s gate T, and 
the controlled-Z gate ^Z. These are expressed in matrix form by 
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V2 



T ■ 
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10 
10 
-1 



(2.1) 



respectively. As H and T generate (up to scalar factors) a dense subgroup of the unitary 
group U(2), and as CNOT = {1 ® H) r\Z (1 (g) H), the set of gates {iJ, T, f\Z} is 
approximately universal for quantum computation. 

To describe computations in the IWQC model, we will also be interested in two 
related sets of unitary gates, involving operators such as J{6) and Zz given by 



J{e) = He-'^o/^ = ^ 

v2 



and 



(2.2a) 
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(2.2b) 



Note in particular that J(^) ^ HT"\ and Zz = e*'^/" (T^ (g) T^) hZ. Therefore, the 
gate-sets 



Sdkp = and Brbb = {JW,^^,^,^ (2.3) 

are also approximately universal for quantum computation]^ 

In unitary circuit models, one usually permits the introduction of fresh qubits into 
the circuit. In order to prevent the preparation of computationally expensive initial 
states, we will require that all fresh qubits be initially in some fixed pure state indepen- 
dent of the other qubits. Without loss of generality, to obtain a simpler correspondence 
to the operations of the 1 WQC model, we will suppose that all fresh qubits are initially 
prepared in the |-|-) oc ( |0) + |1) ) state. 



2.1.2 Parsimonious sets of gates 

As well as being approximately universal for quantum computation, the gate sets de- 
scribed above have another useful property: 

Definition 2.1. A set S of unitary operations is parsimonious if every multi-qubit gate 
in S is diagonal, and if distinct gates in S acting on a common qubit commute if and 
only if they are diagonal]^ 

Each of the gate-sets described above are parsimonious. For instance, it is imme- 
diately clear that {i?, T, A^^} is parsimonious; that Sdkp and Srbb are parsimonious 
follows from the fact that 

^ ^iZa/2^iX(P-a)/2^-~iZl3/2 ^ (2.4) 

which is equal to 1 2 if and only \f a — (3. 

The motivation for this terminology is that, for a "parsimonious" gate-set S and for 
any orthogonal basis {|0o) , If/"!)} G 'H2 other than the standard basis, there is at most 

''The gates J(7r/4) and J(— 7r/4) alone generate a dense subgroup of U(2) up to scalar factors: this is 
in fact shown by 1 17 1 (Section 4.5.2) in the proof that H and T are approximately universal for single-qubit 
operations. We include gates J {9) for each d an integer multiple of 7r/4 for both Bdkp and Srbb in order 
to obtain a closer correspondence with circuits generated over {H, T, AZ} in the construction of IWQC 
procedures in later sections. 

■^The title of |3 1 predates this definition of "parsimonious" by four years, and does not refer to the concept 
we have defined here. However, the sets of unitary operations {J{S), ^Z}g^^ and {J{S), ^Z}g^wz 
described there are both parsimonious in this sense. 
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one gate in S for which the vectors {(pj) are eigenvectors. Requiring that a get-set be 
parsimonious is a strong constraint: by definition, operations in a unitary circuit over 
such a set commute if and only if they are of the same "type" (they are either equal or 
inverse to each other), if they act on disjoint sets of qubits, or if they are both diagonal. 
This is evidently useful in determining when two quantum circuits are equivalent up to 
transposition of gates. 

2.1.3 Stable-index tensor notation for unitary circuits 

For a representation map TZ from unitary circuits to IWQC procedures and a semantic 
map S, to show that S oTZ maps unitary circuits to equivalent circuits of the same (or 
lesser) complexity, we may attempt to show that SoTZ maps each circuit in dom(7?.) to 
one which is congruenj^to the original (possibly up to a set of "trivial" simplifications). 
For this purpose, it will be useful to have a notation for unitary circuits which is easy 
to produce by an algorithm, in which irrelevant details such as the order of commuting 
gates are essentially absent. In general, this is a problematic task: however, we will be 
satisfied with a notation which solves this problem for parsimonious gate-sets, such as 
the sets of gates which we use in this article. 

One solution is to describe a notation whose meaning is invariant under permutation 
of any of the terms. An example of such a notation is Einstein summation notation 
( ifTSl , Section B5), in which (for unitary circuits on qubits) operations U are described 
as tensors, using their coefficients with respect to the standard basis on individual qubits 
involved in the circuit. We use superscripts for column-indices, and subscripts for row- 
indices: thus, we have 

Hi = j=,i-^y\ (2.5a) 
JWt = ^(-l)'^''e*(2f'-i)«/2^ (2.5b) 

V 2 

AZl;: = S,,t6sA-'^y': (2.5c) 

and so on (where 6j^k is the Kronecker delta). In an expression involving multiple 
gates, two such tensors are identified as acting on a common qubit when the column- 
index of one agrees with the row-index of another: the coefficients of the compound 
tensor are computed by summing over the range of the repeated index. The order 
of the multiplication of the tensors is then determined by the location of the index: 
an operation U strictly precedes another operation V if one of the row-indices of U 
matches a column-index of V. To represent tensor products of two operators, it suffices 
to juxtapose them with disjoint sets of indices. For example, given 

C = {l(g>W) {V ^1^1)U {l(g)T^l) {l^\(f>)(g)\ij)) (2.6a) 

for a generic three-qubit operator U, two-qubit W, single-qubit operators V and T, and 

''For two unitary circuits, or two IWQC computations, we will say that tlie computations are congruent 
if they differ only by transpositions of commuting operations (and by a possible relabelling of the qubits 
operated upon). 
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state- vectors \4>) , \tp) G 7^2, we may equivalently write 



lJlb3,C2ira2TTai,b2,CirTibi \ J\bo 

^^b2.ci ^ai '^aoM,co bo IW 



(2.6b) 



Furthermore, the coefficients described on the right-hand side remain the same under an 
arbitrary permutation of the factors, and the equahty remains true after an application 
of any (one-to-one) relabelling of the indices. 

However, while the order of the factors is not significant in Einstein notation, cir- 
cuits which are equivalent up to rearrangements of commuting operations still give rise 
to tensor expressions which are not congruent via permutations of factors and index re- 
labelling. This is because the tensor indices themselves fix an order of multiplication: 
for instance, although the operations ^Za.b, ^Za^c, and AZ^ c commute (where the sub- 
scripts here indicate that the operations are performed on some pair chosen from three 
qubits named a, b, and c), the products 



^Zh r ^Zar ^Zai 



/\Za_r ^Zh r AZ 



AZt'''"^ AZ"-'^-"^ hZ°-^-\} 
/\Z<^2,C2 AZl"^'"' \Z°-^-\} 



and 



(2.7a) 
(2.7b) 



give rise to expressions in Einstein notation which cannot be made the same by a rela- 
beling of indices and re-ordering of terms. 

In the case that two gates commute if and only if they are both diagonal, as in 
circuits generated from parsimonious sets of gates, we may resolve this problem as 
follows. We use a notation similar to Einstein notation, but in which in which an index 
can be repeated multiple times without summation, when the index represents a qubit 
in a part of a unitary circuit in which the standard basis is preserved for that qubit|^ 

Stable index notation. Consider a unitary operator U on N qubits, whose eigenvec- 
tors are all of the form |si)(g) • • • ^ |s„)(8) |^')for si, . . . , s„ e {0, 1} and \^) € TYf", 
where N = n + m. We may represent the non-zero matrix coefficients of U in the 
standard basis by 



U 



Vn+i,--- ,yN 



■^71 + 1 ) 



U, 



Sl, . . . , s„, ?/„+!, .. .,yN\U |si, . 
si,--.,s„,yn+i,...,yN 



(2.8a) 



(By hypothesis, the coefficients of U are zero when the values of the first n row- and 
column-indices differ.) That is, for tuples s e {0, 1}" and a, d e {0, 1}'", we define 
U[s by the operator equality 



U 



se{o,i}" 

a,de{0,l}" 



i)(d|) 



(2.8b) 



■^If one interprets tensor indices as a random variables ranging over {0, 1}, whose value may be realized 
by a standard basis measurement, this con'esponds to re-using indices when the value of that variable would 
be unchanged (or "stable") under some number of subsequent operations. 
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Terminology and composition conventions for stable-index notation. In an ex- 
pression U[s g] as defined by ( |2.8b| l, we call the indices Sj stable, the indices dj dep- 
recated, and the indices aj advanced. A product of such expressions is well-formed if 
each index is advanced at most once and deprecated at most once; there is no limit to 
the number of times an index may occur as a stable index. An index which is advanced 
in one term and deprecated in another is a bound index; otherwise it is free. We sum 
implicitly over all bound indices in a stable-index tensor expression, as in Einstein no- 
tation: we then form composite tensors V o U hy matching advanced or stable indices 
of U to corresponding deprecated or stable indices of V. (If the indices of U and V 
are disjoint, this is a tensor product of U and V.) 

Example 1. The operation T is a single-qubit diagonal operator, and so has one 
stable index and no deprecated or advanced indices; the Hadamard operation H is 
a single-qubit non-diagonal operator, and so has no stable indices and one depre- 
cated/advanced index each. We may then write their coefficients as 

T[a.] =e(2^-i)"/8 and = ^(-1)™. (2.9a) 

The J('V4) operation may be defined by composing these two operators, while the 
inverse of J(~V4) is the other composition of these operators; we may write these 
compositions as 

J(^/4)[|] = H[l\ T[x] and J{-^U)\l\ - T\v] H[1\ . (2.9b) 

Finally, the tensor product P = T ® H has eigenvectors of the form \z) ® \^) for 
z G {0, 1} and |?/') G 'H2 an eigenvector of H: its coefficients may then be written as 

Identified systems. For operators U : Hf^ — > Tif^ acting on distinguishable 
qubits, we describe how U acts on these qubits by fixing an order of the qubits, where 
the advanced indices and deprecated indices of the qubits are presented in the same 
order. 

Extension to general linear operators. For operators A : Hf^ — > Hf^ where 
N and M may differ, a stable-index tensor expression for A will have advanced indices 
which do not have corresponding deprecated indices if AI > N, or deprecated indices 
which do not have corresponding advanced indices if > M. Then A describes an 
operation which either adds qubits to the input space, or discards qubits from the input 
space. (These may occur when A describes an isometric embedding, or a single-qubit 
measurement with a selected outcome, respectively.) We may represent a qubit which 
is added by an advanced index with no corresponding deprecated index (padding the 
sequence of deprecated indices by a space, dot, or similar placeholder), and similarly 
for deprecated indices representing qubits which are discarded. 

Example 2. We may transcribe the operator expression given in ( |2.6b[ ) as 

I0>[^] ■ (2.10) 



C 



a2,fa3,C2 



w 



63, C2 



V 



u 



ao,6i,co 



T 
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The order of the indices is purely conventional (in the same way that the order of the 
indices of C in \2.6h) is purely conventional); we could equally well define C above as 
1 ; the change in the order of the qubits causes a corresponding change in the 



matrix representation ofC. The adjoint of the tensor in \2.\Q\ is given by 



<12 ,63 ,C2 - 



A formal procedure for translating tensor notations such as in ( |2.6a| i, or involving 
operators Ua,b.--- .q labelled with the qubits that they act on, is easy to define: such a 



procedure is presented in Section A. 2. 1 



Example 3. Considering the same tensor C as in ( |2.10[ ), in the special case that U has 
eigenvectors of the form \a) \b) jV') for a, 6 G {0, 1} and for some states G Ti.2, we 
may simplify the expression ofC to 



C 





= W 




V 


02 


U 


aoM^ 


T 


"fcl" 


a-o, ■ ■ ■ 


bi.ci 


ao 


Co 




bo 



'bo] 



(2.11) 



The operators described by the expressions in ( |2.10| l and ( |2.1 1[ ) are identical in this 



The following example illustrates the intuitive purpose of stable-index notation: by 
allowing stable indices to be repeated multiple times while remaining "free" variables, 
we also discard some redundant combinatorial information about the order of the prod- 
uct of operators. 



Example 4. For generic unitary operators U G U(2) 



(gin 



AC/ 



, we may define the operator 
U. (2.12a) 



Let V G U(2) and W G U(2) For qubits a, b, c, and d, consider the circuits 
^Va^d ^Wa^b,c ond ^Wa.b,c ^Va^d (where the subscripts here denote the qubits on which 
the operators act). We may represent these two circuits with stable-index expressions 



AV 



di 


AW 


^buci 


do 




bo, Co 



and 



AW 



bi.ci 
bo, Co 



AV 



di 
do 



(2.12b) 



respectively. Note that these are equivalent expressions, up to permutations of the 
factors: this may be regarded as being due to the fact that the operations /Wa.d <^nd 
/^Wa,b,c commute (i.e. the circuits are congruent). 

In contrast to Einstein notation, while the numerical meaning of a stable-index 
tensor expression is invariant under permutations of the terms, information about the 
order of consecutive operations is not represented when these operations preserve the 
standard basis of qubits on which they both act. Circuits which are congruent up to the 
re-ordering of commuting gates may then give rise to "equivalent" stable-index tensor 
expressions in the following sense: 

Definition 2.2. An isomorphism of stable-index tensor expressions Si and 5*2 is a 
bijective relabelling A of the index labels of one stable-index expression Si to those 
of 5*2, and a bijective mapping r of terms of Si to those of S2, such that for any term 
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U[si,...,s„ ] in iS*! (with n ^ stable indices, and m ^ each of deprecated 

and advanced indices including placeholders), there is a term in 5*2 given by 



r U 



ai, ■ 



di, 



U 



(2.13) 



(where placeholders on the left are mapped by A to placeholders on the right) |^ 

Isomorphic stable-index expressions represent not only equivalent unitary opera- 
tions (which follows from the implicit summation convention), but congruent circuit 
decompositions as well: 

Theorem 2.1. Let Ci and C2 be two sequences of unitary operators over a common 
gate set, and acting on the same number of qubits. If the stable-index tensor represen- 
tations of Ci and C2 are isomorphic, then C'l and C2 are congruent via a re-ordering 
of commuting operators and relabelling of the qubits. Furthermore, if Ci and C2 are 
generated from a parsimonious set of gates, the converse also holds. 

Thus, for unitary circuits constructed from a parsimonious set of gates, isomor- 
phism of stable-index tensor expressions is equivalent to congruence up to rearrange- 
ments of commuting gates. This feature is the motivation both for this tensor notation. 



and for the interest in parsimonious sets of gates as in Definition 2.1 We only explic- 
itly use advantage of this feature of this tensor notation in Section 3.5 however, other 
combinatorial properties of stable-index notation will frequently prove convenient in 
describing constructions of IWQC procedures from the circuit model. The proof of 



Theorem 2.1 is provided in Section A. 3 



While stable-index expressions lack any information about the order of operations 
which preserve the standard basis of the qubits that they act on in common, it is possible 
to retrieve information about the order of non-commuting gates, in order (for example) 
to determine how to actually perform such an operation as a physical transformation. 
As in Einstein notation, the order of non-commuting gates can still be determined from 
the indexing. 

Definition 2.3. For a stable-index tensor expression S, the index successor function f 
for S is the mapping defined on the indices S of S which are deprecated but which have 
a corresponding advanced index a (corresponding conventionally to the same qubit), 
such f{6) — a for each such S and a. 

Note that for a deprecated index v € dom(/), any term U which involves the index 
f{v) must be performed after any operation in the circuit acting on v (excepting U 
itself, if it is the unique term in which v and f{v) both occur). We may define: 

Definition 2.4. For a stable-index tensor expression S, the interaction hyper-graph for 
iS* is the hypergraph G on the set of indices in S, whose edges consist of those sets of 
indices which are acted on by some operation in C. 



'Note that isomorphisms of stable-index representations may change the order of the terms as well as 
the index labels. That is, they are combinatorial homomorphisms, not just a relabelling of the indices in the 
terms of Si kept in the same sequence. 
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For any v S dom(/), the index f{v) is adjacent to v in this hypergraph. We may 
then use / and G to estabhsh a partial order u ^ v on the indices of S, which describes 
when an index u is deprecated before another index v, given an ordering of the terms 
which is consistent with an order in which the operations may be performed. We may 
define ^ as follows. For two indices f{v), /{w) which occur in a common term, each 
must be deprecated strictly after the indices v and w in order for the term to describe a 
well-defined (acyclic) product of operators. If v ^ w denotes the adjacency relation in 
G, we then have 

V ^ f{v) and w ^ f{v) v ^ w (2.14) 

for any v G dom(/). We can then describe a pre-orde:[^ C/ =^ on the operations in 
S, where U ^ V if U has an index u and V has an index v such that u ^ v. We then 
have U ^ V ^ U if and only if either U and V are the same term, or different terms 
which act on common qubits but commute; if U ^ V ^ U, then U must occur strictly 
earlier than V. We may then obtain an order for the operations in S by finding a linear 
order which extends . 



In the case of the gate-sets defined in Section 2.1.1 the operations all act on single 



qubits or pairs of qubits. The interaction hyper-graph G of (a stable-index expression 
for) a circuit over such a gate-set will then be a conventional graph. If we ignore single- 
qubit diagonal operations which add no information to the ordering of the indices, and 
note that the square of the two-qubit operations are in each case a product of local 
unitaries, we may even describe the circuit by a simple graph (without repeated edges 
or loops). This graph structure, and the ordering structures we have described above 
for circuits will prove convenient for analyzing constructions of IWQC procedures. 

2.2 The measurement calculus 

The way in which computation is performed in the IWQC model is very different from 
unitary circuits as we describe them in [Section 2A In the latter case, computation 



is performed with interleaving non-commuting single and two-qubit unitaries without 
any classical control or intermediate measurement. In contrast, a IWQC computation 
consists of an initial stage in which all of the (mutually commuting) two-qubit oper- 
ations are performed, followed by a sequence of classically controlled measurements 
and (for procedures which produce a residual quantum state as output) classically con- 
trolled single-qubit unitaries. In this section, we briefly review the elementary oper- 
ations in the IWQC model and describe standard forms for procedures in the IWQC 
model as a preliminary to describing IWQC constructions for unitary transformations, 
essentially following the terminology and notation introduced by Danos, Kashefi, and 
Panangaden lilOJ . 



8A pre-order =>! is a binary relation whicli is reflexive {x =^ x for all x) and transitive (if x =4 V ™d y ^ z 
then X =^ z), but with no further requirements. Partial orders and equivalence relations are both examples 
of pre-orders. An example of a pre-order which is neither an equivalence relation nor a partial order is the 
relation x =4 V '^==^ Re(x) ^ Re(y) for complex numbers x and y. 
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2.2.1 Conventions for writing CPTP maps 

In the IWQC model, classical information is obtained via measurements which may 
in practise be destructive (e.g. as with photodetectors). Thus, it is usually preferable 
to describe IWQC procedures by completely positive trace-preserving (CPTP) maps 
acting on density operators, where the input and output spaces of these operators may 
not be supported by the same set of qubits. 

In order to better distinguish CPTP maps from linear operators on Hf^, we will 
represent CPTP maps in either sans-serif (as in J ) or fraktur (as in ^ ) typeface. We 
also adopt the following additional notational devices: 

Measurement results. 

When a qubit v has been measured, we will represent the result by a bit s[v\ G 
{0,1} whose value is only defined after the measurement. 

Operands of CPTP maps. 

The qubits on which an operation Op^ ^ acts, whether to introduce it to the sys- 
tem by preparation, discard it by measurements, or to transform it unitarily, are 
represented in the subscript. (The bit-registers s[v] storing measurement results 
are typically omitted from such lists of operands.) 

Parameterization of CPTP maps. 

Operations Qp"'^'^'"''^'™' ■ ■ in a IWQC procedure will often be determined by 
fixed parameters a, /3, . . . and the results slv], s[w], ... of measurements on qubits 
v,w, . . .. (Typically, fixed parameters a denote the values of angles in the range 
(— TT, tt], and operations will depend on the parities s[v] + s[w] + ■ ■ ■ of measure- 
ment results.) Such parameters are always described in the superscript. 

Finally; as there is little risk of confusion, we will typically abbreviate a composite 
CPTP map Op(i) o Qp'^) o • • • o Op(") by Gp^^) Op^^' • • • Op'") . 

2.2.2 Elementary operations of the IWQC model 

As in unitary circuit models, we define IWQC computations as compositions of ele- 
mentary operations. We define these following ||T0|. 

The initial stage of a IWQC computation is the preparation of a large entangled 
state, which may be obtained by performing entangling operations on a collection of 
qubits prepared in the |+) state. Furthermore, the entangling operation which is per- 
formed may be described as a product of operations acting on distinct pairs of 
qubits. It will then be useful to define the operations 

K{ps) ^ Ps<E)\+){+\v and E^^{ps) ^ ^Z^n, ps ^Z„w (2.15) 

for a state ps on a set of qubits 5*, where the subscripts on the right-hand sides indicate 
the qubits on which the operators act. For the map N„ , we require that v ^ S, and 
obtain the system S U {v} as output; for , we require that v,w G S, and obtain 
the same system S as output. We refer to as a preparation map, and Ey^ as an 
entangling map. 
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The single-qubit measurements we require will be with respect to orthonormal 
bases on "the XY plane" (of the Bloch sphere): that is, with respect to uniform su- 
perpositions of |0) and |1) which we denote by 

\+o) = j=^{\0) + o''\l)) and I-,) = J=(|0)-c'^|l)) (2.16) 

for angles — tt < 9 ^ tt. Using these states, we may define the measurement operators 

MliPs) = {+o\ PS \+o)v ® |0)(0|3[^,, + (-,1 PS \-e)v ® , (2.17) 

acting on a density operator ps on a set of qubits S. (Recall that s[v] is a classical bit 
storing the boolean measurement result]^ We require that v G S; the output consists 
of the system S \ {v}, ignoring the newly defined bit s[v]. 

We extend this notation in a simple way to describe how measurements are adapted 
according to previous measurements: we introduce the abbreviation 

M^^'^ = M^-''>'-' (2.18) 

for a classically controlled measurement performing either M" or M depending on 
a boolean expression (3 = s[a] + s[6] + • • • representing the parity of some set of mea- 
surement resultsjl] We refer to 9 as the default (measurement) angle of the operation, 
and (3 a sign-dependency of v; we say that v has a sign-dependency on w when such an 
expression j3 depends non-trivially on s[w\. 

For IWQC procedures which produce a residual quantum state as output (as we 
consider in this article), it is also necessary to adapt the final state depending on the 
outcomes of the measurements by using classically controlled single-qubit operations. 
For boolean expressions /3 describing the parities of some measurement outcomes, we 
may decompose these corrections into X and Z operations with the operators 

Xf(ps) = p5Xf and Z^^ips) = PS , (2.19) 

acting on a density operator ps on a set of qubits S which contains v. 

Finally, a purely classical operation which is useful in descriptions of IWQC pro- 
cedures is a classical bit-flip operation on a measurement result s[?;]: 

SfM(p) = I/3)(0LhP |0)(/^lsH + |l-/?)(l|sM^ |l)(l-/3|sM ' (2.20) 

which acts on the bit s[d] which stores the result of a previous measurement^ This 
operation is also known as a signal shift operation, as it can be used to remove the 

''While we typically do not represent the classical bits used in operations in terms of density operators, 
we may do so when this is convenient for exposition. 

'The calculus defined in 110] also admits the possibility of a measurement angle being modified by 
a possible addition of tt, depending on a boolean expression (as with sign dependencies above). As this 
elaboration is not necessary for measurement-based simulation of unitary circuits, we omit it in order to 
present a streamlined description of how to construct a IWQC procedure in "standard" form. 

J Such "shift" operators are usually not noted in accounts of measurement-based computation (with (lOJ 
being an exception), probably because they are simple operations on classical bits. However, in any imple- 
mentation of quantum computers based on single-qubit measurements, shift operators may in practise be 
used to describe how parity expressions /3 are computed in the midst of computations. Thus, while much 
of the analysis of this article follows the convention of eliminating shift operations where possible, it seems 
productive to extend the sense of "standard forms" found in the literature to include shift operations between 
measurements. 
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dependency on the result (or signal) of the measurement of some qubit a from the 
measurement of a later qubit v; the influence of the measurement on a is then accu- 
mulated into the classical bit s[v]. As a result, shift operators Sfj„j only occur after the 
measurement of the corresponding qubit v. 



2.2.3 Elementary transformations of IWQC procedures 

A typical IWQC procedure is a composition of many operations as described above. 
In composing them, we may simplify them through simple transformations of the oper- 
ations involved. At the same time, we wish to bring these expressions into a canonical 
form in which the entangling and measuring operations are separated into different 
phases of the computation (one of the motivations of the IWQC model being to sep- 
arate the tasks of creating entanglement, and maintaining control of the states of indi- 
vidual qubits). We now describe these transformations. 

A well -formed IWQC procedure prepares each qubit at most once, in which case 
that preparation is the first operation performed. We may then accumulate all prepara- 
tion operations N^, to the beginning of a measurement-based procedure]^ Thus, for an 
operation Opg acting on a set of qubits S, we have 

NyOpg = OpgNy , whenv ^ S. (2.21) 

(When V E S, the expression on the left-hand side is not well-defined.) We may 
collect these preparation maps together and describe them as a single operation on 
many qubits: to this effect, we define the shorthand 

N„„ ••• N,,N,, = N{„^,...,„„} . (2.22) 

Similarly, because every operation aside from the preparation maps and entangling 
maps are either measurements (which are the last operation performed on any qubit), 
Pauli corrections (which we may easily transform by commutation with AZ, a Clifford 
operation), or shift operations (which do not act on qubits), we may apply the following 
commutation operations to accumulate all of the entangling maps to the beginning of a 
IWQC procedure (just after the preparation maps): 

M^''3 = M^''' , when p i {v, w}; (2.23a) 

E«u, Sf[j,j = Sf[p] E„„ , whenp ^ {v,w}; (2.23b) 

^vw = Z£ E^u, , and similarly for E^.^, ; (2.23c) 

E„„ = , and similarly for E^^, Zf, . (2.23d) 

(When p G {v, w}, the expression on the left-hand side of p.23a| i is not well-defined, 
and the bit s[p] on the left-hand side of ( |2.23b| i does not have a defined value.) As 
is a diagonal operation, we also have 

^uv — E^^ . (2.23e) 

''Composition is performed from right to left, as witli standard notations for unitary circuits; tlius, the 
beginning of a IWQC procedure is on the right-hand side of a sequence of operations. 
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As with preparation maps, we may collect these entangling maps together and de- 
scribe them as a single operation on many qubits. Controlled- Z o perators act symmet- 
rically on pairs of qubits: in stable-index notation (Section 2.1.3 1, we have ,iu] = 
AZ[iu,i']. Then the order of the entangling operations and the order of the qubits acted 
on for any single entangling map is unimportant. As l\Z is self-inverse, we have 
Eut) E„„ = 1; then we may represent a product of entangUng maps simply by an 
entanglement graph 



-G = 



uveE{G) 



(2.24) 



where the edge-set E{G) consists of all of the pairs of qubits on which entangling maps 
act. 

We may also consider how to commute correction operations past other operations. 
First, note that commuting and Z'^ operations on a state vector accrues at most a 
sign factor. Then, the corresponding CPTP maps commute: 



X/3Z7 ^ 



(2.25) 



Similarly, X operations on any two qubits commute, as do Z operations. We may 
accumulate corrections of the same type on a common qubit: 



ji 77' — 77+7' 



(2.26) 



With respect to measurement operations, corrections on a qubit v may be interpreted as 
(classically controlled) changes of basis for the measurement on v; then, we can con- 
sider "absorbing" correction operations on qubits to be measured into the measurement 
operations. Consider a measurement in the \+g) , \ —g) basis for some — tt < 9 ^ tt. 
We have 



X\ 



+0)= ^(c^''|0)+|l)) (X -l=(|0) + e-''|l)) = !+(_,)) 



(2.27) 



and similarly X |— 0) cx | — (-e))- We may then incorporate corrections into the basis 
of a subsequent measurement on v via the equation 

M^X(^ = MH)"'' = M^^'^ (2.28) 

for a boolean expression /3, using the notation defined in ( |2.18| l: more generally, we 
have M^'/^t; = M^'''+''. We may describe the effect of Z corrections on measurements 
in a similar fashion: 



Z\ 



1 

V2 



(2.29) 



and similarly Z \ 
we then have 



For an arbitrary state ^5 on a set of qubits S including v, 



ZvPsZl 
ZvPsZl 



IPS 

\PS\ 



and 



(2.30a) 
(2.30b) 



16 



Thus, Z corrections effectively exchange the roles of the states |±6i) for measure- 
ments made in that basis. As this exchange is independent of any classically-controlled 
change of sign in Q, we then obtain 

y\';;^Tl = S^^U'^f" (2.31) 

for arbitrary boolean expressions /3 and 7. (That is, we may toggle the measurement 
result s\v\ instead of performing a Z correction on v prior to measurement.) For such 
a shift operator S^j^j which occurs after a measurement M* and before any operations 
using s,\v\ for classical control, we call 7 a bit-dependency of 

There are some further simplifications to measurement operations which are pos- 
sible when the measurement angle 6* is a multiple of V2. We refer to these as Pauli 
simplifications, as such measurements are with respect to the eigenbases of the X otY 
PauU operators. For e {0, tt}, we have = \+e) and Z \—g) cx |— (-g)); thus, 

we have 

M^tt;/? ^ ^mTT ^ (2.32) 

SO that we may ignore sign-dependencies in this case. For 6 G {""/a, ^M}, we have 
X \+g) cx Z 1+0) and similarly for \~e), so that 

M;^'''^;'^ = Mt^'^X^^ = Mf/^Zf = Sf[,jMf''=. (2.33) 

Thus, measurements in so-called Pauli bases may be performed without dependencies 
on previous measurements. In summary, we have the following relations which allow 
correction operations to be commuted past or absorbed into measurements: 

(2.34a) 
(2.34b) 
(2.34c) 
(2.34d) 
(2.34e) 
(2.34f) 
(2.34g) 

We will occasionally be interested in describing transformations of IWQC proce- 
dures "without Pauli simplifications", in which case we do not use the relations given 
in and p.34g| ). 

Finally, we consider the transformations which are possible by commuting shift 
operators. We may accumulate multiple signal shifts acting on a single measurement 
result: 

c/3 C7 _ c/3+7 (2 3^^ 

'This corresponds essentially to what are called 7r-dependencies in 19J, and to "dependencies caused by 
Z-actions" in [ 10]. 



|\/|e;/3x7 


= M^'^ : 






V u 




when u 




M^'^Z^ 




when u 




V V 








M^'^Z^ 

' 'v V 
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where this sum is evaluated modulo 2. If desired, we may also eliminate shift opera- 
tions entirely from a IWQC procedure, by incorporating the bit-flip which is performed 
on the measurement signal slv] into every expression which depends on s[v]. This has 
the effect of increasing the complexity of the classical control expressions, but when 
performed for all measurement results slv], also makes explicit the dependency of any 
operation on prior measurement results. For an arbitrary operation Op^^ ' which 
is classically controlled by some function / taking s[v] as a parameter, we have 

Op/(-*l.-)sf(„j = Sf(„jOp/(-=M+^--); (2.36) 

that is, we substitute s[?j] + (3 wherever s[?j] occurs, again evaluated modulo 2. We may 
use this commutation relation to accumulate all shift operators at the right-hand side of 
a IWQC procedure, after every operation acting on qubits: we refer to this as signal 
shifting. Conventionally, we may discard any classical measurement results when they 
are no longer required; thus, at the right-hand end of a IWQC procedure, we may 
remove any signal shift operators which remain. 

2.2.4 Standard forms for IWQC procedures 

Using the simplifications above, we may transform IWQC procedures into "standard" 
orderings of the operations which are comparatively convenient for discussion, and are 
also of low operational depth: we now describe two different such forms. 

As we have noted above, in any well-formed IWQC procedure, it is possible to 
commute all preparation and entangling maps to the beginning of the procedure (on 
the right), and to propagate all correction operations to the end (on the left), possibly 
absorbing corrections into measurement operations while doing so. Some of the cor- 
rections will also induce shift operations just after the measurement into which they 
are absorbed; and we may do this in a uniform way without explicit knowledge (or 
despite uncertainties in) the values of the default measurement angles. This motivates 
the following: 

Definition 2.5. A IWQC procedure *p is said to be in standard forn^if it can be 
decomposed as 

*P = CoMoEgoNp, (2.37) 

where Np is a product of preparation maps acting on some set of qubits P, Eq is a 
product of entangling operations acting on pairs uv connected by an edge in the graph 
G, M consists entirely of measurement operations and shift operations (where any shift 
operation S*[„] immediately follows the measurement M*), and C consists entirely of 
correction operations. The procedure of transforming a IWQC procedure into such a 
form using equations ( |2.21| i, ( |2.23| l, and ( |2.34a| l-( |2.34e[ ) is called standardization. 

Note in particular that the process of "standardization" does not apply Pauli sim- 
plifications or signal shifting. For practical applications, it is probably preferable to 
perform Pauli simplifications, at least; but for the purposes of discussion it is useful to 

""This terminology is similar to, but differs slightly from, the usage in 1101 and (9|. 
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consider the case where such simphfications of expressions for IWQC procedures are 
not applied. 

By applying Pauli simplifications on measurements and signal shifting, it is pos- 
sible to have equivalent IWQC procedures which are each in standard form. Rather 
than attempt to characterize when such expressions are equal, it is useful to consider a 
form in which any variation of measurement dependencies has been removed. We may 
define such a form as follows: 

Definition 2.6. A IWQC procedure *p is said to be in normal form if it is in stan- 
dard form, does not contain any shift operations, and any measurement operator with 
an angle 6 E has no measurement dependencies. We call the a procedure *p in 
normal form the normalization of another *P' if we may obtain *p from *P' by standard- 
izing applying Pauli simplifications ( |2.34f| )-(2.34g i, and then performing signal 



shifting ( |2.36[ ), removing the shift operators from the left after doing so. 

While it may not be desirable to compute such expressions for the practical appli- 
cation of IWQC procedures, obtaining the normalizations of IWQC procedures will 
be useful for the analysis of this article. 

2.2.5 Composability of IWQC procedures 

To perform universal computation, it is necessary to compose IWQC computations: 
that is, to apply a second IWQC procedure to the output system produced by an- 
other IWQC procedure Therefore, when composing IWQC procedures *p = 
^2 o we will be interested in restricting the qubits which are allocated in to be 
distinct from any qubits involved in except those which *Pi produces as output. 

It will be convenient at this point to introduce a combinatorial structure which par- 
tially describes a IWQC procedure in order to define composability conditions (as 
well as to visually convey the structures which will be important in later sections). 
For a IWQC procedure *p, we may consider the information contained in *p which is 
independent of ordering or of the particular measurement angles or outcomes: 

Definition 2.7. The geometry underlying a IWQC procedure *p is a triple (G, /, O), 
where G is the entanglement graph on the qubits of *p specified by the entangling maps 
Euuj in *p, / C V{G) is the set of qubits which are not prepared by *p, and O C V{G) 
is the set of qubits which are not measured by *p. We refer to / and O as the input 
subsystem and output subsystem of *p, respectively. (When a given graph G may be 
inferred from context, we may write P = V{G) \ / and 0° = V{G) \ O for the 
complements of these subsystems.) 



We will often illustrate geometries in the style shown in Figure 1 Such diagrams 
give a partial description of a IWQC procedure, in which each qubit u e /"^ is prepared, 
the entangling operations implied by the edges vw G E{G) are performed, and the 
qubits w € O"^ are measured. However, the measurement order, measurement angles 
and dependencies, and any correction/shift operations are left unspecified. 

Using geometries, we can express conditions for the composability of two IWQC 
procedures in terms of a sense in which the geometries of these patterns may be "com- 
posed": 
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o 

Figure 1: Illustration of two geometries for IWQC procedures. Labels for the vertices 
have been omitted. Note in particular that the sets / and O may overlap (indicated by 
the more darkly shaded region in the geometry on the right). 




Figure 2: Illustration of the composition of two generic geometries for IWQC proce- 
dures. Notice that composition in the diagrams is left-to-right, while the composition 
in written notation is right-to-left. 



Definition 2.8. For two geometries Qi — (Gi, /i, Oi) and Q2 = (G2, /2, O2), we say 
that Qi and Q2 are composable if [y(G'i) n y(G2)] C \j2 D Oi] . The composition 
^ = CJ2 o ^1 is then defined by ^ = (G, J, O), where 

V{G) = V{Gi) U V{G2) , (2.38a) 

E{G) = E{Gi)AE{G2), (2.38b) 

/ = h U (/2 \ Oi) , (2.38c) 

O = O2 U (Oi \ /2) , (2.38d) 

where B = [A \ B] U [A \ B] is the symmetric difference of A and B. Two 
IWQC procedures ?Pi and ^2 are composable if their underlying geometries Qi and 
Q2 are composable; the underlying geometry of the procedure *P2 o *Pi is then Q20Q1. 



Figure 2 illustrates the composition of two geometries (and, by way of synecdoche, 
represents the composition of two IWQC procedures). 

In order to obtain a composition of two IWQC procedures ^2 o we may relabel 
the qubits involved in '^2 in order for the conditions of ( |2.38| l to hold. From this point 
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onwards, when we wish to compose two generic IWQC procedures *Pi and ^2, we 
will assume that they are composable in the sense above, eliding over the procedure of 
relabelling the qubits to be prepared which may be necessary. 



2.3 The DKP and simplified RBB constructions 



Having established stable-index tens or notation ([Section 2.13 1 and a calculus of the op 



erations used in a IWQC procedure ( [Section 22 1, we may describe two constructions 



for approximately universal quantum computation in the IWQC model|^The first con- 
struction we will refer to as the DKP construction, and is based on the work of Danos, 
Kashefi, and Panangaden tlOJ for the construction of IWQC procedures in the ab- 
sence of topological constraints on the interactions of qubits. The second construction 
is derived from the elementary constructions among those presented by Raussendorf, 
Browne, and Briegel tZJ for producing IWQC procedures, under the topological con- 
straint that the entanglement graph be an induced subgraph of a suitably large rectan- 
gular grid; we call this the simplified RBB construction^ 

2.3.1 Construction of the gate-sets Bdkp and Brbb 

The first step in both of the constructions we describe is to identify IWQC procedures 
for the gate-sets Sdkp = {>^(^): ^^lee^iz Brbb = {Ji(^), Zz}g^„.^. We shall do 
this in overview, leaving the details for |Section B.l| 

By definition, the CPTP entangling map E„,„ performs controlled-Z operations on 
pairs of qubits: for a system S including qubits v,w,we have 

^vw{ps) = ^Zy^ps^Zy^. (2.39) 

We may implement Zz and J{d) operations in the IWQC model with somewhat more 
elaborate procedures: 

33,,^ = Zf 1 MJ^ E,, N, , (2.40) 
^t/v = X=MM-^E,^N^ , (2.41) 

where the subscript w /v in the latter denotes that it discards v from the input system 
and introduces a new qubit w. It is possible to show that the effects of these two IWQC 
procedures are given by 



^iv.w{Ps) = Zzy^y, PS Zzl^^ , (2.42) 



"The presentation of such construction presented here differs from the approach of 1101 : while the latter 
proceeds by defining the composition of formal denotational semantics; we will adopt a somewhat more 
informal approach involving stable-index tensor notation. 

"The simplification implied in this title is twofold. Firstly, we use the Zz operation given in l |2.2b) to 
describe the construction of arbitrary operations, rather than the CNOT operation as in |2|, page 5. While 
this will lead in practise to a constant-factor reduction in the complexity of some IWQC procedures, this 
construction is clearly implicit in |2|, and so we may attiibute this construction to that work. Secondly, 
the "simplified" RBB construction omits the special-puipose constructions exhibited in Sections IV.A-F 
of |2|, such as compact patterns for reversing an block of consecutive qubits, non-nearest neighbor CNOT 
operations, the quantum Fourier transform over Zjfc for fc J5 0, and addition circuits for Zjfc . 
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{v} {w} 



l = = {v,w} 




Figure 3: Illustration of the geometries of the three elementary IWQC procedures 
5^^^, Evw, and 33u,u, for generic qubits v and w, given by ( |2.39| l, ( |2.40[ ), and ( |2.41| l 
respectively. Note that the input and output subsystems coincide for E^u, and 33t, 



for arbitrary operators ps acting on a system including v and w, and 



(2.43) 



for a system S including v but not including w, where by J{0)^/^ we denote the 
linear operator Tiy Hw which maps state-vectors \tp)v i — > [J{(^) IV') ■ Thus, 
the procedures Eyw, '53v,w^ and d^/^ given in ( |2.39| l, ( |2.40| l, and ( |2.41| l respectively 
may be used to represent each of the operations in Bdkp and Srbb- The geometries 



underlying these three IWQC procedures are illustrated in Figure 3 



2.3.2 Corresponding stable-tensor indices to qubits in the IWQC model 



By the analysis of |Section 2.3.1 we have a simple correspondence between unitary 
circuit operations described in terms of stable-index tensor expressions, and simple 
IWQC procedures for effecting them: we have 



u /v 



Zz 



(2.44a) 



We may extend this correspondence also to account for the preparation of fresh qubits: 
as a preparation map N„ simply prepares a qubit v in the |+) state, we may also write 



-r-] 



(2.44b) 



These parallels between operations on indices in stable-index expressions, and CPTP 
maps operations on qubits in the IWQC model, also extend to syntactical considera- 
tions: 

• If we order terms in a stable-index expression from right to left according to an 
order in which the operations may actually be performed, by construction there 
cannot be any terms acting on an index v to the right of an operation such as 
J(0)[^] which deprecates v; nor operations to the left of v before an operation 
such as in which v is advanced. Similarly, in a IWQC procedure, oper- 

ations acting on a qubit v are not permitted after it has been measured, nor before 
it has been prepared. 
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• In a stable-index expression for unitary circuits, stable indices (which are nei- 
ther advanced nor deprecated) are those indices representing the a qubit whose 
configuration over the standard basis are unchanged by the operator acting on it. 
Similarly, in a 1 WQC procedure, a qubit which is neither added to the system nor 
discarded is in both the input and output subsystems of the procedure, and so is 
only acted on by diagonal operations (setting aside possible X operations), which 
similarly do not affect the diagonal terms of their reduced density operators. 

Thus, the role of tensor indices (in a stable-index expression) in describing quantum 
states as distributions over the standard basis is analogous to the role of qubits in a 
IWQC procedure which actually support quantum states. We may then obtain IWQC 
procedures for unitary embeddings U by simply applying the correspondences in ( |2.44| i 
to a stable-index expression of a unitary circuit. While the the constructions of IdOl 
or ||2j| are not presented in precisely such terms, the DKP construction and simplified 
RBB constructions may be described as transliterating a sequence of gates in a unitary 
circuit into IWQC procedures of the form J^, E, and 33, yielding a IWQC procedure 



for the circuit as a whole. We present an explicit algorithm to do so in Section 2.3 

In the case of the DKP construction, in which the qubits of the resulting proce- 
dures are not subject to any topological constraints, application of the correspondence 
in ( |2.44| i is an essentially complete summary of the construction. The case of the sim- 
plified RBB construction is similar, but is subject to topological constraints which we 
describe below. 



2.3.3 Topological constraints in the simplified RBB construction 

The simplified RBB construction is subject to the constraint that the entangled states 
which are used must be obtained from a cluster state: that is, a state in which qubits 
involved are arranged in a two-dimensional grid, and two qubits ti, w are acted on with 
an entangling map E^^ if and only if v and w are nearest neighbors in the grid. We wiU 
describe circuit constructions in which the indices of the stable-index representation 
satisfies the same constraints as required for the qubits in the cluster-state-based vari- 
ant of the IWQC model, which will allow us to employ the correspondence of ( |2.44| l 
directly. 

In the cluster-state-based (and original) variant of the IWQC model, the only tool 
provided for obtaining graphs other than a rectangular grid is the removal of qubits 
from such a grid, via measurement in the |0) , |1) basis. Temporarily extending the set 



of operations defined in Section 2.2.2 we may define a operation which performs 
a destructive single-qubit measurement of a qubit v in the standard basis, and produces 
a bit s[v] storing the result in the obvious manner. We may then show 

M^, E„,, = Z=i"l M^, , (2.45) 

which follows from \x)s[v] {x\^ I^Z^^ — \x)s\v] for x G {0, 1}. We may use 
this to remove all of the entangling maps incident to a qubit v, propagating the 
measurement to the right until it is performed immediately after the preparation of v. 
If we remove the redundant preparation and measurement of v (and the corrections 
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Figure 4: Illustration of the effect of removing vertices from a two-dimensional grid. 
This corresponds to the effect of performing measurement of select qubits in a cluster 
state, in the standard basis. Vertices correspond to qubits, with vertices in grey corre- 
sponding to qubits which have been measured; edges (or absence of edges) correspond 
to the presence (resp. absence) of entangling maps on adjacent qubits. 



induced by that measurement), this procedure has the effect of removing the vertex v 
from the grid. This is illustrated in [Figure "4| 

We may instead suppose (equivalently) that we may prepare any subset of the qubits 
in a two-dimensional lattice (omitting those qubits which we would subsequently re- 
move as described above), and that two qubits are subject to an entangling map if and 
only if they are nearest neighbors in the grid. That is: we will merely require that 
the entanglement graph G is an induced subgraph of a rectangular lattice. Given this 
restriction on IWQC procedures, we impose the corresponding restriction that the cir- 
cuits which are accepted as input must have a linear nearest-neighbor topology. 

To represent products of single-qubit gates as a IWQC procedure, in the DKP 
construction as well as in the simplified RBB construction, we may compose patterns 
of the forni as in 

... ^ai ^ao ^ (2.46) 

We may refer to such a procedure as a chain pattern. Note that in such a procedure, 
we require that the qubits Vj for < j < to have degree 2 in the corresponding 
entanglement graph, having no neighbors aside from v^-x and Wj+i . In the simplified 
RBB construction, the degrees of these qubits must be "enforced" by a representation 
in the grid, where every neighbor of these qubits Vj in the grid (for < j < N) are 
removed except for the neighbors and t'j+i. An illustration of such an embedding 



of a chain in the grid is illustrated in Figure 5 (Typically, such chain patterns would 



be embedded as a horizontal path through the grid, but as Figure 5 also shows, we may 
also employ more general paths in the grid.) 

Two logical qubits in a circuit which do not interact must be represented by (chains 
of) qubits in the grid which are separated by at least two edges; otherwise there will 
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Figure 5: Illustration of an embedding of a chain in the grid. Grey dots and broken 
lines represent vertices and edges which are removed from the grid by vertex deletions, 
to produce a graph containing the illustrated sequence of vertices with degree 2. 



be entangling relations between them. Thus, we will represent products of single-qubit 
operations on independent qubits as non-adjacent horizontal rows of the grid. We may 
then use a vertical column of three qubits (as illustrated in Figure 3 1 to implement a 
Zz operation between these non-adjacent rows, using a single qubit between the two 
rows to mediate the interaction. However, there are two technicalities which must be 
addressed: 



Adjusting for uneven chain lengths. Consider a decomposition of single-qubit op- 
erators U and V on distinct qubits, as products of J(0) operators. If the number of 
terms in the two products differ, the ends of the corresponding chain patterns may lie 
in different columns. To operate on two chains with a Zz operation as described above, 
we require that the ends of the chains representing each qubit lie in the same column. 
To achieve this, we may extend each chain by two or more vertices, as necessary, rep- 
resenting a decomposition of 1 as a product of J(ff) gates. We may define IWQC 
procedures 

n3 



at) 



'"V3/V2 '^V2/vi '"vi/vo 



(2.47) 



the geometries underlying these procedures are illustrated in Figure 6 In terms of the 



correspondence between stable-index expressions and IWQC procedures described in 
[Section 2.3.2| we have 



(2.48a) 



ao 



J{-k) ^ J{-h) ^ J{-h) 



(2.48b) 



One can verify the stable-index products given on the right evaluate to and 
respectively; then both IWQC procedures above simply transmit the state presented 
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^^l/vo ® • ® 

/ = {vo} O = {V2} 



-V2 "V2 "V2 

/ - {vo} O = {V3} 



Figure 6: Illustration of the geometries underlying the procedures 3d^ and 3d'^ as de- 
fined in ( |2.47| i. The default measurement angles for qubits in O'^ are indicated above 
each qubit. 



as input to its respective output system, using a chain of length either 2 or 3. By 
composing these procedures repeatedly, we can then implement a path of any length 
£ ^ 2 in the grid, representing a sequence of operations which performs the identity 
I2 on a given qubit. Thus, whenever required, we may suppose that the paths in the 
grid corresponding to any two qubits are of the same length, regardless of the number 
of (non-trivial) unitary transformations performed on them. 



"Forbidden" qubits adjacent to the mediating qubit in Zz. Just as we require a 
distance of 2 between qubits in a IWQC procedure representing non-interacting qubits 
in a unitary circuit, we require that the qubit a which mediates a Zzy^^ operation is 
adjacent only to v and w; that is, we require that the qubits immediately to the left 
and right in the grid are not in the entanglement graph. Then we require a distance of 
two edges between any two instances of 33 performed on the same pair of rows. For 
instance, to represent a circuit ^^^[f/,!"/] Zz[vi,wi'^ on a pair of 

logical qubits v and w, in which two Zz operations are separated only by a layer of 
single-qubit operations of depth 1, we may perform the substitution 



Zz[v,.^Of]J{P)[^^ J{a) 
= Zz[vj,Wf]J[Q) 

J(0) 



Zz[vi,w 

J(0) 
J(0) 



m 



Zz\vi 



(2.49) 



Mapping each term to a procedure Z^a/d ^^'^ &^ch Zz[s,t'^ term to a procedure 

33s, t , the distance between the qubits Vi and Vf resulting in the right-hand side is then 
at least 2 (and similarly for Wi and Wf), as required. 



These observations describe ways in which a IWQC procedure for performing 
transformations must be "padded" by operations which perform the identity, in or- 
der to simulate a unitary circuit in the IWQC model. We illustrate these observations 
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in Figure 7 as sample constructions of IWQC procedures. However, in order to ob- 
tain a uniform description of both the DKP and simpHfied RBB constructions, we may 
instead pad the stable-index representations for circuits with the corresponding prod- 
ucts of single-qubit unitaries as a pre-processing stage. By construction, the indices of 
the resulting stable-index expressions will then satisfy the same constraints which are 
required for the qubits in a cluster-based IWQC computation. We may perform such 
padding as part of an alternative procedure for constructing stable-index expressions 
for circuits: such a procedure is described in |Section A.2.2| 

2.3.4 Unified description of tlie DKP and simplified RBB constructions 

Using the above techniques for the simplified RBB construction, we may obtain a uni- 
fied procedure describing both the DKP and simplified RBB constructions, as follows. 

Procedure for botli constructions. Suppose we are given a unitary circuit C over 
the gate-set {iJ, T, ^Z} : as described in the previous section, we further require that C 
have a linear nearest-neighbor topology in the case of the simplified RBB construction. 

1. Obtain a circuit in normal form. 

Given a unitary circuit C generated over the gate-set {H, T, ^Z}, we obtain an 
equivalent circuit as follows. We cancel pairs of gates which act on common 
pairs of wires, and which are separated only by gates with which they commute 
{i.e. T operations and operations acting on other wires). In the case of the sim- 
plified RBB construction, we convert the circuit to the (essentially equivalent) 
gate-set {H, T, Zz} by applying the substitution 

= Zz(T^®T^)2 (2.50) 

to all AZ gates; no such translation is required for the DKP construction. We 
then commute all T and operations to the left (simplifying them using the 
identity = 1), collecting them in each case either at the end of the circuit, or 
immediately preceding a Hadamard operation the same qubit. 

We call the resulting circuit C the normal form of C; this circuit consists of a 
product of operations and HT^ (for various r E Z) for the DKP construction, 
and Zz and HT^ for the simplified RBB construction. 

2. Convert the circuit to a stable-index expression in the gate-set Sdrp or /Brbb- 

Having obtained C", we apply the substitutions 

iJT™ = J(^^) , T" = J(0)J(^) (2.51) 

for each consecutive block of T operations, applying the right-hand equality for 
blocks of T operations at the end of the circuit. These substitutions yield an 
equivalent representation of the same operation using the gate-set ;Bdkp or Srbb, 
depending on whether the circuit makes use of ^Z or Zz gates. As we do so, we 
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Figure 7: Illustration of the constraints on composition of the pattern 33 in the setting 
of an induced subgraph of the grid. Grey dots represent vertices which have been 
removed from the grid. Composition is illustrated with overlapping tiles, whose shapes 
represent spatial constraints imposed by the operations. 

(a) The geometry itself, together with the "forbidden" neighbors of the mediating qubit a. 
(Broken lines between the dots represent the edges removed as a result of removing these 
vertices.) 

(b) A composition of 33 procedures with procedures, corresponding to a circuit in which 
two Zz gates are separated by two J (9) operations. The bottom row of the corresponding 
IWQC procedure is padded with a procedure to synchronize the lengths of the rows. 
In this case, a single vertex is a common forbidden neighbor of each mediating qubit. 

(c) A composition of 33 procedures with 3* procedures, corresponding to a circuit in which 
two Zz gates are separated by one J{6) operation on either wire. The forbidden neighbors 
of the mediating qubits force both rows to be padded with identity operations. 
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convert the resulting circuit into stable-index notation: we use an algorithm such 
as that of Section A. 2.1 to do so for the DKP construction; for the simplified 
RBB construction, we may use instead the procedure of |Section A.2.2| to impose 
additional topological constraints on the resulting circuit. We call the resulting 
stable-index expression for the circuit C" . 



3. Translate individual operations into the IWQC model. 

In C" , note that each single-qubit J{9) operation deprecates one index and ad- 
vances another, and that all of the two-qubit /\Z or Zz operations act "stably" on 
their indices. We may then define a mapping <1> from stable-index expressions of 
operations in the gate-set Sdkp U Srbb (and introduction of fresh qubits in the 
1+) state) to IWQC procedures involving J^, E, and 33 operations as follows: 





) = ; 


(2.52a) 




\ - ■ 


(2.52b) 






(2.52c) 






(2.52d) 



where we identify the labels of indices of the stable-index tensor expression with 
labels for qubits in the IWQC procedure. We then extend this map $ homomor- 
phically to products of such terms. 



4. Obtain a normal form for the resulting IWQC procedure. 



We standardize (Definition 2.5 1 the IWQC procedure <i>(C"), and then apply 
Pauli simplifications and signal shifting to obtain a procedure in normal form 
( [Definition 2.6\ . 

By construction of C" in the above procedure, each index is advanced at most once 
and deprecated at most once; thus, each qubit in <i>(C") is produced in the output of 
an operation (or "allocated") at most once, and removed from the input of an operation 
(or "discarded") at most once. Furthermore, by order of the terms in C" , the operation 
which allocates a qubit v is the first which acts on v, and the operation which discards 
V is the last acting on v. The resulting IWQC procedure is therefore well-formed. 



By the analysis of Section 2.3.1 we may easily show in either case that <i>(C") 
is a CFTP map which performs the same unitary transformation as described by C" , 
modulo an identification of each input qubit with a corresponding output qubit. Then 
the final procedure produced performs the same transformation as the circuit C. We 
refer to this procedure as the simplified RBB construction for IWQC procedures when 
we require that C be linear nearest-neighbor and when we perform the substitution of 
( |2.50| l, and the DKP construction for IWQC procedures when we do not impose the 
added constraint or perform the substitution. (For the sake of brevity, we will often 
refer to this construction simply as "the RBB construction", as we do not consider any 
constructions using the special-purpose IWQC procedures in Section IV of ||2|.) 
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Run-time complexity. Both the DKP construction and the (simpHfied) RBB con- 
struction can be performed efficiently in the size of the input circuit. In both cases, let 
k be the number of qubits which C acts upon, N be the number of one-qubit gates, 
and AI be the number of two-qubit gates: if we consider only those circuits which act 
non-trivially on each qubit, we may assume that k ^ N + M for simplicity. We may 
then show that the run-time complexity of the DKP construction is 0{{N + k)M), and 
of the RBB construction is 0{N AI {N + M)) , in each case bounded by the complexity 
of obtaining a IWQC procedure in normal form. (The difference in complexity be- 
tween the two constructions may be traced to the difference in the procedures to obtain 
the stable-index expression C" for each, and hence to the single-qubit gates which are 
added in the RBB construction in order to impose the desired topological constraints 
on the indices of the stable-index expression.) An explicit analysis of the complexity 



of these constructions is presented in Section B.2 



2.3.5 Partial constructions 

As a final remark on constructions of IWQC procedures, it will be useful to describe 
partial versions of the DKP and (simplified) RBB constructions. We define the DKP 
construction without normalizatior^ and the RBB construction without normalization 
to be, simply, the IWQC procedure constructions which result from omitting the Pauli 
simplifications and signal shifting stages from the phase of obtaining a IWQC proce- 
dure in normal form. A procedure produced by such a partial procedure may then indi- 
cate sign dependencies for a measurements whose default angle is an integral multiple 
of "72, and they may have signal-shift operators immediately following measurements. 

In a practical setting, provided sufficiently fast control of the classical memory stor- 
ing the measurement results s[d] for all qubits w in a IWQC procedure, it is likely that a 
more useful procedure would be an intermediate construction to the constructions with 
or without normalization, in which one performs Pauli simplifications but leaves the 
shift operations in place. The purpose of identifying the partial constructions without 
normalization above is to obtain a simplified analysis of measurement dependencies, 
which may be used as a starting point for the analysis of the "complete" constructions. 



3 Semantic maps 

We may now consider the roles of the DKP and (simplified) RBB constructions as 
defining representations of unitary circuits in the IWQC model. Both constructions 
make use of a simple process of translating unitary circuits from a gate set including 
J{9) gates and a single two-qubit diagonal operation into corresponding IWQC proce- 
dures: this will induce a graph structure which is dominated by vertex-disjoint paths, 
with supplemental vertices and edges linking pairs of paths. After standardizing the 
IWQC procedures in the course of these constructions, the additional structure which 
might arise from measurement dependencies or the grouping of operations is lost or 

PThis construction is refen'ed to in (9) as the "simplified" DKP constiTJCtion; we use different terminology 
here in order to be more uniform with our description of the RBB construction. 
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obfuscated; preparation maps and entangling operations become dissociated from in- 
dividual measurement operations, and corrections will in some cases be absent entirely 
after absorbing them into measurements and performing Pauli simplifications. 

Given a IWQC procedure, in which no additional information is encoded e.g. in 
the labels for the qubits, we may ask whether we can efficiently identify whether a 
IWQC procedure is one that may be produced by the DKP or the RBB constructions. 
In the case where the input and output subsystems of the procedure are of the same size, 
we answer this question in the affirmative. We do so by considering the combinatorial 



structures which are present in the geometries ( Definition 2.7 1 of the IWQC procedures 
which result from these constructions. We then define a map S from IWQC procedures 
to circuits in the {H, T, AZ} model which serves as a semantic map (in the sense of 



Definition 1 . 1 1, for the representation maps 7?.dkp and 7?^rbb corresponding to both the 
DKP and simplified RBB constructions restricted to input circuits which are unitary 
bijections. 

3.1 Measurement dependencies arising in the constructions 

The DKP and RBB constructions produce IWQC procedures from circuits by trans- 
forming gates such as J{0), ^Z, and Zz into simple IWQC subroutines, via the map- 
ping <i> defined in ( |2.52| i. It is instructive to consider the structures described by the 
sign- and bit-dependencies which arise from composing such IWQC subroutines, us- 



ing the partial constructions without normalization described in Section 2.3.5 



3.1.1 Measurement dependencies and index successor functions 

For a unitary circuit C consisting of gates from the set Sdkp = {•^(^)i ^-^lee-Z' 
consider the last gate J{0) performed in the circuit. (If there is more than one such 
gate which may be performed in parallel, we may choose an arbitrary one.) We may 
decompose C into circuits C^C, where C^, consists of this final J{d) gate acting on 
some qubit u, and any controlled- operations acting on u which follow it; C consists 
of the remaining operations of C. We may refer to C\ as a star c/rcM/ij^] Consider the 
translation of a star circuit into the IWQC model via the map (f> defined in p.52| i. 



and standardizing the resulting procedure (Definition 2.5 1. Writing C^, in stable-index 



tensor notation (Section 2.1.3 1, and letting W be the set of qubits w ^ uon which C^, 
acts, we have 



\wew 



These star circuits are similar to, but have a different orientation to, the circuits described in Section III B 
of ll2l . This different choice of construction will facilitate the analysis of combinatorial structures later in 
the article: these are used in much the same way as the circuits arising from "star patterns" do in lT2l . 



31 




(3.1) 



where on the last line we commute the entangling operations to the right. In the final 
IWQC procedure, the qubit u is distinguished as the only qubit measured, and v is 
the only qubit prepared; v is also the only qubit subject to an X^'"' operation, with 
every other qubit (aside from u) being subject to a Z^™' operation. The fact that 
u and V are measured and prepared qubits (respectively) in <i>(C^) above is a direct 
result of the fact that they are deprecated and advanced indices (respectively) in the 
stable-index expression given for C*. In particular we have v = f{u), where / is the 



index successor function (Definition 2.3 1 of C*; thus, the corrections may be described 
in terms of the index successor function of the original circuit, and the entangling 
operations. 

We may make similar observations for the (simplified) RBB construction as we 
have for the DKP construction above. Motivated by the realization of the operator 
Zzvw in the IWQC procedure 33„,«„ we may write 

Zz[v,w] OC P[-] AZ[v,a] AZ[w,a] , (3.2) 

where we define the operator P = {+tt/2 I . and for a constant proportionality factor (as 



illustrated in ( |B.4| i in Section B.1.1 1. In the stable-index expression on the right, the 
index successor function / remains the same (as the indices a are advanced without 
corresponding to any deprecated index, and vice versa), but the interaction graph G 
is transformed essentially by subdividing every edge arising from a Zz gate by intro- 
ducing a new vertex in the middle (representing the mediating qubit with the index a). 
Then, consider the translation of ^Zu.^ to a IWQC procedure via <i>, 

<^(^Zz{e)[u,v] ) = Zf 1 Zfl : (3.3) 

if we extend the successor function / to a function / for which we define /(a) = a, 
while we have noX^"' operations, we have aZ^"' for every qubit w adjacent to a = / (a) 
in the entanglement graph. 

In each case, either the index successor function (or a modest extension of it) for 
a stable-index description of a simple unitary circuit can be used to characterize the 
corrections performed in a corresponding IWQC procedure. We may then consider the 
result of translating a composition of such simple circuits, 

$(0 = <i>(cP...cW) = $(cP)o...oci>(c«) , (3.4) 

where each cl'''' is a star circuit as above, extending the definition to also include 
individual Zz operations as described in (|3.2[). Commuting all of the preparation and 
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entangling operations to the right, we obtain the same dependency of the corrections on 
the successor function (extended for each Zz gate as described above): we obtain Z^"' 
corrections for each qubit in the interaction graph which is adjacent to /(u), either 
due to the commutation of entangUng maps past X^'"^j operations or from the Z^"' 
operations arising from 33 procedures. 



3.1.2 Efficiently verifying consistency of measurement dependencies 

Based on the descriptions of the dependencies arising in the DKP and RBB construc- 
tions without normalization, we may characterize the dependencies of operations which 
arise in the two constructions with normalization, and show that verifying that such de- 
pendencies hold in a IWQC procedure *p can be used as a subroutine to certify that ^ 
performs a unitary transformation. 

Following the remarks of the preceding sections about the constructions without 
normalization, we characterize the dependencies arising from the DKP and RBB con- 
structions as follows. 

Definition 3.1. For a stable-index representation C for a circuit the gate set Sdkp U 
^RBB, let C be the stable-index expression obtained by performing the substitution of 
p.2[ ) for each Zz operation in C. Let / be the index successor function of C, as defined 
in 



Definition 2.3| The extended successor function of C is then the function 



f{v) = / ^^'"^ ' " ™ ^"'^^^ °^ ~ 1 (3 5) 

[ w , if u is an index deprecated in C* by a P operator J ' 

whose domain is the set of indices of C which are deprecated. 

Theorem 3.1. Let *Po = ^(C*), for a unitary circuit C over either the gate set 
Sdkp U Srbb, and let ^ be the result of normalizing *Po. J^t {G,I,0) be the ge- 
ometry underlying *Po. <^nd let f be the extended index successor function of C. For 
the sake of brevity, define Ax Q as the set of qubits to be measured with default 
angles Oori: in *Po. <^nd Ay C O'^ the set of qubits to be measured with default angles 
^'"72 in *Po. define square boolean matrices F, T : V{G) x V{G) — > Z2 in terms 
of their coefficients, as follows^ 



ifv e and w = f{v) Ax U A^ 
otherwise 



(3.6a) 



ifw ^ V E O'^ and w ^ f{v) 
Ty^ = 1 , ifv, w eO'^ and w = f{v) £ Ay } ■ (3.6b) 
otherwise 

Then (1 — T) is invertible modulo 2, and the dependencies in *P between qubits v,w € 
V{G) may be characterized in tenuis of F and T as follows: 



' These operators are the transposes of the operators described in [9J, corresponding to changes in the 
representation in order to simplify the description of the effect of signal shifting; the statement of the result 
is also changed accordingly. 
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' has a sign-dependency on v if and only if [(1 — T) ^F] yyj = 1, for w G O'^; 



• there is an operation conditioned on s[v] if and only if[{l — T) ^-Fjum = 1, 
and a operation conditioned on s[u] if and only if[{1 — T)^^T]yiu = I, for 
w £0. 

Proof. By construction, for qubits u,w e V{G), ^Po contains a Zt,'"' operation if and 
only if t; ^ f{u) — u; as Z^'"' operations commute with entangling maps, we obtain 
no further corrections depending on u if we commute the entangling maps of *Po to the 
right. Similarly, *Po contains an X^'"' operation if and only if u = f{u) ^ u, which 
must arise from a 3^/„ procedure. Any operation on v in *Po which is not a part of the 



/u 



procedure must come after it, every entangling map acting on v must occur after 
the X^'"' operation except for the map occurring as a part of 5^/^. 

Commuting all of the entangling maps in *Po to the right, we then induce Z^,"' 
operations for all w ^ v such that w 7^ it, as described in ( |2.23c[ ). Thus, commuting all 
preparation and entangling maps in *Po to the right yields an equivalent procedure 
in which for there is an X^"' operation for every w = f{v) ^ v, and a Z^"' operation 
for every w f{v) 7^ v, for every qubit v G C^. That is, we have 



n 
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( \ 
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-7S[u] 
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v=f{u) 





(3.7) 



for some linear order ^ satisfying the conditions of ( |2.14| i for the extended successor 
function /. 

Commuting correction operations in *Pi to the left, we obtain another procedure 
*P2i in which the same corrections hold as above for w G O. For w G 0°, there is 
instead a sign-dependency on v if and only if it; = f{v) 7^ v, and a bit-dependency on 
V if and only if w v and w ^ /(«)• Note that the procedure of producing ^2 from 
*Po is precisely that of standardizing ^q: then *p may be obtained by normalizing 

Consider how the process of normalization effects dependencies of qubits w on a 
particular qubit v. 

(i) The effect of Pauli simplifications are to remove sign-dependencies from qubits 
in Ax according to ( 2.34f| i, and to change them to bit-dependencies for qubits in 
Ay according to ( 2.34g| . As signal shifting does not introduce sign-dependencies 
for measurements where none previously exist, the measurement of a qubit in 
Ax U Ay will not have any dependencies on previous measurements in a pro- 
cedure in normal form. 

Applying Pauli simplifications to '^2 will yield a IWQC procedure then 
*P3 will have no sign-dependencies for qubits v G Ax U Ay, and will have 
signal shift operators S^j"] for any qubits v = f{u) G Ay, in addition to the 
shift operators in ^2- We may describe the sign dependencies which remain 
in by a square matrix over V{G), with a 1 in the row v and column w for 
VjW G V{G) when w has a sign-dependency on v. Accounting for the removals 
of sign dependencies of ^2 for qubits w G Ax U Ay, the procedure then 
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has a sign-dependency of w on w when w = f{v) ^ Ax U Ay- The matrix 
described in this manner is then the matrix F in ( |3.6a| i. 

For the sake of simpUcity of discussion, we will assume that the shift operators 
are not accumulated with any pre-existing shift operators using the relation of 
( |2.35| l. However, it is important to note that the shift operators arising from 
Pauli simplifications will not cancel any operators existing originally in *P2: if 
w ~ f{v), it follows that w is not adjacent to f{v), as G contains no loops by 
construction. 

(if) The effect of signal shifting is to propagate a bit-dependency of each qubit w 
on V (represented by a S^j^'j operation just after the measurement of w) to the 
operations which in turn depend on the value of s[w], according to p.36| l. 

We may describe the effect of signal shifting in terms of walks in a directed graph 
D of dependencies, where we have arcs w ^ v when there is a shift operator 
S^l^'j . An operation depending on w may be transformed into one depending on 

w and all qubits w for which D has an arc w ^ u by shifting of operators S^j^j ; 
in turn, the resulting operation may be transformed into one depending also on 
qubits u for which D has arcs w ^ u, ranging over the qubits v of the previous 
step; and so on. The effect is then of transforming bit-dependencies correspond- 
ing to arcs in this digraph D (arising from signal shifts) into measurement and 
correction dependencies corresponding to directed walks in D. 

We may represent this digraph as an adjacency matrix over V{G), with a 1 in the 
row V and column w for w, w e V{G) when contains a shift operator S^j^'j , 
and if not. Considering the shift operations in which are also present in 
*P2, there will be such a shift operator whenever w ^ v and w ^ f{v); and as a 
result of Pauli simplifications, there will be such a shift operator when w — f{v) 
and w S Ay- As there are no other sources of shift operators, these are the only 
positions in which the adjacency matrix will be non-zero. The adjacency matrix 
described is then the matrix T in p.6b| l. 

As the ordering ^ (used in p.7| l to fix the order of operations) is a linear order, 
the digraph D described by the shift operations is acyclic: for any standard basis 
vector e„, e the column- vector Te^, is supported only on indices v ^ w. 

Then we have — for some r ^ |1/(G)|, so that T is nilpotent (and in par- 
ticular has no +1 eigenvectors). Consequently, the operator 1 — T is invertible. 
We may then express the cumulation of dependencies for each operation in 
due signal shifting as follows. For two qubits v and w, the number of walks in 
D from ui to t; of a fixed length ^ ^ is given by e^T^e^ : the total number of 
walks of any length from w to u in is then 



If we represent the dependencies of operations in terms of boolean column vec- 
tors, we may then represent the effect of signal shifting by multiplying a vector 
of dependencies by (1 — T)^^ computed modulo 2. 




(3.8) 
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To determine the classical dependencies of each operation in *p, it suffices then to 
identify the dependencies of the same operations in ^^3, represent these as a vector 
d e 'L2^^\ and compute {l-T)-^d. We may characterize the classically controllable 
operations and their dependencies in as follows: 

• As we have noted above, the sign-dependencies for qubits w E are repre- 
sented by the coefficients of the matrix F, where = 1 if and only if w has 
a sign dependency on v in *P3 . Then w has a sign dependency on w in *p if and 
only if = 

• For w E O, there is no measurement into which corrections on w may be ab- 
sorbed, and so the corrections in *P3 are the same as in these are X^"' 
operations when w = f{v), and and Z?|,"' operations when w ^ f{v), where in 
either case we require w v as well. As O is disjoint from both Ax and Ay 
by definition, we may then note that for w E O there is an X^^' operation in ^^3 
if and only if F-u^, = 1, and a Z^"' operation if and only if T„tu — 1. It follows 
that after signal shifting, the resulting X(^ operation depends on v if and only if 
[(1 — T)~^F]„uj = 1, and the resulting ZJ, operation depends on v if and only 
if[(l-T)-iT]„^ = l. ■ 



Theorem 3.1 [thus characterizes the dependencies which arise in the DKP and RBB 



constructions, in terms of an (extended) index successor function for a circuit C over 
the gate-set Sdkp U Srbb- 



[Algorithm 1 describes an algorithm which, provided a IWQC measurement pat- 
tern *p in normal form and a candidate for the extended successor function / for an 
originating circuit, tests if the dependencies of *p are consistent with the dependency 
conditions described in [Theorem 3.11 

If we are presented with a IWQC procedure *p in normal form which is not known 
to be a result of the DKP or RBB constructions, we may only verify whether these 
conditions hold relative to some "candidate" successor function /. How such a function 



/ may be constructed for a IWQC procedure *p is the subject of Section 3.2 



Run-time analysis of Algorithm 1. We bound the run-time of [Algorithm 1 as fol- 



lows. In the following, we let N be the number of operations in *p, and for the geometry 
(G, /, O) underlying % we let n = \V{G)\, k = \0\, and m = \E{G)\. We assume 
that / can be evaluated in constant time, using an array structure. We assume that the 
dependencies of operations Op^ (for sign-dependencies or other classical control ex- 
pressions) are represented as indicator vectors d E Z^'*^-* representing the presence or 
absence of the term slv] in /3 hy dy — 1 (with d^, = otherwise). 

By definition, *P will contain m entangling operations, n — k measurement opera- 
tions, and at most 2k correction operations; then, the length of *p is e 0{m + n) C 
0{n^). Rather than compute matrix products involving (1 — T)^^, we may verify 
whether (1 — T)d = v for v = Fe^ or v = Te^, as appropriate. The coefficients 
of (1 — T)d may be computed by iterating over the neighbors of f{w) for w S*-^-* d): 
as as Tyjju — 1 only for w adjacent to f{v) (if v ^ w) or for w = f{v), there are at 
most two non-zero coefficients in T for each edge in G, so that this iteration requires 
time 0{m) for each vector d. Similarly, the time required to compute column-vectors 
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Input : *P, a IWQC procedure with geometry (G, /, O); 

/ : O*^ — > P, a candidate successor function for *p. 
Output: A boolean result indicating if the dependency relations in *p are 

consistent with [Theorem 3.1| relative to the function /. 

; procedure TestDependencies(*P, /) : 

Let r be a matrix depending on / as in ( |3.6a[ ). 
Let F be a matrix depending on / as in p.6b| l. 
Letf ^ (l„-r). 

For each measurement or correction operation Op /« *p, do: 
If Op = M^'' for some w, 0, (3, then 

7 Let d e Z^^*^^ be the indicator vector for (3. 

8 If Td ^ Feu], then return false. 

9 else, if Op = Xf^ or Op = for some w, /?, then 

10 Let d e Z^'*^' be the indicator vector for /3. 

11 If (Op = and fd ^ Fe,,) or (Op = Zg and fd ^ Te^) 

12 then return false. 

13 endif 

14 endfor 

15 return true. 



Algorithm 1. An algorithm to test whether the dependencies of a IWQC procedure 
is consistent with those resulting from the DKP and simplified RBB constructions, rel- 
ative to a candidate for the (extended) index successor function for the original circuit. 



Fe^ and Te^ for each w is 0{m). As there are 0{n) measurement and correction 
operations, the cumulative time to perform these matrix computations is 0{nm); as 
the remaining operations in *p require no verification, this is the total time complexity 



of the for loop. The execution time of Algorithm 1 is then 0{nm) 



3.2 Conditions on candidate successor functions 

In the preceding section, we characterized the dependencies in the DKP and RBB con- 
structions in terms of an extended index successor function, from a stable-index repre- 
sentation of the original circuits. This leaves open the question of how such a successor 
function may be obtained. We examine this problem using the observations made above 
about the dependencies arising from circuit decompositions without normalization. We 
show that these dependencies can be captured in each case by a modification of the 
flow conditions formulated by Danos and Kashefi lfT2l . This modification will provide 
the "candidate successor functions" which we require for the verification procedure 



of [Algorithm 1 and which lead to important structural constraints on the geometries 



(G, /, O) arising from the DKP and RBB constructions. 
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3.2.1 Measurement ordering from index successor functions 



As we described in Section 2.3.2 there is a correspondence between the "deprecation 
order" of indices as described in p.l4| l, for a stable-index expression for a circuit C 
over /Sdkp U Srbb- ™d the order of operations in a well-formed IWQC procedure 
<I>(C). This correspondence may be attributed to a sense in which measuring qubits in 
a IWQC procedure simulates the summation over indices in a tensor expression, and 
thus to evaluating a unitary circuit applied to a state as a sum over computational paths. 

When applied to a completely specified pure state, every deprecated index in a 
stable-index expression becomes bound, in which one evaluates the result by summing 
over all deprecated indices. When these indices are represented by qubits in a IWQC 
procedure, one may simulate this summation by measurement of qubits v with respect 
to the l+e) , |— e) basis. 

• Consider a measurement arising from a procedure, corresponding to summing 
over the deprecated index of a J {6) operation. If the \+g) state is measured (i.e. if 
s[u] = 0), the resulting transformation corresponds to summing over the |0) and 
|1) components, weighted with the with the relative amplitudes phases as those 
associated with the J {6) operation which deprecates the corresponding index. 

• For measurements arising from a 33t, ^ procedure, we measure a mediating qubit 
a in a fixed basis |±^^). Again, a measurement result of I+tt/s) corresponds to 
summing over an auxiliary index a in a stable-index expression: the relative 
phases of the |0) and |1) components correspond to those of the projection oper- 
ator P. 

In either case, if the measurement yields the |— e) result instead, the corrections and 
measurement adaptations are chosen precisely to effect the same transformation which 
would have occurred for the result, by influences on the measurement angles due 
to sign-dependencies or the measurement results due to bit-dependencies on 

Thus, we may interpret the sign- and bit-dependencies in such IWQC computations 
as potential "influences" in evaluating the summation of deprecated indices. By con- 
struction, the measurements are chosen in such a way that the tensor which is produced 
at the output is related to the tensor at the input by a unitary circuit: the amplitudes 
arising from each measurement are a function of at most two boolean indices, where 
the two-index dependencies are signs encoded by the entanglement graph, and whose 
single-index dependencies are given by the choice of measurement basis. This may be 
regarded as indicating that a IWQC procedure where such an ordering is possible has 
a "circuit-like structure". 

A measurement order which corresponds to a "deprecation order" for indices in a 
circuit C as in ( |2.14| l is governed by the index successor function of C. An extended 



index successor function as defined in Definition 3.1 has the effect of imposing con- 



straints on the deprecation of the mediating indices a arising from the substitution of 



"Bit-dependencies, as represented by sliift operators S*|^j immediately after measurements, may be in- 
formally interpreted as flipping the bit s[v\ to produce "the result which would have been measured if an 
appropriate Z correction had been performed prior to measurement". As the measurement results are of- 
ten maximally random whether or not such corrections are performed, this is a somewhat counterfactual 
interpretation, but one which could perhaps be given an ontological foundation in terms of hidden variables. 
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Figure 8: Examples of geometries with flows. Arrows indicate the action of a function 
/ : 0° — > P along otherwise undirected edges. Corresponding partial orders =^ for 
each example are given by Hasse diagrams below the graphs (with minimal elements 
on the left and maximal elements on the right). In the right-most example, the two 
vertices a and b are incomparable, i.e. there is no order relation between them. 



p.2[ ) for Zz gates, which may be interpreted as creating a well-defined notion of when 
the Zz operation is performed. By construction, IWQC procedures arising from the 
DKP and RBB constructions do admit extended successor functions: this can be used 
to define the domain of a semantic map for these constructions. 

3.2.2 Flow conditions 

The ordering in ( |2.14| i was first formulated by Danos and Kashefi fT2l, for qubits in the 
IWQC model rather than tensor indices, to describe dependencies arising in the DKP 
construction without normalization. We may formulate these properties in terms of the 
geometry of a IWQC procedure without reference to the DKP construction as follows: 

Definition 3.2. For a geometry (G, /, O), let u ^ v denote the adjacency relation in 
G. A flow is a pair (/, =4) consisting of an injective function / : 0° — > P, and a 
partial order ^ on V{G), such that the conditions 

V ^ f{v) , (3.9a) 

V 4 f{v) , and (3.9b) 
w ~ f{v) =4« V ^ w (3.9c) 

hold for all v eO^mdw e V{G). 



Examples of geometries with and without flows are illustrated in Figure 8 and Fig 



ure 9 respectively. The function / and the partial order =^ capture the essential structure 
of the dependencies in the DKP construction without normalization: / represents the 
mapping of qubits to their successors (corresponding to an index successor function for 
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Figure 9: Example of a geometry with no flow {c.f. the second geometry from the left 
Figure 8 i. Arrows indicate the action of an injective function / : 0° — * P along 



m 



otherwise undirected edges. Also given is a reflexive and transitive binary relation 
=<; which satisfies conditions p.9b| i and p.9c| l for this function /, but which is not 
antisymmetric. 



a circuit), and the partial order ^ represents a suggested order in which the qubits may 
be measured (or the tensor indices deprecated). 

Danos and Kashefi show that any geometry which has a flow underlies some 1 WQC 
measurement pattern which can be obtained by the DKP constructionj^ and which 
therefore performs a unitary transformation: 

Lemma 3.2 (fT2l, Theorem 1). Suppose (/, ^) is a flow for {G,I,0). Let v ^ w 
denote the adjacency relation of G: then for any linear order ^ extending =4> the 
measurement procedure 



( 



n 



n 



EgN/c 



(3.10) 



performs a unitary transformation. 

The proof presented in [12 | is essentially by showing that applying the map $ 
defined in p.52| i to a suitably constructed circuit yields the IWQC procedure above. 
If we can construct a flow for a geometry (G, /, O), we may then use the same circuit 
construction as in 1 12 1 to construct a candidate circuit whose representation as a IWQC 
procedure (via the DKP construction) has the same geometry. We may then verify 
whether the dependencies of a IWQC procedure are consistent with one arising from 



the DKP construction, using Algorithm 1 



'The converse of this statement does not necessarily hold: by omitting or manipulating the dependencies 
in a IWQC procedure, a variant procedure which does not perform a unitary may be obtained without 
changing the underlying geometry. 
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3.2.3 Extending to modified flows 



By considering the dependencies described in the proof of Theorem 3. 1 following ( |3.7| l, 
we may formulate a similar combinatorial condition which generalizes flows, and also 
captures the dependencies arising in the (simplified) RBB construction without nor- 
maUzationj3 

Definition 3.3. Let (G, /, O, M) consist of a geometry (G, /, O) and a subset M C 
C^, and let u ~ u denote the adjacency relation in G. A modified flow for (G, /, O, M) 
is an ordered pair (f,^) consisting of an injective function / : O*^ — > and a partial 
order =<; on V{G), such that the conditions 



f{v) V , for i; e M° ; 
either f{v) — v or f{v) ^ v , for w e Af ; 
V 4 f{v) ; and 
w ^ f{v) =4> V 4 w 



(3.11a) 
(3.11b) 
(3.11c) 
(3.1 Id) 



hold for all v eO^andw e V{G). 



The function / again corresponds essentially to an (extended) successor function 
for stable-index expression of a unitary circuit, and =<; again represents an order of 
measurement corresponding to an order in which the indices of stable-index expression 
are deprecated. The set M represents qubits whose measurement angles are "72 : these 
are qubits which might be mediating qubits a in the 33i, w procedures. We may recover 
the definition for flows in Definition 3.2 if we either set M = 0, or simply require that 
V 7^ f{v) for all v. 

The following generalization of Theorem 1 of fT2l to modified flows provides the 
basis for the approach we will take to define semantic maps for the DKP and RBB 
constructions: 

Lemma 3.3. Suppose (/, =<;) is a modified flow for (G, /, O, M). Let v ^ w denote 
the adjacency relation of G: then for any linear order ^ extending ^, the measurement 
procedure 





( ] 
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n 




n ' 






iu~/(tt) 
\ w^u / 









(3.12) 



performs a unitary transformation, provided 9u = for every u £ M. Furthermore, 
for IWQC procedures *P arising from the DKP or (simplified) RBB constructions, the 
geometry (G, /, O) underlying *p (together with the set M of qubits with measurement 
angle "/a ) has a modified fiow. 



"A similar generalization of the flow conditions was anticipated in ll21 . but does explicitly specify com- 
binatorial conditions for this generalization. 
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Proof. We may obtain a procedure as above using a map <& extending $, as follows. 

Let B = {J{0)^ ^Z, Zzzd}g^2Lz d>2' where we define the d qubit operator 



d times 

We may define the IWQC procedure 333'* analogously to 33, as follows: 



333-ui,...,i 




(3.13) 



(3.14) 



by an analysis deferred to Section B. 1. 1 (see also Section 3.8 of 1 19 1), we may show 

5d5tu...,vAPs) = [^^d],^,...^,_^Ps [Zzzdll^^ ^^^^ , (3.15) 

for a system S including the qubits Vj, but excluding a. For a stable-index expression 
for gates in B, we may then define: 

(3.16a) 
(3.16b) 
(3.16c) 

<^>( Zzzjvu....v,]] ^ y,^£ (3.16d) 

extending this to unitary circuits over B in the usual manner For a star circuit = 
AZ[v.w,] ■■■AZ[v,w^] Ji^u)[l] as in ([TT), we then have 



i+)[^] ; 


) = N„ ; 












333iii,...,t)d ) 



u=i 



(3.17) 



after commuting the entangling maps to the right, where G* is a star graph with center 
V and edges from v to each Wj and to u; we may define a similar star graph for the 
entangUng maps of ■■■ vd^° ^^^^ write 



<i>{Zzzd) = 




(3.18) 



for such a star-graph. 

Given a tuple (G, /, O, M) with a modified flow (/, =4), define the circuits G„ for 
u G O"^ via the following stable-index expressions: 



\ iu~/(m) / 



/(") 



foru^/(u); (3.19a) 
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C„ = (^2Zd)^^_ , for u = /(u) with neighbors wi, ... ,'u;„. (3.19b) 



We also define a circuit Cj by 



(3.20) 



uveE{G) 
u,v<fimg{f) 



we may describe the operations in this circuit by the induced subgraph G[I], where 
/ = V{G) \ img(/). Consider the circuit C described by the stable-index expression 
consisting of the concatenation of these circuits. Placing the circuit Cj at the beginning, 
and ordering the the remaining circuits C„ in a linear order ^ extending we obtain 
an expression in which operators with index u ^ f{u) occur no later than any operator 
involving f{u) or w ^ f{u), and any operator with index u = f{u) occurs before the 
indices it; ~ u are deprecated. Thus, applying the homomorphism $ to each sub-circuit 
C„ in C yields a well-formed 1 WQC procedure 



n 

ueo" 



n z 



\ ( 



n 



1 



J 



-G[i] 



(3.21) 



where Gu is the star graph for each qubit u as described above, and the congruence 
arises from commutation of preparation maps and entangling maps to the right of each 
term in the ordered product. 

Note that for every edge vw e E{G), we have either v,w ^ img(/), or one of the 
vertices (without loss of generality, v) is of the form f{u) for some u £ 0°. In the 
former case, vw corresponds to a /\Z operation in Cj, and thus to an entangling oper- 
ation in $(Cj); in the latter case, we have w ^ /(«), in which case vw corresponds 
to an entangling operation in the star-graph G„ and thus to an entangling operation in 
$(C„). Conversely, every entangling operation in *p above corresponds to some edge 
in G, by construction; the geometry underlying *p is then {G, I, O). Note also that the 
individual operations in *p and the procedure described in ( |3.12| i are identical, except 
that corrections of the form Z^'" do not occur in *P; commuting all the preparation 
and entangling maps to the right of the expression, these are induced by commuting op- 
erations Ef^u)j^{u) past the corrections X^'"^)- Thus, p.l2| i is the result of performing 
this commutation to *P; the IWQC procedure in ( |3.12| i thus performs the same unitary 
transformation as C. 

Finally, we may show that every geometry underlying a IWQC procedure obtained 
from the DKP or (simplified) RBB constructions has a modified flow. Consider a 
stable-index expression C for a circuit over the gate-set Bdkp U ,Brbb, whose terms 
of C are ordered according to the order of performance in the circuit. Let / be the 
extended successor function of C, and ^ be the order in which the indices are dep- 
recated in C. Then / will be a mapping from the deprecated indices to the advanced 
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indices, and (/, ^) satisfies ( |2.14| l by construction. If {G,I,0) is the geometry un- 
derlying the procedure <&(C), and M is the set of qubits with measurement angle "/a, 
by consti-uction (/, ^) will be a modified flow for (G, /, O, M); as the DKP and RBB 
constructions both involve producing such a stable-index expression C, and the process 
of normalizing a IWQC procedure leaves the geometry invariant, the result holds. ■ 

The foregoing proof not only demonstrates that a IWQC procedure as in ( |3.12[ l 
performs a unitary transformation, but gives a description of a unitary circuit perform- 
ing the same transformation. This will form the basis of a similar construction in 
[Section 34] for constructing circuits without reference to corrections or measurement 
dependencies. 

The proof above also describes a construction of IWQC procedures based on a 
gate-set which subsumes ,Bdkp U Srbb, containing operations acting on arbitrarily 
many qubits. The extension to this gate-set B from Sdkp and Srbb corresponds to the 
way in which modified flows extend the pattern of dependencies found in either the 
DKP or RBB constructions without normalization; the inclusion of the operators Zzzd 
for d > 2 is essentially a byproduct of the lack of further constraints on the definition 
of modified flows. 

3.2.4 The natural pre-order, and uniqueness in the case |/| = \0\ 

Following flSl, for any geometry (with or without a modified flow) and for any function 
/ : 0° — > P, we may consider an ordering satisfying the conditions of p.9| l or 
as follows: 

Definition 3.4. Let (G, /, O) be a geometry and / : 0° — > P. The natural pre-order 
=4 for / is the reflexive and transitive closure of the conditions 

V 4 f{v) (3.22a) 
w ^ f{v) =4> V ^ w (3.22b) 

for all w G O" and w e V{G). 

Note that the conditions of p.22[ ) are exactly the conditions p.llc[ ) and p.lld| i for 
modified flows. The definition of the natural pre-order then relaxes the condition that 
^ be a partial order, i.e. that =^ is antisymmetric. 

If the natural pre-order =<; happens also to be a partial order, then (/, is a mod- 
ified flow. By the definition of the transitive closure, any pre-order satisfying the con- 
ditions of p.22| i must also contain the relations of the natural pre-order =<;; then if ^ 
is not antisymmetric, there can be no antisymmetric relation ^ satisfying the modified 
flow conditions with /. 

Because of the distinguished nature of the natural pre-order, we may prove a useful 
uniqueness result for modified flows: 

Theorem 3.4. Let (G, /, O) be a geometry and M C O'^. If \I\ = \0\, then there is 
at most one function f such that (/, is a modified flow for (G, /, O, M), for =^ the 
natural pre-order of f. 
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Proof. Suppose that |/| = \0\. If / = O, then a tuple (/, ^) can only be a flow if 
f{v) = V for every v G O^, as otherwise we obtain cycles v ~< f{v) -<■■■-< f^{v) = 
V for some v. Then the only function / for which (/, may be a modified flow is the 
identity function on O'^. 

Otherwise, suppose that I ^ O, and consider two functions f,g: 0° — > P such 
that f{v) = V only for v G M, and similarly for g. Note that the natural pre-order =^ 
of / can only be anti-symmetric if there are no cycles of the form v -< f{v) 
f''{v) = v; thus, we require that / have no such cycles for A; > 1, and again similarly 
for 5. 

Suppose / and g differ, and let uq G 0° be some vertex for which /(?io) 7^ .9(mo)- 
Let g : P — > O'^ be the inverse function of g: we may recursively construct an 
unbounded sequence of vertices Uj for j e N by defining m = f{vj), and 



for j ^ 1. Note that dom{g) = P = img(/) and dom(/) = O*^ = img(g); then this 
sequence is well defined for all j ^ 0. We may also show that u^+i 7^ uj^i for each 
J ^ 1 by induction, as follows. Suppose that Wj+i 7^ Uj-i for some particular j ^ 1: 

• Suppose that Wj+i = g{uj) ^ .f{uj)' it follows that uj = f{uj-\), so that 
Uj+2 = f{uj+i) ^ f{uj-i) = Uj by the injectivity of /. 

• Suppose that Mj+i = f{uj) ^ g{uj): it follows that = g{uj-i) ^ f{uj-i), 
in which case we have Uj+2 = 9{uj-i) = Uj, by the injectivity of g. 

• Suppose that Wj+i = f{uj) = g{uj). Because / is injective, we have Wj+i = 
f{uj) Uj if Uj 7^ Uj-i; and as 7^ Uj^i, at least one of the inequalities 
Uj+i Uj or Uj 7^ Uj-i must hold. Thus, we have Uj+i — g{uj) 7^ Uj in any 
case; by the injectivity of g, we then have Uj+2 = g{uj+i) = g{g{uj)) ^ Uj by 
the acyclic property of g. 

It then follows by induction that «j+2 7^ Uj for all j ^ 0. 

We may construct a subsequence 1x^(0) -< "^(i) ~< ■ ■ ■ , with a : N — > N a 
monotonicaUy increasing function, as follows. We let cr(0) = and for each j ^ 0, 
and define 



for all j > 0. Suppose for some j > 1 that ^^(jj+i = f{u„(^j)). 

• If Ua(i)+2 = as weU, we have u„(j) < UaO)+i = ■"^o+i). and we 
also have = /(^^(j+i)). 

• Otherwise, we have ^/^(j+i) = '"(t(j)+2 = 9{'^a{j)+i)7 which implies that either 
UaU)+2 = = f{uaU)) or Ua{j)+2 ~ = f{uaU))- In either 
case, we again have UaQ) -< ^^(j+i), and once more we have = 




(3.23) 



(t(j) + 1 , if Wa(i)+2 = + = /^KO)) 

a{j) + 2 , otherwise 



(3.24) 
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As Ma(o)+i = ui = f{uo), by induction we then have ^^(A;) ^ "cr(fe+i) for all k ^ j. 
As V^(G') is finite, there must be integers N, N' such that N < N' and u„(n) = Wcr(Af')" 
then we have 'u^(^) =^ w<x(Ar+i) =^ • • • =^ "^^(Ar') = Uo-(7V), where ^^.(Ar) 7^ Uo-(7V+i) as 
noted above. Then ^ is not antisymmetric. 

If on the contrary (/, ^) is a modified flow, and {g, ^) is also a modified flow for 
some partial order ^, it follows that / = g. Thus (G, /, O) and M have at most one 
modified flow (/, =<;), for =<; the natural pre-order of /. ■ 

The uniqueness of a modified flow (/, =^), when we require that =^ be the natural 
pre-order of /, will play an important role in describing a semantic map for the DKP 
and RBB constructions. Different "flow-functions" / correspond to different sets of 
dependencies in a IWQC procedure: provided a procedure *P in which some of the 
dependencies have been obfuscated or eliminated altogether by Pauli simplifications 
and signal shifting, it may be difficult to determine whether there exists a successor 
function which describes the actual dependencies which are present. This problem is 
simplified if we are promised that there is at most one successor function / yielding a 
modified flow (/, ^): we may simply attempt to construct /, and determine whether 
the dependendices in *p are consistent with those described by /. 



3.3 Decompositions of geometries 

The proof of jLemma 3.3| suggests an approach to determining when a tuple (G, /, O, M) 
has a modified flow, by means of a decomposition of the geometry (G, /, O) into sub- 
geometries. Examining such decompositions will lead to an efficient algorithm for 
obtaining the candidate successor function / associated with a modified flow, which 
we may use to implement a semantic map for the DKP and RBB constructions. 



3.3.1 Star geometries and star decompositions 

For a geometry (G, /, O), note that the natural pre-order of a injective function / : 
— > P is characterized entirely by the transitive closure of local relations of ver- 
tices in G. Similar local relations are used in the proof of |Lemma 3.3| to define star 
graphs which are used to define both star circuits, and their representations as proce- 
dures in the IWQC model. For an injective function / : 0° — > P, we may character- 
ize these local relations in terms of the following structures: 

Definition 3.5. Let (G, /, O) be a geometry, and / : 0° — > P an injective function. 
For V e V{G) define the vertex-set 

= {w e V{G) \w(.{v, f{v)} or [w ^ f{v) and w ^ /(/(«))] } . (3.25) 

We then define the star geometry rooted at v (or centered at f{v)) is the geometry 
) given by 

G,^ = G[V,,] - {xw\x,w^ fiv)} , (3.26a) 
Iv. = \ Uiv)} , (3.26b) 
O.* = \ {v} . (3.26c) 
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Star geometries correspond to star circuits (as described in Section 3.1.1 1 and Zz 
operations, via their representations as 1 WQC procedures. Consider the 1 WQC proce- 
dure <p = $(C^[ «!,...,«„ —1 ) obtained from a stable-index expression for star circuit 

as in p.l| i. Then *p has a geometry (G, /, O) consisting of a star graph G with 
center w E O \ I, and where v G I \ O. If / is the index successor function of 
the stable-index expression, we may then show that (G, /, O) is the star geometry with 
root V and center w as in Definition 3.5. We may similarly show that ^{Zz) also gives 
rise to a star geometry, where the mediating qubit a = /(a) is both the center and the 
root of the star geometry. (The same also holds for ^{Zzzd), for the mapping $ and 



unitary Zzzj_ defined in the proof of Lemma 3.3 ) 



Star geometries (G„*, 0„*) are thus precisely those geometries which under- 
lie the 1 WQC procedures <i>(G*), where G* are the featured circuits in decompositions 
described in Section 3.1.1| Extending the results of lil2i . we may consider how such 



structures for a geometry (G, /, O) may be used to describe "candidate" circuits G for 
IWQC procedures *p, where G may correspond to *p via the DKP or RBB construc- 
tions: we consider how to construct such circuits in [Section 3.4 Motivated by the 



connection with circuit constructions, we may consider a decomposition of geometries 



similar to the compositions of simple circuits in Section 3.1.1 



Lemma 3.5. Let (G, /, O) be a geometry and M C O^. Then (G, /, O, M) has a mod- 
ifiedflow (/, ^) if and only if there is a decomposition of{G, I, O) into star geometries 
rooted at each of the vertices v € O^, together with the geometry [G [/] ,1,1^, where 
i = V{G) \ img(/). 

Proof. We proceed by induction on the size of img(/). Ifimg(/) = 0, thendom(/) = 
as well, in which case / C V{G) O; we then have (0,1,0) = (G[/],/,/). 
The natural pre-order in this case is the equality relation, which is antisymmetric; then 
(/, =) is a modified flow. The Lemma then holds for img(/) = 0. 

For the induction step, suppose that |img(/)| = iV + 1 for some N E N, and that 
the claim holds for geometries (G, /, O) with flows (/', ^) in which |img(/')| ^ N. 

• Suppose that (G, I, O, M) has a modified flow: in particular, the natural pre- 
order =<; is a partial order. Because V{0) is finite, there is a vertex v € V{0) 
such that V ^ w only for w maximal in the natural pre-order ^. Let Qy^, = 
(G,. , , O,.), and g' = (G', /, O'), where 

G' = G \ /(w) , and (3.27a) 

0'=(n./^fru ^^"^^("^V (3.27b) 
[ OA {w, /(«)}, Otherwise J 

(Again, A^B denotes the symmetric difference of A and B.) Note that Oy^, C 
O by definition, and either = A {v, f{v)} C O A {v, f{v)} (if v and 
f{v) differ) or = C O (if f{v) = v). Because ¥(0,,) = U {f{v)} 
in both cases, we then have V{0') n V{Ov*) C n O' : then the composition 



47 



(G,/,0) 
G 



Gv* ° G' is well-defined, and we have 

U G' 

{xw€G\w = f{v)} U [g \ f{v) 



G; 



u / 



o 



0„ U (O' \ = 0,„, U (O \ 0„.) 



(3.28a) 



(3.28b) 



(3.28c) 



Then, we have (G, /, O) — Gv* ° G' ■ 

The restriction /' of the flow-function / to V{G) \ O' does not include f{v) in 
its image, and therefore has cardinality at most N ; by recursively applying the 
decomposition above, G' has a decomposition into star geometries and the geom- 
etry (G' \ img(/') , /, V{G') \ img(/')), where in particular V{G) \ img(/') = 
V{G) \ img(/) = /, by the definition of /' . Thus, for any tuple (G, /, O, M) 
with a modified flow, {G,I,0) has a decomposition into star geometries as 
above. 

For the converse, suppose that {G,I,0) ~ Gv* ° G' for some star geometry 
Gv* = {Gvi,, Iv*,Ov*) and a residual geometry C?' = {G',I,0'). If the latter has 
a star decomposition, and if /' is the restriction of / to V{G'), then (/', is a 
modified flow for G' by hypothesis as img(/') = img(/) \ {f{v)}. By assump- 
tion, Gv* and G' are composable, in which case we have V{Gv*) H V{Gv*) = 
Iv* n O'. If V ^ w for some w e V{G), we then have either w = f{v), 
w ^ O' w, or w w' for either w' = f{v) or w' ^ f{v)- By construction, we 
have either O = O' w in the case that v — f{v), or O — {O' \ {v}) U {/(«)} 
otherwise; then for any w' = f{v) or w' ^ f{v), we have w' £ O. As such w' 
lie outside of the domain of /, they are maximal in the natural pre-order; then 
w > V only if w e O. In particular, if the restriction of the natural pre-order 
to G' is such that (/', =4) is a modified flow, then (/, =<;) is also a modified flow. 
Thus, the existence of a star decomposition for (G, /, O) implies that it has a 
modified flow, by induction. ■ 



Thus, we may decompose a circuit with a modified flow into star geometries. Fig 



ure 10 illustrates a decomposition of a geometry with a flow into star geometries in 



the manner described above; Figure 1 1 illustrates a more general case for a modified 
flow. Also illustrated in these figures are "candidate circuit" constructions for unitary 
circuits, which we describe in |Section 3.4| 



3.3.2 Structural characterization of star decompositions 

In order to define a semantic map S for the DKP and RBB constructions, it will be 
necessary to attempt to obtain a star decomposition of a geometry which is initially 
not known to have a modified flow. We may facilitate this by showing that we may 
characterize star decompositions independently of modified flows, as follows; 
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Figure 10: Illustration of a decomposition of a geometry {G, I, O) with a flow (/, 
into star geometries for each v e 0°, and a reduced geometry (G[/],/, /) for 
I = y(G)\img(/). In the upper three panels, the action of / is represented by arrows; 
different styles of edges indicate edges belonging to different star geometries. The mid- 
dle panel illustrates each star geometry in the decomposition separately (with vertices 
repeated between different star geometries), adjacent to the geometries with which they 
may be composed. The bottom panel illustrates the unitary circuits corresponding to 
these star geometries (including a labelling of wire-segments corresponding to stable 
tensor indices), and the complete circuit obtained from composing them. 
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Figure 1 1 : Illustration of a decomposition of a geometry {G, I, O) with a modified 
flow (/, into star geometries for each v G 0°, and a reduced geometry (G [/] , /, /) 
for I = V{G) \ img(/). In the upper three panels, the action of / is represented by ar- 
rows. Different styles of edges indicate edges belonging to different star geometries (as 
are edges of the same style, but which are disconnected from each other). Arrows are 
omitted for f{v) = v. The middle panel illustrates each star geometry in the decom- 
position separately (with vertices repeated between different star geometries), adjacent 
to the geometries with which they may be composed. The bottom panel illustrates the 
composition of the unitary circuits corresponding to the star geometries (including a 
labelling of wire-segments corresponding to stable tensor indices). 
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Lemma 3.6. Let (G, /, O) be a geometry and M C O'^. Suppose that {G, I, O) can 
be expressed as Qh ° Q(i-i)* o ■ ■ ■ o Gi-k ° Qj, where Gj — /, I) for some subset 

i C V{G) with \i\ \0\ = \V{G)\ - £, and where for each j < i the geometry 
Qji, = {Gj, Ij, Oj) satisfies either 

(0 Ij = Oj = V{Gj) \ {vj}, for some Vj £ V(G) D M which is adjacent to every 
other element ofV{Gj); 

(ii) V{G) = Ij U Oj, Ij = {Ij n Oj) U {vj}, and Oj = {Ij n Oj) U {wj},for 
some distinct Vj, Wj G V{G) such that Wj is adjacent to every other element of 
V{Gj). 

Let f : O'^ — > I'^ be defined so that f{vj) = Vj for any Vj as in the first case above, 
and f{vj) = Wj for any Vj , Wj as in the second case; and let =^ be the natural pre-order 
for f. Then (/, is a modified fiow for (G, /, O, M). 

Proof. Let / be defined as above, and let =^ be its natural pre-order. By construction, 



each geometry Qj is then a star geometry with respect to / in the sense of Definition 3.5[ 
whether or not (/, is a modified flow. We may show by induction on ^ = | V^(G) | \ 
\0\ that =^ is a partial order 

• If £ = 0, then O = V{G), in which case dom(/) = 0. Then the natural 
pre-order is the equality relation, which is antisymmetric and therefore a partial 
order. 

• Suppose the Lemma holds for all geometries with fewer vertices than (G, /, O), 
but with an output subsystem of the same size. By hypothesis, we have (G, /, O) - 
{Gi, Ii, Oe) o Q', where we may write 

g' = (G', /, O') = Gi^t-l). o • • • o o ; (3.29) 

in particular, we have V{G') n V{Gi) C n O' . If /' is the restriction of / 
to O' , and ^' is the restriction of =<; to V{G), the pair (/', =<;') is then a modified 
flow for Q' . 

(0 Suppose f{vt) — vg: then V{G) \ V{G') — {vi}, and we have w e O' = 
O for every w ~ vg. Then the only relations that =<; contains that =<;' does 
not are u ^ vi for qubits such that f{u) £ O, and vg ^ w for w E O; that 
is, u ^ Vi ^ w for some w maximal and u bounded above in only by 
maximal elements. 

(//) Suppose f{ve) ^ wg ^ vf. then ve E O' \ Oi, and G \ O' . Then 
the only relations that =<; contains that does not are u =<; for qubits 
such that f{u) — w such that w ~ we, and vg ^ w for qubits w ^ wg. 
In either case, we have wi £ O'; thus we only add additional (maximal) 
upper bounds for u bounded only bounded by maximal elements, and new 
upper bounds for V( which are also maximal. 

In either case, =<; is a partial order; thus (/, =<;) is a modified flow. ■ 
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We will call a decomposition as in Lemma 3.6 a star decomposition of a geome- 



try (G, /, O), and the associated geometries Qj in the decomposition star geometries, 
whether or not they are defined with respect to a corresponding modified flow (/, =^). 
More generally, we will call a geometry (G^,, J*, O^,) a star geometry when the input 
and output subsystems satisfies either property (0 or property (//) in the statement of 



Lemma 3.6 above. In the case that = O*, we say that unique vertex v G V{G^) \ 
is both the root and the center of (G*, O*); otherwise, we say that the unique vertex 
u e /* \ is the root of (G*, O*), and the unique vertex w € \ is the 



center. 



The proof of Lemma 3.6 suggests a recursive decomposition of a geometry (G, /, O) 



into star geometries using "maximal" star geometries, in the following sense: 

Definition 3.6. For a geometry (G, /, O), a maximal star geometry is a star geometry 
(G^, O,,) such that (G, /, O) = (G*, I^,, O*) o (G', /, O') for some sub-geometry 
(G'j /, O'). We call the latter the residual geometry of this decomposition. 

We may characterize a maximal star geometry as follows: 

Lemma 3.7. Let (G^,, I^,, O^,) be a maximal star geometry, with root v and center w, 
for some geometry (G, /, O). Then one of the following two properties hold: 

(J) V ^ w ^ C] O'^, and is adjacent in G only to elements ofO; 

(ii) V ^ w, w Cz O, and v G O'^ is the only neighbor of w in G which is not in O. 

Conversely, for a geometry (G, /, O) and vertices w, w G V{G) for which either of the 
above hold, the geometry {Gyi,, ly^,, Ot,*) such that Gyi, = G[W]for W consisting of 
w and its neighbors in G, with input and output subsystems = V{Gyi,) \ {w} and 
Oy-k = V{Gyi,) \ {«}, is a maximal star geometry for (G, /, O). 

Proof. We first show that any maximal star geometry has either property (/) or prop- 
erty (;;). By the definition of a maximal star geometry, we have C O; as v ^ O^, 
we have v G C^. If v = w, by definition v is adjacent to every element of — O-i,, 
and every neighbor of i; in G is an element of thus v is adjacent only to elements of 
O. Also, in the case that w = w, we have v ^ I^^hy definition, in which case v G P 
as well. Otherwise, we have w G 0^, C O and w G \ O*, so that v E O^. Also, 
if w 7^ w by definition, all edges incident to w in G are of the form uw, where either 
It = u or tt G C O; then v is the unique element of 0° to which w is adjacent in G. 

Conversely, suppose that v and w are vertices satisfying either property (/) or prop- 
erty (ii), and let (Gt,*, Oi,*) be the geometry described. Let G' — G \ w and O' = 
O \ {w}. In either case, we may show that (G, /, O) ~ {Gy^,, ly^,, Oyi,) o (G', /, O'), 
and that (Gt,*, Oyi,) is also a star geometry. The Lemma then holds. ■ 

By definition, any geometry with a star decomposition has a maximal star geom- 
etry; we may then attempt to obtain a star decomposition by recursively constructing 
maximal star geometries for successive residual geometries. In the next section, we 
show how this leads to an efficient algorithm for obtaining a star decomposition (and 
a modified flow) for an arbitrary geometry (G, /, O) and subset M C O^, when one 
exists. 
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3.3.3 Efficiently decomposing geometries 

In order to be able to test the consistency of measurement dependencies of a IWQC 
procedure *p with those arising in the DKP or (simplified) RBB construction, we may 
attempt to determine whether the tuple (G, /, O, Af ) associated to Cp has a modified 
flow. As a special case, we may efficiently determine when a geometry (G, /, O) has 
a flow, and construct one if one exists. The most efficient known algorithm to do this 
is that of Mhalla and Perdrix |14|. We show how their algorithm may be adapted to 
produce modified flows when they are present, using an analysis in terms of maximal 
star geometries. 

Maximally delayed fiows. We may summarize the algorithm of lfT4ll as follows. To 
find a flow (/, or a geometry (G, /, O), we consider the length of the longest chains 
of each qubit from the output set, and construct "layer sets" consisting of vertices of 
similar depth: 

Definition 3.7 (IHl, Definitions 4 & 5). For a geometry (G, /, O) and a flow (/, s^), 
let Vf^{G) be the set of vertices of V{G) which are maximal in ^, and let V^+i (G) be 
the maximal set of vertices of the set 

W< = ViG) \ V<^ (3.30) 

in the order ^, for any j > 0. A flow (/, =4) is more delayed than another flow (/', ^) 
if we have | | ^ | | for all j € N , where this inequality is strict for at least one 
such j ; and (/, ^) is maximally delayed if no flow for (G, J, O) is more delayed. 

A maximally delayed flow is also a flow of minimum depth, i.e. that the order ^ 
of such a flow has the same depth as the natural pre-order for the given flow-function 
/. For a maximally delayed flow (/, =^), let £ : V{G) — > N be the function mapping 
each vertex v E V{G) to the integer j such that v G . Then, the partial order ^ 
characterized hy v ^ w [v — w qj: £{v) > £{w)] also yields a minimum-depth 
flow by construction. 

Mhalla and Perdrix show for a maximally delayed flow (/, =4) for (G, /, O) that 
Vf^ = O, and is the set of elements of such that, for some w £ O, v is the 
only neighbor of w which is contained in . By induction, augmenting the set O to 
include successive "layers" Vj^, one may show that 

3w e V{G) \ Wr< : Ndw) n VK/ = , (3.31) 

where Ng{w) is the set of vertices adjacent to w in G. By a depth-first search from O, 
we may find in a finite number of iterations of j those qubits w £ V{G) \ which 
have such a neighbor v, and accumulate them into new layers V^^i- 

Relation to maximal star geometries. This decomposition of the geometry (G, /, O) 
into layers V^^ corresponds strongly to the star decomposition of a geometry (G, /, O). 



^1 = [vew. 
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The maximal elements of V{G) in the natural pre-order are exactly the elements of O 
(as z =4 f{z) for any z E 0°): thus, the root w of a maximal star geometry in (G, /, O) 
is one such that all neighbors z ~ f{v) are elements of O or equal to v itself. The 
set as described above is the set of precisely those vertices which are the root of a 
maximal geometry. Thus, the algorithm of 1 14] may be described as finding flows by 
performing a star decomposition of {G,I,0) when it has a flow. 

We may easily extend the maximal delay principle of jl4| for flows to one which 
identifies a maximal star geometry for (G, /, O) in the more general case where this 
geometry (together with a set M) has a modified flow. For such a tripl e (G, /, O, M) 
with a modified flow (/, =^), we may define the sets Vj^ and as in 



Definition 3.7 



Then for each j e N, V^+i consists of those vertices v £ Wj such that, for any 
w G V{G) such that w = f{v) or w ~ f{v), we have w — v or w £ V{G) \ W, 



-< 

For V ^ M , this implies that (as in the analysis of 11141 ) the only neighbor of f{v) 
which Ues in is v itself; for v e M, it may also be possible that v has no neighbors 
in WJ^. The flow-finding algorithm of [14J can then be easily extended to test for 
the latter condition, yielding a procedure to produce a star decomposition of (G, /, O) 
when one exists, as we describe below. 



Bounding the search for roots of maximal star geometries. The flow-finding algo- 
rithm of Mhalla and Perdrix ([14 |, Algorithm 1) operates on the basis of maintaining a 
set of "correcting vertices" w which, for each new layer V^+i, may potentially be the 
successor of some vertex v e V^"^!. These correcting vertices are precisely the centers 
of maximal geometries, in which case the corresponding u is a root. We thus take a 
similar approach to obtain a star decomposition by "removing" layers of maximal star 
geometries, identified with their centers. 

For modified flows, we must accommodate vertices v E M which may be simul- 
taneously the root and the center of a geometry. To ensure that we do not fruitlessly 
examine vertices of M at each layer, we maintain a digraph whose arcs point in either 
direction along each edge of G, and from which we delete any arc ending at a vertex 
which has already been indicated as the root of a star geometry. Then, if a potential 
center w has a single neighbor v which has not been matched to a star geometry, v is 
the root and w the center of a maximal star geometry; and if a vertex w has not yet been 
identified as the center of a star geometry only neighbors vertices which have already 
been identified as centers, w is both the center and the root of a maximal star geometry. 
As a result, we may easily determine whether a potential center w has a single neighbor 
V which has not been matched to a star geometry (in which case v is the root and w 
the center of a maximal star geometry), or whether a vertex which has not yet been 
identified as the center of a star geometry only neighbors vertices which have already 
been identified as centers (in which case w is both the center and the root of a maximal 
star geometry). 

[Algorithm 2| defines a subroutine RemoveStarGeom to perform the necessary 
arc deletions for the digraph described above whenever we define the layer of a vertex, 
as well as several other operations useful to an iterative procedure for discovering the 
maximal star geometries of a star decomposition. 
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Data : M, a subset of specially designated vertices 

W, the set of vertices not yet marked as roots of star geometries 
emptyLayer, a flag indicating if the present layer is empty 
D, a digraph in which vertices v e V{G) \ W have in-degree 
/ : V{G) — ^ y{G) a mapping from "roots" to "centers" 
S' , a set of possible centers of star geometries for the next layer 
V E V{G), a vertex to mark as the root of a star geometry 
w e V{G), a vertex to mark as the center associated to v 

Effect: Assigns f{v) <— w, adds v to the set of potential centers for the next 
layer and to the ordered Ust of qubits with identified star geometries, 
removes all of the inbound-arcs to v in the digraph D, and checks the 
neighbors of w in £> for vertices in M with out-degree to mark as a 
center for the next layer. 

/ subroutine RemoveStarGeom('i;, w) : 

2 Set emptyLayer <— false. 

3 Remove v from the set W, and set f{v) ^ w. 

4 Insert v to the beginning of L. 

5 For every neighbor z with an outbound arc to v in D, do: 

6 Remove the outbound arc z ^ u from D. 

7 It z Cz M n W and z has no more outbound arcs, then add z into 5'. 

8 endfor 



Algorithm 2. A subroutine to attribute a vertex to the current layer, for use in an 
iterative procedure for finding the roots of maximal sub-geometries with a modified 
flow. In particular, this subroutine maintains the arcs of a digraph D (a "sub-digraph" 
of the graph G of a geometry, interpreted as a symmetric digraph) from which we 
remove any arc inbound to a vertex with a well-defined layer. 

Accumulating roots of star geometries. Using RemoveStarGeom as a subrou- 
tine, we may easily define a procedure similar to Algorithm 1 of lfT4ll which either 
(a) determines that a star decomposition of the geometry does not exist, implying the 
non-existence of a modified flow, or (b) constructs a modified flow, represented as a 
tuple (/, L) consisting of a partial function / : V{G) — ^ V{G) mapping the root of 
each star geometry to the corresponding center, and a list L of the roots corresponding 
to a composition order for the star geometries f [Algorithm 3 [ presents such a procedure, 
which we describe below. 

We begin by creating the digraph D which governs the search for roots: we copy 
all of the neighborhood relations of G, omitting those arcs ending in O, as these cannot 
be roots of star geometries. Similarly, as no w G O is a root, it must lie outside the 

*The list L constructed by|Ai gorithm 3] is a linear order, i.e. a partial order of maximum depth. To find 
one of minimum depth as in Algorithm 1 of |14| , it suffices to maintain a layer counter and to record the 
layer for which each vertex v £ O*^ is the root of a maximal star geometry. However, as we are ultimately 
interested in producing a (linearly ordered) circuit expression, it is not necessary to obtain a partial order of 
smaller depth. 
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Input : G, a graph with subsets /, O, M, where M C V{G) \ O consists 






entirely of vertices with default measurement angles Va. 




Output: Either nil, if (G, /, O, M) does not have a modified flow; or an 






ordered pair (/, L) consisting of a partial function 






/ : V{G) ~^ V{G) mapping roots of star geometries to their 






respective centers, and a sequence L of the roots representing a 






composition order for star geometries. 


; 


procedure FindModStarDecomp(G', /, O, M) : 


2 


Let / : V{G) — ^ V{G) be a mapping from "roots" to "centers". 


3 


Let L be a sequence of vertices (initially empty). 


4 


Let ^ be a set of potential roots of star geometries. 


5 


Let ^ be a set of potential centers of maximal star geometries. 


6 


Let D be the completely disconnected digraph on V{G). 


7 


For all w e V{G),do: 


8 




For every w e Ng(v) \ O: add the arc v ^ w to D . 


9 




If w e 0,then 


10 




Set f{v) <- nil. 


11 




If t; ^ /, then add v into S. 


12 




else 


13 




Add V into W. 


14 




endif 


15 


endfor 


16 


Let emptyLayer be a flag indicating failure to find any new roots. 


17 


Repeat: 


18 




Set emptyLayer ^ true, initially. 


19 




Let S" ^ be the set of potential centers for the next iteration. 


20 




For each w e S, do: 


21 




If w e M^,then 


22 




Call RemoveStarGeom(z(;, w). 


23 




else, if D has only one outbound arc w — > w, then 


24 




Call RemoveStarGeom(?;, w). 


25 




If u ^ /, then add v into S' . 


26 




else 


27 




Add w back into S' . 


28 




endif 


29 




endfor 


30 




Set S ^ S' for the next iteration. 


31 


until emptyLayer . 


32 


liW = 0, then return (/, L); else return nil. 



Algorithm 3. An iterative algorithm to discover a star decomposition inspired by 
the analysis of 1141 . using the subroutine RemoveStarGeom of |Algorithm 2| 
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domain of /, so we set f{w) ^ nil; but if w £ O is not an element of I, it will be 
the center of some star geometry, so we insert it to a set S of candidates for centers of 
maximal star geometries. Any vertex w ^ O is the root of some star geometry in a star 
decomposition of (G, /, O), provided such a decomposition exists; we then insert such 
w into a set W of candidates for roots of maximal star geometries. 

We then search for layers of new roots, corresponding to the layer-sets in the anal- 
ysis of 114 1 . For each layer, we iterate through the set of potential centers for maximal 
star geometries, which initially includes no elements of W. We maintain the invariant 
that D contains all (directed) adjacency relations of G, except that it contains no arcs 
ending at vertices w ^ W. We also maintain the invariant that v G W becomes an 
element of S only if it is an element of M and has out-degree zero in D, in which case 
it is both the root and the center of a maximal star geometry. We indicate the removal 
of vertices (as a result of decomposing the geometry into star geometries and an as-yet 
decomposed "residue") by absence from both the sets S and W; thus S represents at 
each stage the output subsystem of the residual geometry. 

In any iteration, for the center w of any maximal star geometry, either w has a 
unique neighbor v ^ W in G (in which case w ^ v is the only outbound arc from w 
in D), OT w G M and has no neighbors in W (in which case it has out-degree in D). 

• In the latter case, unless w is isolated in G, there was some final neighbor u which 
was removed from W in an earlier iteration; immediately after the removal of u 
from W, we may then identify that w is the root (and center) of a maximal star 
geometry. In the call to RemoveStarGeom(M, m'), we thus examine each of 
the neighbors w G M D W (restricting our search to qubits in AI which have 
not yet been assigned as the root to a previously removed star geometry) which 
have an outbound arc w ^ u, to see if that arc is the last outbound arc from w. 
If this is the case for some w, we add w into the set of potential centers S' for 
the next iteration. 



In Algorithm 3 this is the only condition under which an element of W may be 
inserted into S"; the only other condition under which S' accumulates elements is 
for a vertex v which has already been removed from W. Thus, in the subsequent 
iteration, any element of 5 n Vt^ is necessarily the root (and center) of a maximal 
star geometry, which we may then decompose from the parent geometry. This 
decomposition is represented by the removal of w from W, and the omission of 
w from S". 

• In the case where w has a unique outbound arc w ^ u in D, the vertex v is nec- 
essarily the root of a maximal star geometry for which w is the center. We may 
then decompose this star geometry from the parent geometry, which we represent 
by the removal of v from W (and the insertion of v into S'), and the omission 
of w from S'. In the next iteration, as an element of the output subsystem of the 
residual geometry, v is then a potential center of a maximal star geometry. 

For any w G S which is not the center of a maximal geometry as above, we carry w 
forward for the next layer by including w in S'. We assign the list of new candidate 
successors to S before the end of the iteration; and we repeat until we cannot find any 
more roots for some layer. 
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If we cannot assign any vertices to new layers, this is either because we have ex- 
hausted the possible roots (in which case W — 0), or because there are no vertices 
which are reachable from the existing layers which have viable successors. In the latter 
case, there is no modified flow; and in the former, we return the tuple (/, L). 



Run-time analysis of Algoritlini 3. The run-time of Algorithm 3 is similar to that 
of Algorithm 1 of |14|, which we may show as follows. In the following, we let 
n — \V{G)\, k — \0\, and m = \E{G)\. We assume that the graph G and digraph 
D are represented with adjacency lists for each vertex, in the latter case using lists 
of neighbors from inbound and outbound arcs separately. The partial function / is 
represented with an array structure taking the value nil (indicating a vertex not in the 
domain) or vertices in V{G). We may represent the sets M and W as arrays of pointers 
to nodes of a linked-list structure, enabling constant-time membership checking and 
element insertion, as well as linear-time traversal of the elements of the set. The sets S 
and S' we represent simply as linked lists, which also provides us with constant-time 
initialization to the empty set. 

In the subroutine Remove St arGeom(ti, w), the run-time is dominated by the loop 
over the neighbors zofv, which requires time at most 0{dcgQ{v)) . As this subroutine 
is called at most once for each vertex v, the total run-time of Remove StarGeom(t;, w) 
over the execution of FindModStarDecomp is 0( degQ(i;)) = 0{m). 

V 

Initially, 5" has k elements, none of which are elements of W ; and in each iteration, 
as every element w E S \ W gives rise to the removal of some vertex v from W and 
the potential insertion of v into 5", there always remain at most k elements of 5 \ 
at each iteration of the loop on lines 17-31. The comparisons on lines 21 and 23 can 
both be performed in constant time, as can the set-inclusions on lines 25 and 27. The 
assignment on line 30 can be performed simply by re-attributing the element list of S' 
to S. There are at most n iterations of the loop on lines 17 - 31, as it only repeats as 
long as at least one vertex has been assigned to a new layer Then, setting aside the 
work performed for vertices w E S O W in the loop on lines 17-31 and the work 
performed by RemoveStarGeom, the work performed by that loop is 0{kn). 

As elements of W are only included into S' when they are guaranteed to be at- 
tributed to the next layer, the total work performed by that loop (again aside from work 
performed by RemoveStarGeom) for elements of M is then 0{n). The run-time of 
FindFlow(G, /, O, M) is then 0{kn + m), dominated by the work performed by 
RemoveStarGeom(v, -w) for each w £ V{G), and by the iteration of the loop on 
lines 17-31 for at most n layers. 

3.3.4 Extremal analysis of geometries with modified flows 



The run-time of O^kn+m) for Algorithm 3 is precisely the complexity of Algorithm 1 



of lfT4l . However, Mhalla and Perdrix also use an extremal result presented in fT5| 
which bounds the number of edges m in a geometry with a flow by kn — C^^^)- This 
allows any geometry with edges in excess of this to be rejected at the outset, so that 
the run-time for geometries satisfying this bound may be simplified to 0{kn). By 



observing the correspondence drawn in Section 3.1.1 between successor functions and 
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the ordering described in ( |2.14[ ) for indices of stable-index expressions, we may extend 
this extremal result as follows: 



Lemma 3.8. Let (G, /, O) be a geometry such that (G, /, O, M) has a modified flow, 
for some vertex set M C V{G). If\V{G)\ = n, \E{G)\ = m, and \0\ = k, then 
m ^ nk — C^^^)- Furthermore, there exists a geometry with a modifled flow which 
saturates this bound. 

Proof. Consider a decomposition {G,I,0) = {Gyi,, Iy^,,Oy^,) o {G',I,0') into a 
maximal geometry and a residual geometry. If u is both the root and the center of the 
geometry, then it is adjacent to at most k vertices, as it is adjacent only to elements of 
O' = O. Otherwise, let w — f{v) be the center of (G„*, Oy*)', then w is again 
adjacent only at most k vertices, as it is adjacent to v, and otherwise only to elements 
of O \ {w}. In either case, {Gy^,, lyi,, Oy*) contributes at most k edges to (G, /, O). 
In either case, we have \0'\ = \0\ = k. 

By induction, each star geometry in a decomposition of {G,I, O) contributes at 
most k edges, and has an output subsystem of size k. Then the output subsystem 
/ = V{G) \ img(/) of the initial geometry {G[I],I,I) in the star decomposition 
also contains k vertices; this geometry then has at most ((j) = ifc(fc — 1) edges. 
Furthermore, as {G[I],I,I) does not contain the center of any star geometry, and as 
each star geometry has a center distinct from the others, there are at most n — k star 
geometries in the decomposition, each with a distinct center The number of edges in 
(G, /, O) is then bounded by 

{n - k)k + ^k{k ~ 1) ^ nk-\k^-\k ^ nk-{^+^), (3.32) 

as required; this bound is saturated by the extremal construction of ifTSll . ■ 



As a result, we may simplify the run-time of Algorithm 3 by 0{kn) in this case 
as well, again by a pre-processing stage which returns nil for geometries with more 
than kn - (^X^) edges. 



3.4 Candidate circuit constructions 

Following the approach in the proof of |Lemma 3.3[ we may consider how to construct 
"candidate" circuits G for IWQC procedures *p, by decomposing of the underlying 
geometry into star geometries and mapping each star geometry to a corresponding 
unitary circuit. This will give rise to a circuit G over the gate set Sdkp or Srbb, such 
that ^{G) is another IWQC procedure with the same underlying geometry. If we 
can prove that the circuit G performs the same operation as then G describes the 
semantics of *p in the unitary circuit model. We may use this as the final step towards 
defining a semantic map S for the DKP and (simplified) RBB constructions. 



3.4.1 Star circuits corresponding to star geometries 

For a star geometry, we may consider a corresponding star circuit, possibly parameter- 
ized (in the case of a star geometry with root v and center w — f{v) ^ v) by an angle 
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related to the measurement angle 0^ of v: 



( 



\ Uy^V 



\ 



J 



(3.33) 



For star geometries whose center and root are the same (in the case that f{v) = v), cor- 
responding to the mediating qubits of a Zzzd operation for some d > 2 as in p.l3| ), the 
corresponding star circuit is just the operator Z2z^ux,...,ud^ acting on the input/output 
subsystems: this corresponds to the decomposition of Ztz^ as 



(3.34) 



generalizing p.2[ ). We may also correspond the geometry (G[/], /, /) to a circuit con- 
sisting solely of f\Z gates, as described with p.20| i: 



(3.35) 



as described in Section 3.1.1 such a circuit arises in star decompositions of circuits 
as well. We may then construct a candidate circuit C as in the preceding section by 
re-composing the circuits corresponding to the geometries in the star decomposition of 



(G, /, O). Such re-composed circuits are illustrated in the bottom panels of Figure 10 



and [Figure 11 For instance, in Figure 10 we may recover the geometry in the upper 
left-hand panel by composing 



(G,/,0) = {Gc.,Ic*,Oc.)o{Gd.,Id*,Od.)o{Ga.,Ia*,Oa.)o 



(3.36) 



(Gfc,,/b.,OfcJo (G [/],/,/) 
and similarly, the larger circuit C in [Figure 10| is obtained by composing 

C ^ Cc.o Cd* o Ca. o Gb. o G| , (3.37) 

where the circuits Cu^, are those corresponding to each star geometry, and Gj corre- 
sponds to the geometry (G [/],/,/) . 



3.4.2 The gate model for candidate circuit constructions 



Neither of the circuits described in [Figure 10| or [Figure ll| are generated over the gate 
set {H, T, AZ}; and the latter is not defined over eit her /BpKp or Brsb, but instead over 
the more general gate-set B described in the proof of Lemma 3.3 (As we suggested just 



after the proof of that Lemma 3.3, this gate-set corresponds naturally to the structure 
of tuples (G, /, O, M) with modified flows; we may refine this statement and describe 
this set of gates as that which naturally corresponds to star geometries. In order to use 
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star decompositions to describe a semantic map corresponding to representations of 
circuits over {H, T, ^Z}, we must first translate the circuits from the gate-set B. 

There is one apparent problem: the circuits arising from such a decomposition are 
more general than those in the domain of <!>. This may conceivably complicate the task 
of comparing the circuits arising from a semantic map S based on star decompositions, 
and the circuits which lie in the domain of either 7?.dkp or T^^rbb ■ We may show that 
this more general sort of construction will not arise in practise for IWQC procedures 
resulting from either the DKP construction or the RBB construction when |/| = \0\. 
In this case, a geometry has at most one star decomposition up to commuting terms, 



by the uniqueness of the modified flow (/, =<;) shown in Theorem 3.4 for this case. In 
particular: 

• for (G, /, O) which underly procedures Cp produced by the DKP construction, 
there will be no qubits v E O'^ such that f{v) — v, because (/, will be a flow; 

• for (G, /, O) which underly procedures *p produced by the RBB construction, 
the only edges vw in (G, /, O) for which v ^ f{w) and w ^ f{v) are those such 
that either v = f{v) oiw = f{w): i.e. where one of the qubits is a mediator qubit 
for a Zz operation between two logical qubits. 

Thus, while the candidate circuits for arbitrary geometries with modified flows may be 
more general, the candidate circuits for IWQC procedures *p arising from a circuit G 
over Bdkp or Srbb will be restricted to the same gate-set as G itself when |/| = \0\. 

For a semantic map defined relative to a particular construction, we may also ensure 
that we obtain candidate circuits over the appropriate gate-set by imposing restrictions 
on the type of star geometries which we permit in the decomposition of a particular 
geometry. However, as we considering only the special case |/| — \0\, this will not be 
necessary for our analysis. 



3.4.3 Relationships between IWQC procedures and candidate circuits 

Properly speaking, a candidate circuit is a construction from a IWQC procedure *P 
rather than from a star decomposition of a geometry (G, /, O), as the angular parame- 
ters of the J{9) gates in a candidate circuit are determined by the measurement angles 
in ^p. In order to use candidate circuits as a means of defining a semantic map S, we 
must consider the relationship between IWQC procedures and their candidate circuits. 

By [Theorem 3.1| we may characterize the dependencies of operations in procedures 
*P obtained from either the DKP or RBB constructions, relative to a candidate for the 
extended successor function for a corresponding circuit. If the geometry (G, /, O) 
underlying ^ (together with the set M of mediating qubits) has a modified flow (/, ^), 
this defines a candidate circuit G, in which the gates J{au) arising from star-circuits 
Cu* (for a qubit u ^ f{u) at the root of a star geometry) have their angles fixed by 
a„ = —9u, where 6^ is the default measurement angle of u in *p. 

By the observations made above, the circuit G will be defined over the gate set 
Bdkp or Srbb, depending on whether *p arises from the DKP or RBB construction: we 
may then consider the IWQC procedure *P = ^{C), which will coincide with ^{C). 
The procedure *p will have the same geometry as *p, and the same default measurement 
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angles: if it also has the same correction and measurement dependencies, this implies 
that *p performs the same transformation as the circuit C. 

For procedures *p arising from the DKP or RBB constructions in the case that |/| = 
\0\, there is only one modified flow for (G, /, O), which must then be the extended 
successor function for the circuit over the Kdkp or Srbb gate sets produced in the 
course of constructing *p. It follows then that C is also precisely the circuit arising 
in the construction of *p, in which case the normal form of $(C) will be congruent 
to Thus, to obtain a semantic map for the DKP construction when |/| = it 
suffices to translate C into the gate-set {H, T, ^Z} . For the RBB construction, to 
obtain a circuit with complexity no greater than the original circuit from which *p was 
produced, we must also simplify C by removing products of the form J(0) J(0) and 
JC'Vz) J^Va) J(V2) introduced by the construction of Section A.2.2 

In the case that |/| ^ \0\, even conditioned on *p arising from the DKP or RBB 
constructions, it may be difficult to obtain the modified flow function / corresponding 
to the circuit arising in the construction of *p, and therefore to produce a suitable candi- 
date circuit. Without an algorithm that can construct such candidate circuits in the case 
that |/| — \0\, for instance by finding the modified flow (/, for which *p satisfies 
[Theorem 3.1| we must therefore restrict ourselves to the DKP and RBB constructions 
in the case where no fresh qubits are introduced, i.e. which perform unitary bijections. 

3.4.4 Efficient construction of candidate circuits over {H, T, /\Z} 

For the sake of completeness, we present a procedure BuildCircuit (extending the 
remarks made in Section III B of |12| ) for constructing a stable-index representation 
for a candidate circuit for *p, represented in the gate set {H, T, A^} (rather than the 
sets Bdkp or Srbb)- The procedure BuildCircuit takes as input a representation 
of a star decomposition of a geometry (G, /, O): we may represent this in terms of 
the graph G of the geometry, the successor function / from the modified flow which 
defines the star decomposition, and a sequence L of the roots of the component star 



geometries. This procedure is presented in Algorithm 4 



For each root u in the list L, BuildCircuit determines whether it is the root of 
a geometry whose root and center are both u (corresponding to a Zzzd gate acting on 
the neighbors of u in G), or where the center v — f{u) differs from u (corresponding 
to a star circuit as in p.33| l. It then appends to a circuit expression G the appropriate 
sequence of operations: 

• By expanding Z = 1— 2|l)(l|in the tensor product Z ^ ■ ■ ■ ^ Z, and con- 
sidering the expansion of Zzzd — e^*^^* in terms of multi-qubit operations 
fQj. j < d, we may show that 



Zzzd[ 



Wi,...,Wd OC 



(3.38) 



We may construct such a gate for Wj ranging over all neighbors of a qubit u = 
f{u) by obtaining the set N of these neighbors, and for each w E N, adding 
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Input : G, the a graph for the geometry Q of a IWQC procedure 
/ : V{G) — ^ V{G), a modified flow function for G 
L, a list of roots of star geometries in a decomposition of G 
t : V{G) — ^ Z an array denoting multiples of ''^U 

Output: A candidate circuit generated over {H, T, ^Z}, corresponding to 
the star decomposition described by G, /, and L, with powers of T 
gates given by the parameters stored in t{u). 

1 procedure BuildCircuit(G, /, L, : 

2 Let G a stable-index representation of a circuit, initially empty. 

3 For each u G L, do: 

4 Itu= f{u) , then 

5 LetN^Nciu). 

6 For each w G N, do: 

7 For each w' G TV: add AZ[ to C. 

8 AddT[w] T[w] toC. 

9 Remove w from N. 
10 endfor 

// else 

12 Let V ^ f{u). 

13 Let N ^ Ng{u) \ {u, f{v)}. 

14 Foreach w G N: add A^[w,-!«] to G. 

15 Addi?[^] toC. 

16 lft(w)<0 

17 then for j = 1 to -t{u) : add rt[«] to G ; 

18 else for j = 1 to ^(m) : add T[m] to C . 

19 endif 

20 endfor 

21 Leti ^V{G) \img{f). 

22 Foreach u € I, do: 

23 Foreach w e Ng{u) , do: 

24 If w e /, then add AZ[«,ii'] to G. 

25 endfor 

26 If u 1^ I, then add |+)[^^] to G. 

27 Remove u from I. 

28 endfor 

29 return G. 



Algorithm 4. A procedure to construct a candidate circuit corresponding to a star 
decomposition of a geometry, realized over the gate-set {H, T, AZ}. 
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T operations on w to C, as well as controlled- Z operations between w and the 
other elements of N . Removing w from N after doing so will then prevent 
controlled- Z operations from being added twice for each pair of neighbors. 

• For each qubit u such that u ^ f{u), we may obtain the set of N neigh- 
bors of V = f{u) belonging to the output subsystem 0^* of the star geometry 
{Gu*, lu*, Ou^,)■ N will then consist of the neighbors of v excluding u and (if 
it is defined) f{v) — /(/(u)). We may then add the star circuit Cu* to C by 
adding operations between each w ^ N and v, followed by the components 
of the gate J{6) for the appropriate angle 9. For angles 6 constrained to inte- 
ger multiples of "74, this will consist of H[^] T*'"-'[m] for some integer array t 
given as input: if t{u) is negative, we may add products of T^[«] to C. 

After having added the star circuits for each star geometry, we then add the circuit 
corresponding to the geometry {G[I] , /, /) for / = V{G) \ img(/) . After constructing 
the set /, we iterate through the qubits u ^ I, adding to C the AZ operations acting 
on u and (if necessary) representing the fact that it ^ / by adding a preparation of a 
fresh qubit |+) in order to advance u as a fresh index in the circuit Cj^ Subsequently 
removing u from / prevents duplicate operations from being added to C. 



Run-time analysis of Algorithm 4. We may bound the run-time ofBuildCircuit 
as follows. Letn= \V{G)\,m~ |_E(G)|, and fc = |0|. We assume a list-like structure 
for the stable-index expression C and for the sets Ng{u), and an array-structure for /, 
where f{v) = nil denotes that v ^ dom(/). 

(/) For each qubit u = f{u), the number of gates to add corresponding to the ex- 
pression for ^z2(jcgg(M) above is 0(dcgQ(u)^). As the output subsystem of each 
star geometry in a star decomposition has the same size input and output subsys- 
tem, we have degQ(u) ^ \0\ = k for each such u; then the total work number 
of gates arising from these qubits is 0{k^n). 

(ii) For each qubit u ^ f{u), the number of gates to add corresponding to the circuit 
Cu* is 0(degQ(/(u)); summing over all such u, we have a total of 0{m) gates 
arising from these qubits. 

(Hi) For each qubit u E I = V{G) \ img(/), we have degQ{u) gates to add, plus 
possibly an initial state |+)[-] . Summed over all it e /, we again have at most 
0{m) operations arising in this case. 



As each operation requires time 0(1) to add, the run-time of Algorithm 4 is then 

0{k'^n + m). 

"While the semantic map we will consider does not involve the introduction of fresh qubits, such fresh 
qubits do not complicate the process of constructing candidate circuits, and so we include such preparations 
in |Algorithm 4| 
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3.5 Complete algorithm for the semantic map when \I\ = \0\ 



We may now describe a single, complete procedure Semantic (presented in Algo 



|rithm 5| l to perform a semantic map S with respect to the representations 7?.dkp and 



TZrbb represented by the constructions of Section 2.3 We may describe the operations 



Input : *P, a IWQC procedure in normal form. 
Output: Either nil, if *p is not (congruent to) the normal form of a 

procedure as in [Lemma 3.3 or a stable-index representation of a 



circuit C generated over the gate-set {H, T, /\Z, Zzzd\ci>2 which 
performs the same transformation as *p. 

; procedure Semantic(*P) : 

2 Let (G, /, O) be the geometry underlying 

3 Letn<- |y(G)|,TO^ |£;(G)|,fc^ |0|. 

4 Let M be a set of possible mediator qubits, initially empty. 

5 Let t : O*^ — > M be an array of multiples of ""A. 

6 Foreach o/7eraf/o« M^''^ *p, do: 

7 Set<(u)< 46l/7r. 

8 If 6* = Va , then add u into M. 

9 endfor 

10 If m> nk — {^X^) , then return nil. 

/; Let starDecomp ^ FindModStarDecomp(G, /, O, M). 

12 If StarDecomp = nil 

13 then return nil; 

14 else let (/, i) ^ starDecomp . 

15 If TestDependencies(^, /) = false, then return nil. 

16 LetG ^ BuildCircuit(G, i). 

17 return Remove I dops(G). 



Algorithm 5. A complete algorithm for a semantic map with respect to either the 
DKP or simplified RBB constructions, using the procedures defined in [Algorithm 1[ 
[Algorithm 4[ and [Algorithm 3| as well as a procedure Remove I dops which performs 
trivial simplifications of products of single-qubit gates. 

performed by Semant ic in three different phases, which we may describe as follows. 

Obtain the underlying structure of *p. We first obtain the geometry (G, /, O) un- 
derlying *p, as well as the set M of qubits measured with default angles ^A, in or- 
der to obtain a star decomposition. We also obtain the cardinalities n = \V{G)\, 
m = \E{G)\, and k = \0\. For each qubit u measured, we record the integers 
t{u) ^ ~A0/tt which corresponds to the power of T for which i7T*(") = J(-6l„); 
if the measurement at u corresponds to the execution of a J{6) gate in the candi- 
date circuit, we may use t{u) to determine the expansion of this gate in the gate-set 
{H,T,f\Z). 
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Rather than perform these operations Hterally as presented in Algorithm 5 we may 
obtain this data with a single pass of the operations of *p, constructing G, P, and 0° 
by inserting qubits and edges as each qubit is acted on by preparation maps, entangling 
operations, and measurements; we may obtain n, m, and k similarly. Using arrays of 
pointers to linked list nodes to represent P and 0°, we may then construct a similar 
representation of both / and O with a single traversal ofV{G). As we do so, we also 
assign values to the array t and insert qubits into the set M. As Cp has at most 0{n) 
preparation, measurement, and correction operations, and 0{m) entangling operations, 
it has size 0{m + n): this is then the time required for this phase of the procedure. 



Examine the geometry and the measurement dependencies of *p, 

in 



As we noted 



Section 3.3.4 a IWQC procedure cannot have a modified flow if it has more than 
nk — C"^^) edges, in which case it lies outside the range of both the DKP and RBB 
constructions by [Lemma 3.3[ we then return nil. We then attempt to obtain a star 
decomposition for (G, /, O) with the set M of potential mediating qubits: if this fails, 
we also return nil. 



By Theorem 3.1 we may characterize the dependencies of operations in procedures 
obtained from either the DKP or RBB constructions, relative to a candidate for the suc- 
cessor function for a corresponding circuit. For the star decomposition described by 
the ordered pair (/, L) returned in the previous step, there is a corresponding candidate 
circuit G for which / is an extended index successor function, as described in |Sec-| 
tion 3.4 In the case that |/| = \0\, there is no other candidate successor function for 
(G, /, O, M), by Theorem 3.4 Then the candidate circuit G described by (/, L) is the 



only one which describes a unitary transformation which may be performed by *p. 
The IWQC procedure <i>(G) has underlying geometry (G, /, O) by construction; if 
has the same operations (with the same dependencies) as the normalization of <&(G), 
they then perform the same transformation, and thus *p is a representation of G. Again 



by construction, the normalization of $(G) satisfies the conditions of Theorem 3.1 with 
respect to /. Then if TestDependencies(^P, /) returns true, we may provide the 
semantics for *p in terms of unitary circuits by constructing G; otherwise, we have no 
such guarantee, and so we return nil. 

We may compare m to nk — (''^^) in time polylog(n); and the respective run- 
times of FindModStarDecomp and T estPepende nc ies are 0{kn ) C 0{n^) 



and 0{kn ) C O(n^), by the analyses of Section 3.3.3 and Section 3.1.2 respectively. 



Constructing a circuit for *p. Provided that the procedure does not return n i 1 in the 
previous phase, the procedure *P performs a unitary transformation, and in particular is 
a representation of the "candidate" circuit G obtained from the star decomposition of 
(G, /, O). We therefore construct the circuit G, using the star decomposition described 
by (/, L). We produce a stable-index repre sentation for G, representing it in the gate- 
set {if, T, AZ} rather than the set B used in Lemma 3.3 this takes time 0{k^n + m) = 
0{k^n), as m ^ fcn — C^^^) by hypothesis. 

In the case of procedures *p obtained by applying the RBB construction to some 
original circuit G, we must perform some final simplifications. We call a final sub- 
routine Remove I dops(G) in order to remove superfluous powers of T which may 
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arise from the conversion of operations Zz = (T^ ® T'^)^Z, and which also removes 
products of the form = J(0)J(0) and HT^HT^HT^ = J C^M) J {%) J C^M) on 
individual qubits such as those introduced in the stable-index construction presented in 



Section A. 2. 2 As the dependencies of *p satisfies the constaints of Theorem 3. 1 rela 



tive to the successor function /, this circuit provides the semantics for *p as a unitary 



circuit, by the analysis of page 61 



The list L produced by FindModStarDecomp contains roots of star-geometries 
in a non-increasing order (in terms of the natural pre-order =^ for the successor function 
/ it produces), and the operations of C are in a similar order from the construction of 
BuildCircuit(G', /, L, t). We may therefore perform this simplification of C in 
a single pass from right to left, using buffers for each logical qubit of the circuit to 
maintain accumulated sequences of single-qubit operations performed between pairs 
of two-qubit operations. We may then cancel products of single-qubit unitaries that 
occur at the end of the buffer, as appropriate. The amount of time required to do so is 
then a constant factor time the size of the circuit, which by the analysis of page 62 is 
also 0(fc^n + m) = 0{k'^n). 

Remarks on circuit con gruence and total run-time. In order to consider the cor- 
rectness of Algorithm 5 we may consider without loss of generality circuits C rep- 
resented in the gate-sets Sdkp or Srbb which do not contain redundant operations 
J(0)J(0) or JCVa) JC"^) JC/z), and which do not introduce fresh qubits. We may 
proceed as follows. 

If we apply the DKP construction to such a circuit C over the gate-set ,8dkp, we 
obtain a procedure *P for which the qubits of 0° are labelled by deprecated indices in 
C, and P by advanced indices; where the entanglement graph is the same as the inter- 
action graph of C\ and where the default measurement angles are the additive inverses 
of the angles of the J{0) gates. Furthermore, the geometry (G, /, O) has a modified 
flow (/, =^), in which / is the index successor function of C. Correspondingly, the can- 
didate circuit C for *p has an interaction graph given by the entanglement graph of *p, 
and a successor function which agrees with /. The single-qubit gates corresponding to 
the edges uf{u) in the interaction graph correspond to operations J{—Oy), where 0^ is 
the default measurement angle of u; and all other edges correspond to /\Z operations. 
We may then define a bijective mapping between the terms of C and the terms of C, 
via the structures induced in and recreated from *p. 

For the RBB construction, the analysis is similar: the only changes are the fact 
that the edges of the interaction graph of C which correspond to operations Zz must 
be subdivided by a single vertex to yield the entanglement graph of *p, and paths of 
the form J(0) J(0) or J(V2) J(^/2) J^Va) may be inserted to achieve the topological 
constraint of being an induced subgraph of a two-dimensional grid. In constructing 
the candidate circuit C, the subdivided edges corresponding to Zz operations are re- 
merged, and the superfluous products of single-qubit unitaries are once more removed. 
In this case, we may then also define a bijective mapping between the terms of C and 
the terms of C, via the structures induced in and recreated from *p. 

Thus, Semantic defines a semantic map for both the DKP and (simplified) RBB 
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constructions, restricted to unitary bijections. The run-time of Semantic is domi- 
nated by the time required to compare the dependencies of *p against the criteria of 
Theorem 3.1 which is 0{nm) C 0{kn?). 



4 Extensions and limitations 



The analysis of IWQC procedures presented in Section 3 in terms of modified flows 



and star decompositions, may be extended in a natural way to more general formula- 
tions of the IWQC model. At the same time, this approach has limitations in its ability 
to describe even some of the earliest known constructions for IWQC computation. In 
this section, we will describe in brief some of these extensions and limitations. 



4.1 Extending beyond measurement in the XY plane 

The model of measurement-based computation which we have focused on in this article 
(following the treatments in 1 1 1 and 0|) is one in which all single-qubit measurements 
are performed in the XY plane, i.e. with respect to uniform superpositions of the eigen- 
states of the Pauli Z operator. A natural extension of this is to allow measurements in 
any basis in a "Pauli plane", i.e. for any uniform superposition of the eigenstates of X 
and Y as well. We may define the notation |±xyg) = \±g) (denoting the fact that these 
states lie in the XY plane), and then define 



l±yze) = 
l±xz«) = 



1) ± 



^(|0> 



-ID 



\o)-^\^) 



(4.1a) 
(4.1b) 



which denote bases lying in the YZ and XZ planes respectively. Then, we define the 
measurement operators 



l(YZ,( 



l(XZ,( 



'(p) = (+yze|p |+yze),®|0)(0|,[^, 

+ {-yZg\p\-)/Ze)a®\l){l\, 



(+xze|p 



hxze>a® |0)(0|,[^, 

+ {~XZe\p |-XZ(,)a® 



(4.2a) 



(4.2b) 



(4.2c) 



We may then consider extensions of the analysis of Section 3 in measurement-based 
models which include operations of the form above. 



4.1.1 Stabilizers, certificates of unitarity, and local graph structure 



Recafl from Lemma 3.3 that we may construct procedures performing unitary trans- 
formations in the IWQC model by performing the following operation after each mea- 
surement IVli^^'^': 
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' If "Ya, we perform an X operation on some specially designated qubit v ^ u 
(which we denote f{u) for each u), and a Z operation on each of the neighbors 
of V except for u itself, if the state of u is projected onto |-xyg). 

• If 6' = "/a, either we perform the operation described above, or we perform a Z 
operation on each neighbor of u, if the state of u is projected onto l-xy^). 

In either case, the operation which is performed after each measurement may be de- 
scribed by an operator 

Bu = I II zA , (4.3) 

up to an operation acting on u itself (which we describe as being destroyed by the 
measurement and which thus cannot be carried out). This unitary operation is the PauU 
group on V{G), and in particular stabilizes the state \'^g,i) defined by 

I^-C/X^-gjI = EgN,c(|^)(V'|) (4.4) 

for any ji/i) e Tif^ . In particular, as the preparation maps N^, adjoin qubits prepared in 
the + 1 eigenstate of the X operator, and the Ey^^ maps perform two-qubit AZ operations 
(which lie in the Clifford group) on pairs of qubits, the state of the system just prior to 
any measurements is a stabilizer code |5 1 §g.i generated by the operators 

= X^Yl Z^ = Eg(X,) . (4.5) 

The modified flow conditions may then be considered as a certificate that there 
exists an order in which the qubits of may be measured so that one may use the 
stabilizer formalism |5| to simulate the post-selection of the state |+xyg) for each mea- 
surement (by applying suitable corrections should the measurement yield j-xy^) in- 
stead). 

We may observe that the analysis of measurements by a single-qubit observable 
Ov in the stabilizer formalism requires only that Oy either commutes or anticommutes 
with every generator of the subgroup of the PauU group which stabilizes the state of 
the system. We may observe that XY-plane measurements such as those above can be 
described by measurements with respect to observables 

O^Y'^ = l+xyeX+xy^l - |-xy,)(-xy,| = cos(0)X + sin(0)y ; (4.6) 

such operators always anticommute with the Pauli Z operator, and also anticommute 
with X when = ^/2. Then, a modified flow function / : 0° — > P for a geometry 
{G, I, O, A/) describes a mapping from u (with an associated observable 0^^'^") to 
some qubit v — f{u) such that Ky anticommutes with 0^^'^" . The constraints on the 
partial order =<; ensure that, for any linear measurement order extending the gen- 
erators which stabilize the state just prior to the measurement of u consist of Kfi^y^, 
and operators which all commute with u. Then, the evolution of the state-space system 



69 



under these measurements may be efficiently described via the stabilizer formalism as 
a transformation of stabilizer codes (although how individual states within these codes 
transform may be difficult to simulate). As the modified flow conditions ensure that 
each measurement observable has a corresponding generator with which it anticom- 
mutes, the overall transformation will then be unitary. 



A straightforward extension of modified flows. The fact that qubits u e M may 
have either u ~ f{u) or u = f{u) stems from the fact that 0^^'^" — Y in this case: 
that is, the measurement basis for such qubits lies in the intersection of the XY and 
YZ planes. The condition u ^ f{u) in this case is common to all other measurement 
angles, as this condition suffices for treatment via the stabilizer formalism for observ- 
ables which anticommute with Z. We may then extend the notion of modified flows 
as follows to also include other measurement observables which, like Y, anticommute 
with X: 

Definition 4.1. For a tuple (G, J, O, Mxy, Myz) consisting of a geometry (G, /, O) 
and sets Mxy, Myz ^ O'^ such that 0'^^ = Mxy U Myz, an extended flovJ^is an 
ordered pair (/, consisting of a function / : 0° — > P, and a partial order =<; on 
V{G), such that the conditions 

f{v) ^ V , for w e Mxy \ Myz ; (4.7a) 

f{v) =v, for w e Myz \ Mxy ; (4.7b) 

f{v) r^v or f{v) = v, for w e Myz n Mxy ; (4.7c) 

V 4 fiv) ; and (4.7d) 

w ^ f{v) =4> V ^ w (4.7e) 

hold for all v eO^andw e V{G). 

In the above definition, the sets AIxy and Myz correspond to qubits whose default 
measurement basis lies in the XY-plane or the YZ -plane, respectively; the intersection 
are those qubits whose measurement basis is the eigenbasis of the Y operator, which 
is common to both. Modified flows then correspond to the special case where Myz C 
Mxy. 

As with modified flows, an extended flow (/, =4) for a tuple (G, /, O, Mxy, Myz) 
certifies that any procedure on the geometry {G,I,0) which performs the appropri- 
ate types of measurement will perform a unitary transformation, provided also that it 
performs the appropriate corrections (or measurement adaptations) in response to the 



measurements. To this effect, we may prove an extension of Lemma 3.3 as follows: 



Lemma 4.1. Suppose (/, =<;) is a modified flow for (G, I, O, Mxy, Myz). Let v ^ w 
denote the adjacency relation ofG, and let A : V{G) — > {XY, YZ} indicate the plane 
of measurement for each qubit, with X{u) — XY or A(u) = ^T-for u G Mxy H Myz, 



"The definition of extended flows here is a refinement of the definition presented in (9]: the definition 
there can be recovered by requiring Mxy and Afyz to be disjoint, in which case the set T in Definition 3-17 
of (21 coiTesponds to -A/yz ■ 
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Figure 12: Illustration of the geometry for the procedure 333d defined in ( |4.10| l, to- 
gether with the "default state" l+yZg) to be measured for the qubit a in that procedure. 
(In the case that the measurement of a yields |-yZg) instead, each element of / = O is 
to be acted on with a Z operation.) 



and \{u) = P for u G Mp otherwise. Then for any linear order ^ extending ^, the 
measurement procedure 
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(4.8) 



performs a unitary transformation. 

Again, as we described above for modified flows, this may be proven for extended 
flows by simulating the evolution of the state-space via the stabilizer formalism, steer- 
ing the state after each measurement to the one which would arise upon selection of 
the l+xyg) or l+yzg) result as appropriate. 



Generalization of Z2Zd operations using extended flows. An alternative proof would 
be to extend the analysis of the proof of Lemma 3.3 by considering the operation 



d times 



(4.9) 



generalizing the operation Zzzd defined in <|4.9[). We may define a IWQC procedure 



n 



l(YZ,( 



d 



Na : 



(4.10) 



the geometry for this procedure is illustrated in Figure 12 Following the analysis of 
Section 3.8 of fT9l, we then have 



(4.11) 
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The procedure of ( |4.10| l is a straightforward generalization of the procedure 333'' = 
333'''"'^^ • In the same way that we map qubits u such that u = f(u) to a circuit 
Cu = ZzZd acting on the neighbors of u for a modified flow as a part of a circuit 
construction, we may map qubits such that u = f{u) to a circuit C'„ — Zzz^^, where 
|±yZe^ ) is the default basis for the measurement of u. 



Candidate circuits and star decompositions for extended flows. Just as geometries 
with modified flows correspond to circuits over the gate-set {J{0), ^Z, Zzzd}g^iL-z 11^2 
via Lemma 3.3 we may define a map from circuits over {J{0), ^Z, ^zz^lee^z d>2 
to the IWQC niodel. For an extended flow, we may apply the same definition for star 
geometries as for modified flows: the only extension required in the correspondence 
between star decompositions and extended flows is that qubits u g Myz \ AfxY are nec- 
essarily both the root and the center of a star geometry with respect to an extended flow. 
It is then possible to show that star decompositions with respect to extended flows may 
be obtained with only minor modifications to the procedure FindModStarDecomp 



defined in Algorithm 3 and its subroutine RemoveStarGeom defined in Algorithm 2 



In the former case, the set M on line 7 plays the role of Myz C Afxy of qubits which 
may be both the root and center of a star geometry: it suffices them to replace M with 
Myz. In FindModStarDecomp, we must refine the condition on [line 23j as there 
may not be qubits w e il/yz \ ^/xY which cannot be the center of a star geometry 
without also being the root: we must then revise the condition to require that 



D has only one outbound arc w v and v G M) 



XY, 



(4.12) 



the second condition of which was guaranteed for the special case of modified flows. 



Applying the same analysis as for Algorithm 3 this modified algorithm then finds a 
star decomposition for (G, /, O), with respect to the constraints described by Mxy and 
Myz on the possible relationship between the roots and centers of each star geometry. 
In particular, as the proof of the extremal result of [Lemma 3. 8| carries forward for 
extended flows, the run-time of the modified algorithm is also 0{kn). 



4.1.2 Generalized certificates of unitarity based on stabiUzers 

In the preceding section we observed one sense in which modified flows may be gen- 
eralized, by observing that they essentially provide sufficient information to certify via 
the stabilizer formalism that a IWQC procedure performs a unitary transformation. 
This is done by attributing some generator to each qubit m e to be measured, 
where the measurement on u may be described by some observable 0„ which anticom- 
mutes with A'„, but commutes with for all w ^ u. 

Similar, but more general, criteria that extended flows were defined in [4 1 in gener- 
alized flows and Pauli flows, which extend further to arbitrary generators of the initial 
stabilizer group {Ku,)^^jc (and thus to non-trivial products of the operators K^) when 
XY, YZ, and XZ-plane measurements are all allowed. We may paraphrase the defini- 
tion of Pauli flows as follows: 

Definition 4.2. Let (G, /, O) be a geometry, and let Afxy, Myz, Mxz C 0° such that 
O'^ = Mxy U M\z U AIxz- Let V denote the power-set function, and oddG(S') the 
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set of vertices in G which are adjacent to an odd number of elements of a set 5 C 
V{G). Then a Pauli flow for (G, /, O, Mxy, Myz, A/xz) is a tuple (g, consisting 
of a function g : — > ^(^°) and a partial order ^ which satisfy the following 
conditions^] 



w e g{v) U odd{g{v)) =^ v ^ w, 
V e odd{g{v)) \ (7(w) f G A/xY 
f G g{v) \ odd((7('i;)) v G Afyz 

u G g{v) n odd((jf(w)) =4> w G Afxz 



(4.13a) 
(4.13b) 
(4.13c) 
(4.13d) 



From this definition, we recover extended flows if we require that the sets g{u) be 
singleton sets for all u. (The definition of generalized flows may be recovered by 
requiring that A/xy> Myz, and A/xz be disjoint.) 



Based on these conditions, we may obtain a further extension of Lemma 3.3 which 
we paraphrase from Theorem 4 of (4'\ as follows: 

Lemma 4.2. Suppose {g, =4) is a Pauli flow for (G, J, O, Mxy, Afyz, Afxz)- Let v ^ w 

denote the adjacency relation of G, and A : V{G) > {XY, YZ, XZ} indicate the 

plane of measurement for each qubit so that for each qubit, X{u) = P only if u G Mp. 
Then for any family {f^ti}„gQc of measurement angles satisfying 



u G XZ n YZ, or 
M G XY n XZ and \{u) = XY 



=^ e^e {0,^} , 

0„G {-^/2,V2} 

and for any linear order ^ extending =<!, the measurement procedure 



u G XY n YZ, or 
u G XY n XZ and \{u) = XZ 



(4.14a) 
(4.14b) 
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(4.15) 



performs a unitary transformation. 

This represents a further extension of the possible range of IWQC procedures 
which may be identified as performing unitary transformations; the proof of the above 



^We may recover the definition of Pauli flows found in 14] by defining a function A from O"^ to 
{X, Y, Z, XY, YZ, XZ}, sucli that 

• X{v) = Xifv& Mxy n A/xz, 

• X(v) =Yifv e Mxy n J\/yz, 

• \{v) =Zifv e Mxz n Mxz, and 

• \(v) = P if V £ Mp otherwise; 

the fact that O'^ = Mxy U Afyz U A/xz then implies that the conditions of |Definition 4.2| and Definition 5 
of f4\ are equivalent. 
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Figure 13: Illustration of the geometry for a qubit-reversal procedure of Section IV.B 
of 121, and a partial decomposition into star geometries. Each qubit in V{G) \ O is to 
be measured in the XY plane with default angle 0: when this yields the measurement 
result s[(t] = for each u G V{G) \ O, the result is to perform a reversal of the vertical 
ordering of the logical qubits: input states a,b.c.d of the input subsystem (with a-d 
ordered from top to bottom) are mapped to jV') w,v,u,t of the output subsystem (with t- 
w ordered from top to bottom). Notice that the residual geometry in the decomposition 
on the right has no maximal star geometries, as every element of O' is adjacent to more 
than one element of V{G) \ O', and every element of V{G) \ O' is adjacent to at least 
one other element of V{G) \ O'. Thus, this geometry has no star decomposition, and 
therefore no modified flow. 



Lemma may again be obtained by considering how each measurement may be certified 
to transform the state-space unitarily conditioned on a preferred measurement result, 
where the transformation of the state-space and the corrections to apply are described 
by the stabilizer formalism. 

In order to use this to provide semantics for such IWQC procedures in terms 
of unitary circuits, we may consider how the concept of star geometries generalize 
from extended flows to Pauli flows. The natural approach, rather than attempt to 
retain the star-graph structure, would be to consider geometries with a single root 
u, multiple "centers" v e g{u), and whose output subsystem consists of the qubits 
[g{u) U odd{g{u))] \ {u}. However, it is not as clear what the circuit corresponding 
to such a geometry would be, nor what results may be shown for the decomposition of 
a geometry (G, /, O) into such geometries, as these lack the simple local structure of 
star geometries. How a Pauli flow structure may be exploited to produce more general 
semantic maps for IWQC procedures in unitary circuit models is an open question. 



4.2 A simple IWQC construction lacking a star decomposition 



Lemma 4~2] above presents a IWQC procedure which performs a unitary embedding. 



but for which the techniques of Section 3 do not apply. Another example is the pro- 
cedure for performing reversals of a linear array of logical qubits presented in Sec- 
tion IV.A of ||2|. Figure 13 illustrates the geometry underlying this IWQC procedure: 
the entire procedure may be obtained from the specification that each qubit in 0° is 
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to be measured in the XY plane with the default angle 0, and noting that as this rep- 
resents selection of the +1 eigenstate of the observables X„ for u E O*^, from which 
suitable correction operations may be obtained by the stabilizer formalism. However, 
as Figure 13 also illustrates, this procedure has no star decomposition, and so cannot 
be mapped to its meaning as a reversal of a linear array of qubits using the techniques 
of lSection 31 

Despite this, we may gain some information about the qubit-reversal pattern using 
modified flows by an appropriate extension of the input and output systems, as exhib- 
ited in ifTTl . It is easy to see that if we include the left-most element of each row into 
/, and the right-most into O, the resulting altered geometry has a flow with a successor 
function which maps each vertex to the one immediately to the right. Chains in the 
natural pre-order =^ have monotonically increasing column numbers: and therefore =<; 
is antisymmetric. Given that the default measurement angle for each qubit is zero, we 
may obtain a candidate circuit for this procedure, as illustrated in [Figure 14| The final 
circuit illustrated there can easily be described as a reversible classical computation: 
for seven input bits {vq, . . . , wg), the computation performed by that circuit is 

{vo, Vi, V2, V3, V4, W5, Ve) 

I > ("6, W4 © W5 ® V6, Vi, V2(BV3(B V4, V2, Vq ® Vi (B V2, Vq) . (4.16) 

If we apply this circuit instead to qubits, with the odd-indexed qubits initialized to the 
1+) state (i.e. a uniform superposition of the bit-values described above), those qubits 
remain in the |+) state, and the even-indexed qubits are reversed, as required. 

We may then understand the qubit-reversal pattern in |2| as corresponding to a 
circuit which prepares auxiliary qubits in the |+) state, and subsequently removes them 
by an X observable measurement. The fact that we are able to obtain the circuit in this 
case is possible because the highly structured geometry makes it easier to guess which 
qubit measurements correspond to terminal measurements in a unitary circuit. Without 
exploiting this structure, however, there are no known approaches to algorithmically 
provide semantics for such a measurement pattern in terms of unitary circuits. 



5 Review and open problems 

In this article, we have presented an algorithm for obtaining the "meaning" of IWQC 
procedures *p in terms of unitary circuit models, using a combinatorial analysis of the 
structures underlying procedure ?p arising from standard constructions of IWQC pro- 
cedures. In doing so, we have identified a correspondence between elementary combi- 
natorial structures {star geometries) and elementary circuits {star circuits) which, given 
a IWQC procedure whose geometry which may be decomposed into these elementary 
structures, allows the construction of a corresponding circuit. Together with a charac- 
terization of the dependencies of operations on classical measurement outcomes, this 
enables us to efficiently construct a circuit C for a IWQC procedure *p, and to de- 
termine whether C and *p perform the same operation, when *p has input and output 
subsystems of the same size. 

We have used the gate model {H, T, ^Z} as a reference model of circuits, using 
the related gate-sets Bdkp = {JiS),^Z}g^^^ and Srbb = {J{S), Zz}g^^j^ to de- 
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Figure 14: Illustration of the qubit-reversal procedure of Section IV.B of f2\, with 
augmented input/output subsystems, and corresponding candidate circuits (represented 
in the middle in the gate-set {H, T, ^Z}, and in the bottom with CNOT gates). 
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scribe constructions of IWQC procedures. However, essentially all of the results of 
this article generalize to {H, R{9), ^Z}g^^, where 



m 



e^«/2 



(5.1) 



we then have J{e) = HR{e), and the gate sets {H, R {0),/\Z}f,^^ , { J(6I), AZj^p, , 
and {J{9), Zz}g^g^ are all parsimonious in the sense of Definition 2.1 

There are three natural open problems: how this analysis may be extended to the 
case of unitary embeddings in general, how this may be extended to admit an analysis 
for quantum circuits with ancillas (qubits which are prepared, operated on, and subse- 
quently removed without creating mixtures of states in the other qubits), and how we 
may extend this analysis to include stabilizer certificates such as Pauli flows. 



General unitary embeddings. The problem presented by the case of a unitary em- 
bedding with |/| < |0| is that a star decomposition of the geometry may not be unique. 
Consequently, there may be many modified flow functions against which the dependen- 
cies of a IWQC procedure must in principle be tested. We may relieve the situation 
if we can use the dependencies themselves to narrow the search space for a single 
modified flow; this seems a likely approach to extending the analysis of this article 
to a semantic map for the DKP and (simplified) RBB constructions, without further 
constraints. 



Quantum circuits with ancillas. The qubit-reversal procedure described in Sec- 



tion 4.2 presents a situation where the introduction and removal of ancillas can be 
used to obtain a unitary circuit which corresponds to a IWQC procedure. Combina- 
torial tools for the "detection" of instances in which a qubit may represent the state of 
an ancilla qubit would yield more general tools for presenting the semantics of IWQC 
procedures in terms of unitary circuits, possibly permitting a semantic map to be de- 
fined for a construction employing all of the sub-procedures described in Section IV of 
|l2|, and for IWQC procedures in general arising from circuits involving partial trace 
operations on qubits which are disentangled from the rest of the system. 



Extension to more general stabilizer certificates. As noted in [Section 4.1.1| and 
Section 4. 1 .2 modified flows essentially represent a means via the stabilizer formal- 



ism to certify that, for a given geometry (G, /, O), there is an efficiently describable 
IWQC procedure having that geometry (G, /, O) and which performs a unitary trans- 
formation. We can then consider how to compare an arbitrary procedure with geometry 
(G, /, O) to this canonical geometry. In Section 4. 1 . 1 we showed how such certificates 
may be extended to more general gate-sets: how we might exploit the definition of PauU 
flows in Definition 4.2 to obtain more general semantic maps is an open problem. 
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A Stable-index tensor notation 



A.l Sum-over-paths interpretation of stable-index tensors 

The informal semantics of stable-index tensor notation is to provide a simplified de- 
scription of circuits in terms of computational paths, as follows. 

Stable-index notation for circuits is partially inspired by the path-labelling of cir- 
cuits by Dawson et al. |21 1, which is explicitly concerned with descriptions of circuits 
in terms of computational paths: the tensor indices of the stable-index expressions cor- 
respond to labels of "wire segments" of [21], which are separated by Hadamard gates 
(or more generally, other operations which do not preserve the standard basis). For a 
given operator, a stable index s then represents a qubit for which the standard basis 
states |0) and |1) remain unchanged (i.e. the density operators \s) {s\ are unaffected) by 
the operation for s e {0, 1}. Over a single "wire-segment", the reduced state of such a 
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qubit will then remain unchanged, even if one measures this state in the standard basis 
at the beginning of the wire-segment. 

As with tensors in Einstein notation, we may describe a "path" as a mapping of 
each index (whether bound or free) to some fixed value; the amplitude associated with 
a path is the product of the corresponding coefficient of each factor involved in the 
expression. The convention of summing over bound indices (whose possible values 
may be thought to represent intermediate configurations which were not observed, and 
therefore not distinguished from one another) then represents a summation over in- 
distinguishable paths. The simplification introduced by stable-index expressions cor- 
responds to explicit recognition when operations leave some of the the "dynamical 
variables" associated to a path (represented by tensor indices) unaffected. 

A.2 Constructing stable-index representations 

We may automatically generate a stable-index representation S for a unitary circuit 
from a representation as a composition of operators C on specified qubits in an efficient 
manner. We consider two distinct methods: one which is appropriate for a network 
without any topological constraints on the operators, and one which is motivated by 
the constraints described in [Section 2.3.31 

A.2.1 Construction for topologically unconstrained circuits 

Let L be the set of qubit labels in C. We may then transliterate the terms in C, defining 
a set of tensor indices incrementally as we do so, as follows: 

1 . For each qubit u e C, we define the tensor index vq ■ We designate vq as the cur- 
rent tensor index for v. (Throughout the transliteration of the circuit, the index 
designated as "current" for a qubit will be updated.) We also start with an ex- 
pression 5' consisting of the null sequence of operators: this may be regarded as 
representing a stable-index representation for the the empty product (the identity 
operator). 

2. We translate each operator Ua.b,... in C taken in sequence, as follows. For each 
qubit V G {a, b, . . .} on which U acts (again, in sequence), let vj be the current 
index for v, and perform the following: 

(0 If Ua,b.... introduces w as a "fresh" qubit (or a new qubit which is possibly 
entangled with the others), then t^o is still the current index for v. We then 
use Wo as the advanced index corresponding to v, without a corresponding 
deprecated index. 

(//) If Ua,b.... not preserve each of the standard basis states on v (that is, if 

Ua,b,... C^Ib,... 7^ |s)('S|t, for some s € {0,1}), then we define a 

new tensor index Uj+i for v. We then use Vj as the deprecated index corre- 
sponding to V, and Vj+i as the corresponding advanced index. We update 
the current index for v to Wj+i. 
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(Hi) If Ua.b,... preserves each of the standard basis states on v (that is, for each 
s e {0, 1} we have Ua,b,... \s){s\,^ ^ = |s)(s|,y), then we use vj as a 
stable index corresponding to v, without defining any new indices]^ 

(iv) If Ua,b,... removes a qubit v (e.g. describing an operation conditioned on a 
a particular measurement outcome), then we use Vj as a deprecated index 
corresponding to i;, without a corresponding advanced index. 

This produces sequences of stable indices s and sequences of matched advanced 
and deprecated indices a and d (possibly with placeholders); we then translate 
Ua,b.... as U[s , and append it to S. 

3. Having translated each of the terms in C, the sequence S' is a stable-index tensor 
expression corresponding to C. 

Note that in particular, every new index which is advanced by a gate U[s must be 
distinct from any other index being advanced, as well as from all of the indices which 
have occurred to that point. This ensures that the resulting product is acyclic if the 
original network C was acyclic. 

Note that the index successor function / ( [Definition 2.3 1 can be easily described 



from this construction: / has as a domain all those indices vj such that vj^i is well- 
defined, and Vj+i = f{vj) for all such indices. 

A.2.2 Construction suitable for use in the simplified RBB construction 



Motivated by the constraints described in Section 2.3.3|on IWQC procedures whose 



entanglement graph is an induced subgraph of the grid, we also present a construction 
for stable-index expressions for circuits with a linear-nearest neighbor topology with 
the gate set Srbb = {•^(^): 'S^legmz' which the two-qubit gates 

(0 act "stably" on both their indices, 

(//) can only act on indices of equal "depth" from the input — that is, can only act 
on pairs of indices Vj and Wj, where Vj — /■' (wq) ™d Wj — /^ (u'o) for some 
indices vo,wo ^ dom(/), where / is the successor function of the stable-index 
expression; and 

(///) can only act on a pair of indices Vj , Wj if there is no two-qubit gate acting on 
f j_i and Wj^i. 

Given a circuit C using the gates of Srbb, we may produce such a stable-index rep- 
resentation (and at the same time explicitly construct its interaction graph, applying 



Definition 2.4 in the special case of interactions of at most two qubits) as follows. 



adapting the algorithm of Section A. 2.1 above 



^In the definition of stable-index notation in |Section 2.1.3| all stable indices must occur as a contiguous 
block at the beginning of the linear ordering of the qubits operated on. Generalizing this definition, we may 
easily alternate between sequences of stable indices and sequences of advanced and deprecated indices. If 
we choose not to generalize stable-index notation in this way, we may only represent a qubit v with a stable 
index if all preceding qubits on which U operates are also represented by stable indices; otherwise, we must 
apply case (n') above regardless of how U acts on v. 
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1. For each qubit v acted on by C, we set the current index to vq. We also initialize 
the graph G to the totally disconnected graph whose vertices are the indices vq . 
We start with an expression S consisting of the null sequence of operators: this 
may be regarded as representing a stable-index representation for the the empty 
product (the identity operator). 

2. We translate each operator J{9)a or Zzab in C taken in sequence, as follows. 
For each qubit v on which the gate acts, let Vj be the current index for v. For a 
gate J{6), we append to S, add the index Vj^i and the edge VjVj+i 
to G, and set the current index to Vj+i. Otherwise, for a gate Zzy^, we perform 
the following operations: 

(0 We ensure that the constraints on successive applications of two-qubit gates 
is satisfied as follows. If the current indices for v and w are Vj and Wj for 
a common value of j, and there is an edge of the form vj-iWj-i in G, we 
append the expression 



J(0) ^ J(0) J(0) J(0) 



(A.1) 



to S, and add the edges VjVj^i , Wj-|_iWj+2 , ''^jWj+i , and Wj^iWj^i (along 
with the indicated indices) to G. We then set the current indices of v and w 
to and respectively. 

(ii) We then recursively perform the following, letting Vj and Wk stand for the 
current indices of v and w at each iteration: 

• If |fc — j| ^ 3 or |fc — j| = 2, we append either 



J(0) 



"3 + 1 



J(0) 



"3 + 1 



or 



J(0) 



Wk-\-2 



Wk + 1 



J(0) 



Wfc + 1 



(A.2) 



to S, according to whether k > j or j > k respectively, adding the 
appropriate indices and edges to G. We set the current index either of 
u or w to or Wk+2 accordingly, and repeat for the new current 
indices. 

If 1^ ~ il =3 or \k — j| = 1, we append either 



V I + 2 ^ I "j + 1 ^ ' I Vj 



Wk + 2 



Wk + 2 
Wk + 1 



Wk + 1 
Wk 



(A3) 



to S, according to whether k > j or j > k respectively, adding the 
appropriate indices and edges to G. We set the current index either of 
u or w to Vj+3 or Wk+a accordingly, and repeat for the new current 
indices. 

• If j = k, we append Zz[vj,Wj] to S, add the edge vjWj to G, and 
stop. 

3. Having translated each of the terms in C, the sequence 5 is a stable-index tensor 
expression corresponding to C, and G its interaction graph. 
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In the above construction, note that the final index vi for each qubit v may not have 
the same value of i; and that such stable-index expressions may both begin and end 
with operations Zz[vj,wj'^ , These degrees of freedom may complicate the procedure 
of composing two such circuits. These issues may both be mitigated by appending ap- 
propriate products of J (6) operators which evaluate to the identity, of differing lengths 
for each qubit; this can be done to both uniformize the value £ of the final indices ve, 
and also ensure that any Zz operation is a distance of at lest 2 from the end of the 
circuit. However, for the purposes of this article, such operations will not be necessary, 
as we do not consider the composition of two such expressions. 



A.3 Congruence of circuits and isomorphic stable-index expres- 
sions 



In this section, we prove Theorem 2.1 (page 11), and expand on the role of parsimo- 



nious sets of gates in this result. 

Lemma A.l. Let Ci and C2 be two sequences of unitary operators on a common 
set of qubits, composed from a common set of unitaries. If the stable[index tensor 
representations of Ci and C2 are isomorphic, then Ci and C2 are congruent up to 
re-ordering of commuting operators. 

Proof. Consider the coarsest equivalence relation = on circuits which consist of per- 
mutations of the gates of Ci, such that C = C if C differs from C by a transposition 
of commuting gates. By definition, the equivalence classes of = are sets of circuits 
which are congruent up to re-ordering of commuting operators. Consider stable-index 
tensor expressions S and S' arising from two unitary circuits C and C, where S acts 
on the same index set as S', and S differs from S' only by a transposition of two terms 
Ui[si ^] and [/2[s2 g|] , where without loss of generality the former precedes the lat- 
ter in C, and vice-versa in C . Then, the corresponding gates Ui and U2 in C occur 
with Ui preceding U2, and occur in C" in the opposite order. Consider the index sets 
for these two operations in the stable-index tensor expression: 

• If the sequences of indices (si ; ai ; di ) and (s2 ; a2 ; do not overlap, the oper- 
ations Ui and U2 act on disjoint sets of qubits, and so they commute. 

• Suppose the sequences of indices (si; ai; di) and (S2; a2; d2) overlap, and let 
V be an index occurring in both sequences. In the construction of S, because v 
occurs in f/2[s2 g^] , it cannot be a deprecated index in Ui[si g^] ; and because 
it occurs in Ui[si g^] . it cannot be an advanced index in J72[s2 g^] . By the 
construction of S', we similarly have that v cannot be a deprecated index in 
C/2[s2 g^] or an advanced index in Ui[si g^] ■ Then, w is a stable index of both 
terms. 

Therefore, the only qubits v on which Ui and U2 both act are those whose standard 
basis states are preserved by both Ui and U2. Then, Ui and U2 commute, in which 
case C ^ C". 

Let Si and ^2 be the stable-index tensor representations of Ci and C2, where the 
isomorphism is given by an index relabelling map A and mapping of the terms t. Let 
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5*2 be the representation obtained by applying A to all of the indices in S2 '■ then 5*1 
and ^2 consist of a product of the same terms in different orders. Without loss of 
generality, we will then suppose that ^2 and Si differ only by a reordering of terms. 
Because the terms of 5*1 differ from those of 5*2 by a permutation, by induction on the 
decomposition of this permutation into transpositions, we therefore have Ci = C2. ■ 

Whatever set of gates is used in a stable-index tensor expressions, rearranging the 
terms in the expression not only preserves the meaning (by summation over indices 
which are both advanced and deprecated) as a unitary operator, but also preserves 
the relevant structural information about the circuit — preserving the necessary or- 



der information in the way described following Definition 2.3 while discarding such 
information for each pair of commuting operators. As a consequence, we may regard 
permuting terms of such an expression as a way of exploring the elements of the con- 
gruence class of a circuit. 

However, in the case where two gates commute without preserving the standard 
basis states of their common operands, we may obtain congruent circuits which do not 
have isomorphic tensor index expansions: 

Example 5. Let U G U(2) be a non-diagonal operator, and define the operation NJ 
as in ( |2.12a| l. Consider the two circuits MJc^b ^Uafi Ha and Wa,b ^Uc.b Ha on 
three qubits a, h, c (where the subscripts here denote the linear order of the qubits as 
operands of NJ and H): these have the respective stable-index expressions 





m 




H 




f)2 


bi 




aa 



and 



NJ 





NJ 




H 




b2 


bi 




aa 



(A.4) 



These two expressions are non-isomorphic: in the first expression, the leftmost factor 
contains the indices c ( which is neither advanced nor deprecated in the entire expres- 
sion) and 63 (which is not deprecated in the entire expression); the second expression 
does not contain any factor with these properties. 

In the example above, the inequivalence of the two stable-index expressions is due 
to the fact that the two operations ^Ua,b and ^Ub.c preserve the eigenbasis of U on b, 
rather than the basis |0) , The restrictions imposed on parsimonious sets of unitary 
operations prevent precisely this problem from arising, which allows us to prove the 
following: 

Lemma A.2. Let Ci and C2 be two unitary circuits composed from a common par- 
simonious set of unitaries. If Ci is congruent to C2 up to re-ordering of commuting 
operators, then the stable-index representations ofCi and C'2 will be isomorphic. 

Proof. Let Kq ^ Ki ^ ■ ■ ■ = Kg he a sequence of circuits which differ only by a 
transposition of two commuting gates, where Ci = Kq and C2 = Kg. We may show 
by induction on i that the stable-index tensor expressions for Ci and C2 are isomorphic; 
it suffices in this case to prove the case £ — 1. 

Let 5*0 and 5*1 be the stable-index representations of A'o and Ki . Up to a relabeling, 
we may assume that the indices of Sq and 5*1 are all of the form vo,vi, . . . for different 
qubits V acted on by Kq and Ki : we set vq to be the first index corresponding to a 
given qubit in either Sq or ^i, and Wj+i to be an advanced index corresponding to each 
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deprecated index Vj for j G N. Let Ui and U2 be the pair of gates whose order differs 
between the two circuits; and consider the truncation of Sq just prior to the terms 
corresponding io Ui and U2, and similarly for S'l. Because Kq and A'l differ only 
by a transposition of commuting gates Ui and U2, S'q = S[. Let Sq and S"" be the 
truncations of Sq and Si just after the terms corresponding to Ui and U2 ■ 

• If Ui and U2 act on disjoint sets of qubits, the terms corresponding to C/i and 
U2 in S'q may differ from those in S" only by the indices which are advanced 
in each; and as the indexing scheme we have fixed is the same for each, and the 
indices which will be advanced for Ui and U2 are independent of each other in 
both expressions, Sq differs from S" only by a transposition of those two terms. 

• If Ui and U2 both have the standard basis as their eigenbases, their corresponding 
terms in Sq and S'l only have stable indices, and therefore do not advance any 
new indices. Then the term in Sq and in Si corresponding to Ui are the same, 
and similarly for U2 ; the difference between Sq and 5" is a transposition of 
those two terms, and so they are isomorphic. 

• If C/i and U2 perform the same single-qubit operation, there will be some com- 
mon tensor index Uj which is current for that qubit prior to the application of 
Ui and U2 in both Sq and S[. Then, the terms corresponding to Ui and U2 
in S'q will be C^2[i^^^] ' ™d 'he corresponding terms in Si will be 
'^ilu^] for some indices v" and w". Given that Ui and U2 perform 
the same operation, we then have Sq = S". 



Because the operations of Kq and Ki are identical after Ui and U2, the differences 
between Sq and S'l after the terms corresponding to Ui and U2 can only arise from the 
indices advanced by each term; because we have fixed the indexing scheme, the only 
difference between Sq and Si is then the (possibly trivial) transposition of those two 
terms. Thus, Sq and Si are isomorphic. ■ 



Theorem 2.1 then follows from Lemma A.l and A. 2. Using this, we may identify 



congruent circuits simply by relabelling their indices in a uniform way (e.g. as in the 



construction for stable-index expressions in Section A.2i, and determining whether 



there is a perfect matching between the two expressions defined by equality of terms. 



B Details for constructions of IWQC procedures 

B.l Construction of the gate-sets Bd^p and i^RBB 

In this section, we provide details on the construction of the IWQC procedures 3^ 
and 33, as well as a generalization of the latter operation to a many-qubit diagonal 
operation. 
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B.1.1 Two-qubit operations 

As described in ( |2.40| l, we may implement Zz in the IWQC model with the procedure 

•2, _ VSW VSW |W|V2 ^ P M (V,\\ 

^dv,w — '-V '-w '"'a '-av'-aw'^a ■ \'->-^) 



For the analysis of Section 3 it will be convenient to generahze this to a procedure 
333^j «^ acting symmetrically on many qubits: 



(B.2) 



We may compute the effect of this procedures as follows. For standard basis vectors 
, \x'^ for each 1 ^ j ^ d, ignoring the measurement result s[a] after the correction 
operations, we may obtain 



^a;iH Vxi |_|_^ ^_|_| ^^^i+'-'+^d 



. i=i 



where we define the non-negative scalars ct by 

I r^.Ti -i Vt.a 

'■it/2 



G = 

= 



X, H \-x. 



(B.3a) 



(B.3b) 
(B.3c) 



We may consider the terms of ( |B.3a| i as outer products of standard-basis vectors: we 
may verify that 



1 

V2 



(B.4a) 
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Zv\xi)vi (81 • • • (gl Zw\xd)v 



2 

7^ 



(B.4b) 



By taking outer products, and extending ( |B.3a[ ) linearly over tensor products of the 
rank-1 projectors \xi){x'i\(^ ■ ■ ■ ® \xd){x'j^\ with linear operators on an external system 
S', we then obtain 



(B.5) 



for arbitrary operators ps, acting on a system 5 which includes the qubits vj but not 
the qubit a. In the special case c? = 2, we then recover 



5}..Ups) = e— PS e-[^"«^™l/4 = Zz,^ PS Zz,^ , (B.6) 
as claimed in p.42| i. 
B.1.2 Single qubit operations 

We may describe IWQC procedures for the operators J{9) with a similar analysis as 
for the Zz operator above. For an arbitrary angle —it < 9 ^ it and qubits v and w, 
define the CPTP map by 



(B.7) 



for p a state supported on a set of qubits S which includes v but not w. We may 
characterize the effect of 5^/^ as follows. For {tp) E H2 arbitrary, we have 



(B.8) 



Ignoring the measurement result slv] after the correction operation, we may characterize 
the compound operation X^"' on a joint pure state \'^) by 



X!, 



+ (^»(-^IJ^).,»)((^l,, 



(B.9) 
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In the case — ^Z^^^ \'i>)v |+) the terms in this equation may again be ex- 

pressed as outer products of state vectors: 



(B.lOa) 



J{-'P)w/v 1-0) t> , 



(B.lOb) 



where by J{—'f)w/v we denote the linear operator Tiy Hw which maps state-vectors 
\^/j)v I — > [J{~f) linearity, we may extend this to arbitrary entangled 

states \tp)v ^ 141)3' with a system S" which does not include v. Then, for a system S 
which includes v and excludes w, we have 



(B.ll) 



B.2 Efficiency of the DKP and simpUfied RBB constructions 

Both the DKP construction and the (simplified) RBB construction, as described in 



Section 2.3.4 can be performed efficiently in the size of the input circuit. We may 
illustrate this for both procedures, as follows. In both cases, let k be the number of 
qubits which C acts upon, N be the number of one-qubit gates, and Af be the number 
of two-qubit gates: we may assume that k ^ N + M if we suppose that C acts non- 
trivially on each qubit. 



B.2.1 Obtaining a circuit in normal form 

We may easily construct a normal form in each case in a single pass of the expression 
of the circuit C. For each qubit v acted on by C, we may maintain a buffer By of 
single-qubit operations which are performed. These are initially empty for each qubit; 
we accumulate T and operations (cancelling them as appropriate) of each qubit as 
we encounter them; and upon encountering a Hadamard operation on a qubit v, we 
insert the contents of the buffer immediately before the Hadamard, resetting the buffer 
to zero after the Hadamard operation has been traversed. In both cases, we may also 
cancel AZ operations by storing in the buffers By indications of which qubits w the 
qubit V has interacted with via a AZ: initially this is set to zero, and we toggle the 
presence of each qubit in the set with each A^ operation between the two qubits. 

When we empty a buffer for a qubit v, we also produce the corresponding ^Z gates 
between v and other qubits w. When we clear a buffer By prior to a Hadamard, we 
may count the number of AZ gates between v and each qubit w, and produce a ^Z 
operation when the total is odd. For each ^Z gate between v and some other qubit 
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w, we also remove a gate in the buffer Bu, between w and v to avoid double- 
counting the /^Z gates. For the RBB construction, we may elaborate this further by 
converting operations to Zz operations by performing the substitution i — > 
Zz{T^ (g) T^)^ when we encounter such an operation, accumulating the resulting T^^ 
gates in the buffers By and B^, when we do so. 

The total amount of work performed by this routine per operation in C (counting 
both insertions and removals of terms corresponding to these operations, to or from the 
buffers) is 0(1). Using list-structures to implement the buffers By and adjacency-list 
representations for the ^Z operations, we may initialize the buffers in time 0{k). The 
amount of work required by this phase is then 0{k + N + M). 

B.2.2 Conversion to stable-index representation 

This phase may be performed simultaneously with the previous phase: with each gate 
produced, we may produce a stable-index representation of the gate, and we produce 
J{6) gates rather than H, T, and gates when we clear the buffers. In the case of the 
RBB construction, an additional amount of work may be required for the insertion of 
J {9) gates to achieve uniform wire lengths: the amount of additional work per Zz gate 
is bounded to within a constant by the number of single-qubit operations performed 
on either wire, which is at most N. The amount of work required to do this is then 
0{k + N + M) for the DKP construction, and 0{k + NM) for the RBB construction. 

B.2.3 Translation of individual operations to the 1 WQC model 

This phase can also be performed simultaneously with the previous phases: rather 
than produce a stable-index expression C", we may produce the IWQC procedure 
<!>([/[•■■] ) for each term [/[■■■] produced for the stable-index expansion. As each such 
IWQC procedure is of constant length, this requires no additional overhead asymptot- 
ically. 

B.2.4 Obtaining a IWQC procedure in normal form 

Let ^0 be the resulting IWQC procedure from the previous phase. Let n be the number 
of qubits acted on by ^Po. and and let m be the number of operations Z^, E , or 33 
in *Po- (We may easily show that n = k + N for the DKP construction and n = 
k+N + M £ 0{N + M) for theRBB construction; we have m = -h M f or the DKP 
construction and m ^ N + NM e 0{NM) for the RBB construction, essentially by 
the arguments given for the run-time for the production of the stable-index expression 
C" .) As each operation 3^. E, and 33 consists of a constant number of elementary 
operations and contains only one or two entangling operations, the procedure also 
has 0(m) elementary operations and 0{m) entanghng operations. 

Standardizing *Po • Let G be the entanglement graph associated to • In the process 
of standardizing *Po. for each qubit v which is part of a Zv/u operation for some u, we 
must conamute a single operation X^'"' past entanghng operations Eyyj for all w adjacent 
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to V in the graph G. Rather than expHcitly commuting operations, we may standardize 
^0 by preparing a new IWQC procedure *Pi as follows: 

1. We initialize *Pi by assigning to it the procedure consisting of all preparation 
maps N„ in ^q, followed by all entangling maps E^w (Note that the sets of 
prepared qubits and the graph G of entangling operators can be obtained with no 
additional asymptotic cost, during the phase of translation to the IWQC model 
in which *Po is produced.) 

2. Note that every deprecated index Vj in the stable-index expansion of C" has a 
corresponding advanced index Wj+i- As the procedures are the only op- 
erations in the DKP construction which allocate or discard qubits, each qubit u 
which is discarded in has a corresponding qubit v which is allocated, and 
this correspondence arises precisely from the presence of an operation in 
*Po. If / is the index successor function of C", we then have v = f{u) for 
such qubits u and v. For the RBB construction, we similarly have a gate 

for qubits with f{u) — v; but we also have qubits a which are allocated and 
discarded in the operations 33, which lie outside of the domain of /. We may 
then extend / to a function / such that /(a) — a for such mediating qubits, and 
/('^i) = fi'^^j) — "^j+i foi" ''^j ^ dom(/). (Such a function / corresponds to an 
extended successor function as defined in |Definition 3.1| ) For the DKP construc- 
tion, we will use / = /. Note that we may again construct / at no additional 
cost during the phase in which is produced from the stable-index expression 
C". 

We may use the function / to append the measurements of *Po to the left-hand 
side of *Pi in order, with the correct sign- and bit-dependencies included. We 
may make a single pass of and in doing so maintain a list of correction 
operations for each of the qubits in *Po> as follows. Initially, no qubit has any 
corrections. For any measurement M* encountered, we perform the following: 

(/) If u does not have any accumulated corrections, we append to the 
left-hand side of *Pi without any dependencies; 

(//) If u has accumulated corrections and Z^, we append S^j^j M^''' to the 
left-hand side of (If 7 = 0, we may omit the shift operation.) 

(Hi) We add anX^'"' correction to v = f{u), provided that v ^ u, corresponding 
to the X correction arising from an operation (Note that the condition 

u 7^ f{u) is satisfied for all u in the DKP construction; it is satisfied if and 
only if u g dom(/) in the RBB construction.) 

(iV) We add Z^"' corrections for every w adjacent to u = /(u) in the graph G, 
provided w ^ u. (In the DKP construction, and for qubits u e dom(/) 
in the RBB construction, these corrections arise from commuting the X^'"' 
operation past the entangling maps E^^ for w ^ u; for qubits a ^ dom(/) 
in the latter construction, this corresponds to the Z^"' corrections in the 33 
procedure, in which case the condition w ^ u will be satisfied for all w.) 
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We may simplify the expressions for any accumulated corrections using the re- 
lations X^XT^ = X^^'"' and Zgzt^ = 2^^+='"'. 

The amount of work performed for each measurement M* in *Po is at most 
0( deg[/(u)]), where deg(w) is the number of neighbors of v in G. The total 
number of operations performed in this stage is then 0(J2 deg{v)) = 
0{\EiG)\)=Oim). 

3. After we have read all of the measurements, we append the corrections for each 
unmeasured qubit (computed in the previous step) to to the end of 

The resulting procedure *Pi is then the same that would result from standardizing *Po> 
and can be obtained in time 0{m). 



Normalizing Obtaining a normal form *p from *Pi can be done by checking 

for Pauli simplifications for each measurement, and then performing signal shifting by 
accumulating shift operators to the left (without commuting shift operators past each 
other). We may check for and apply PauU simplifications as a part of the standardiza- 
tion process above, with no additional overhead asymptotically. We may also perform 
signal shifting simultaneously with the standardization process, by immediately com- 



muting the shift operation S^^^^ 



past any corrections Xi,'"' 
been introduced by the measurement on a given qubit u: 
cency relation in G, we have 
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After performing this commutation, there are no classically controlled operations to be 
inserted to the left of the shift operator which will depend on s[m]; we may then com- 
mute it past all additional operations and remove it from the expression. Then, rather 
than inserting shift operations and the corrections described, we may accumulate the 
corrections described on the right-hand side after the measurement on u. The amount 
of work required for each measurement in this alternative procedure is then dominated 
by the work required to update the corrections, which is 0{ndeg\f(u)\); summed over 
all vertices u, obtaining a normal form may then be done in time 0{nm). 



B.2.5 Total run-time complexity 

Comparing the complexity of the phases above (and recalling the relationship between 



n and m with N , M, and k at the beginning of Section B.2.4 the run-time of both the 
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DKP and RBB constructions are then dominated by the phase of obtaining a IWQC 
procedure in normal form. This requires time 0{nm) for both constructions, which 
evaluates to 0{{k + N)M) for the DKP construction and 0{{N + M)NM) for the 
RBB construction. This bound on the run-time coincides with that of the semantic map 
S described in lSection 3.51 
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